Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,087,214 C→A P107T (CCA→ACA)  SAUSA300_RS10530 → chemotaxis inhibitory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,087,2140CA99.8% 3689.4 / NA 1261P107T (CCA→ACA) SAUSA300_RS10530chemotaxis inhibitory protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (962/296);  total (963/298)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1326 total aligned reads displayed.

AAATTATTTAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACACCACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGTAATTA  >  NC_007793/2087117‑2087310
                                                                                                 |                                                                                                 
 aaTTATTTAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATAcaacac                                                                                                >  1:1577762/1‑100 (MQ=255)
       ttAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGccg                                                                                          <  1:452334/100‑3 (MQ=255)
       ttAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTg                                                                                          <  1:1545272/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGa                                                                                         <  1:250575/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGa                                                                                         <  1:110942/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGa                                                                                         <  1:3533154/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGa                                                                                         <  1:4188398/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGa                                                                                         <  1:1954102/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGa                                                                                         <  1:1872148/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTACTGa                                                                                         <  1:1471682/100‑1 (MQ=255)
        tAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTCGCTGa                                                                                         <  1:2063899/100‑1 (MQ=255)
         aaagaaagGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGag                                                                                        <  1:4243507/100‑2 (MQ=255)
             caagGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATGCAGTATATCTGAATACAACACTTGCTGAAGATa                                                                                    >  1:694308/2‑100 (MQ=255)
             caagGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATa                                                                                    >  1:5397423/2‑100 (MQ=255)
              aagGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGAttg                                                                                   >  1:590543/1‑98 (MQ=255)
              aagGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGACAg                                                                                   <  1:2812959/100‑1 (MQ=255)
              aagGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTAGCTGAAGAtcg                                                                                   >  1:2111882/1‑98 (MQ=255)
               ctgCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑a                                                                                 >  1:1275469/3‑100 (MQ=255)
               agGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTACTGAAGAtcgg                                                                                  >  1:1006651/1‑97 (MQ=255)
                ggCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑aa                                                                                >  1:1911594/1‑100 (MQ=255)
                ggCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑aa                                                                                >  1:4905009/1‑100 (MQ=255)
                         aaTTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTg                                                                       <  1:380658/100‑1 (MQ=255)
                         aaTTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTg                                                                       <  1:2152711/100‑1 (MQ=255)
                         aaTTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGCTg                                                                       <  1:5692843/100‑1 (MQ=255)
                           ctCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAACGTTGAg                                                                     >  1:525759/2‑100 (MQ=255)
                             cTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAGAAAAAAATGTTGAGt                                                                    >  1:3654840/1‑100 (MQ=255)
                             cTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGtt                                                                   >  1:3239700/1‑100 (MQ=255)
                             cTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGtt                                                                   >  1:2332130/1‑100 (MQ=255)
                              tGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTa                                                                  >  1:866251/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAGAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:4476628/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGCTAG‑TAAAAATGTTGAGTTAc                                                                 >  1:2714403/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAgatcg                                                                 >  1:3047775/1‑96 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:2464364/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:4059947/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:3984769/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:4957690/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:1665436/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:2905652/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:3124862/1‑100 (MQ=255)
                               gCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAc                                                                 >  1:5309663/1‑100 (MQ=255)
                                cTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTAcc                                                                <  1:3705956/100‑2 (MQ=255)
                                  cAATTTGTAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGATGATAG‑AAAAAATGTTGAGTTACTAg                                                              >  1:3477038/1‑100 (MQ=255)
                                  cAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAg                                                              >  1:2009341/1‑100 (MQ=255)
                                  cAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAg                                                              >  1:4821124/1‑100 (MQ=255)
                                  cAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAg                                                              >  1:5009642/1‑100 (MQ=255)
                                  cAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAg                                                              >  1:3007356/1‑100 (MQ=255)
                                   aaTTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAgg                                                             <  1:4383198/100‑1 (MQ=255)
                                            aaaTGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAgatcggaagag                                                    >  1:1369905/1‑90 (MQ=255)
                                              aTGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGtata                                                  >  1:5746920/1‑100 (MQ=255)
                                               tGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATaa                                                 >  1:1182515/1‑100 (MQ=255)
                                                     gcttAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATa                                           >  1:1803350/3‑100 (MQ=255)
                                                     gcttAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATa                                           >  1:5447064/3‑100 (MQ=255)
                                                       ttAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATAct                                         >  1:64960/1‑100 (MQ=255)
                                                           cTGAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGCAAAATGTATAAAACATActtctt                                     >  1:1427356/1‑100 (MQ=255)
                                                             gAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTa                                   >  1:1229971/1‑100 (MQ=255)
                                                             gAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTa                                   >  1:5649648/1‑100 (MQ=255)
                                                             gAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTa                                   >  1:5570272/1‑100 (MQ=255)
                                                             gAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTa                                   >  1:2552993/1‑100 (MQ=255)
                                                             gAAAATAAAGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGCTACTAGGTAAAATGTATAAAACATACTTCTTTa                                   >  1:4418284/1‑100 (MQ=255)
                                                                   aaaGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAg                             >  1:5073448/1‑100 (MQ=255)
                                                                     aGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAGAAAAGga                           >  1:4724681/1‑100 (MQ=255)
                                                                     aGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGga                           >  1:2862933/1‑100 (MQ=255)
                                                                     aGGTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGga                           >  1:2192914/1‑100 (MQ=255)
                                                                      ggTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAATGgag                          <  1:3188607/100‑1 (MQ=255)
                                                                      ggTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGgag                          <  1:3749562/100‑1 (MQ=255)
                                                                      ggTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGgag                          <  1:4578265/100‑1 (MQ=255)
                                                                      ggTTACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGCATAAAACATACTTCTTTAAAAgaggag                          <  1:2415165/100‑1 (MQ=255)
                                                                         tACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAgggagag                       >  1:228661/1‑95 (MQ=255)
                                                                         tACTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGt                       >  1:1071114/1‑100 (MQ=255)
                                                                          aCTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGagagat                      >  1:4304104/1‑98 (MQ=255)
                                                                          aCTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGagagat                      >  1:1144098/1‑98 (MQ=255)
                                                                          aCTATACAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGagagat                      >  1:5344787/1‑98 (MQ=255)
                                                                             ataCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAGAGGagagatcgg                   >  1:131212/1‑95 (MQ=255)
                                                                             ataCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGagagatcgg                   >  1:954759/1‑95 (MQ=255)
                                                                             ataCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGagagatcgg                   >  1:3938847/1‑95 (MQ=255)
                                                                             ataCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTaa                   >  1:3397309/1‑100 (MQ=255)
                                                                             ataCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTaa                   >  1:5501849/1‑100 (MQ=255)
                                                                             ataCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTT‑AAAAAGGAGGGGCTaaa                  >  1:5626734/1‑100 (MQ=255)
                                                                              taCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTaaa                  >  1:2758813/1‑100 (MQ=255)
                                                                               aCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGagagatcggaa                 >  1:2230656/1‑93 (MQ=255)
                                                                               aCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAat                 >  1:2649473/1‑100 (MQ=255)
                                                                               aCAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAat                 >  1:3076057/1‑100 (MQ=255)
                                                                                cAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAatc                >  1:475449/1‑100 (MQ=255)
                                                                                cAGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAatc                >  1:5162471/1‑100 (MQ=255)
                                                                                 aGTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGagagatcggaaga               >  1:1629741/1‑91 (MQ=255)
                                                                                  gTATATCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAatcat              >  1:3324645/1‑100 (MQ=255)
                                                                                    gtatCTGAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCt            >  1:3866407/2‑100 (MQ=255)
                                                                                         cgAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGt       >  1:2954021/2‑100 (MQ=255)
                                                                                          gAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGta      >  1:4857069/1‑100 (MQ=255)
                                                                                          gAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGta      >  1:336149/1‑100 (MQ=255)
                                                                                          gAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGta      >  1:5217560/1‑100 (MQ=255)
                                                                                          gAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGta      >  1:1727683/1‑100 (MQ=255)
                                                                                          gAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGta      >  1:760010/1‑100 (MQ=255)
                                                                                          gAATACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGta      >  1:907904/1‑100 (MQ=255)
                                                                                             tACAACACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGtaatt   >  1:3577983/1‑100 (MQ=255)
                                                                                              acaacaCTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGGGAGTCTAAATCATCTTATGtaatta  >  1:798848/1‑100 (MQ=255)
                                                                                              acaacaCTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGtaatta  >  1:2590731/1‑100 (MQ=255)
                                                                                              acaacaCTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGtaatta  >  1:5131302/1‑100 (MQ=255)
                                                                                              acaacaCTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGtaatta  >  1:4138516/1‑100 (MQ=255)
                                                                                              acaacaCTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGtaatta  >  1:4670429/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
AAATTATTTAAAGAAAGGCACAAAAAATTCTGCTCAATTTGAAAAAATGGTTATTTTAACTGAAAATAAAGGTTACTATACAGTATATCTGAATACACCACTTGCTGAAGATAG‑AAAAAATGTTGAGTTACTAGGTAAAATGTATAAAACATACTTCTTTAAAAAAGGAGAGTCTAAATCATCTTATGTAATTA  >  NC_007793/2087117‑2087310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑