Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,208,201 G→A 6.7% V60M (GTG→ATG)  ttdT → L‑tartrate/succinate antiporte

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,208,2010GA 6.7% 284.9 / 9.2 89V60M (GTG→ATG) ttdTL‑tartrate/succinate antiporte
Reads supporting (aligned to +/- strand):  ref base G (44/39);  new base A (2/4);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 103 total aligned reads displayed.

TGCGCCGCTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGC  >  NC_000913/3208056‑3208349
                                                                                                                                                 |                                                                                                                                                    
tGCGCCGCTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   <  3:173750/149‑1 (MQ=255)
  cgccgcTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGt                                                                                                                                                 <  3:68009/149‑1 (MQ=255)
      gcTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTatat                                                                                                                                             <  3:1299957/149‑1 (MQ=255)
             gTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                      <  4:299125/149‑1 (MQ=255)
             gTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCcaccat                                                                                                                                      <  2:1273951/149‑1 (MQ=255)
              tCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCatcatc                                                                                                                                     <  1:1004646/149‑1 (MQ=255)
                aTCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGc                                                                                                                                   >  2:1103938/1‑149 (MQ=255)
                aTCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGc                                                                                                                                   >  1:34111/1‑149 (MQ=255)
                   gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   <  3:648538/130‑1 (MQ=255)
                   gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   >  2:1299610/1‑130 (MQ=255)
                   gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   >  4:648538/1‑130 (MQ=255)
                   gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   <  1:1299610/130‑1 (MQ=255)
                    ccATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTggg                                                                                                                                                  <  2:650699/130‑1 (MQ=255)
                    ccATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTggg                                                                                                                                                  >  1:650699/1‑130 (MQ=255)
                    ccATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                               <  2:34111/149‑1 (MQ=255)
                    ccATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                               >  1:1100734/1‑149 (MQ=255)
                          ttGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATGCCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCatcatc                                                                                                                                     <  1:591638/137‑1 (MQ=255)
                          ttGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGACCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCatcatc                                                                                                                                     >  2:591638/1‑137 (MQ=255)
                          ttGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcat                                                                                                                                         <  2:892443/133‑1 (MQ=255)
                          ttGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcat                                                                                                                                         >  1:892443/1‑133 (MQ=255)
                             ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   <  2:1236503/120‑1 (MQ=255)
                             ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   >  1:1236503/1‑120 (MQ=255)
                             ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTccc                                                                                                                      >  4:1412395/1‑149 (MQ=255)
                              tctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCt                                                                                                                     >  1:1076202/1‑149 (MQ=255)
                                   ttCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   <  4:5341/114‑1 (MQ=255)
                                   ttCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   >  3:5341/1‑114 (MQ=255)
                                   ttCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                               <  1:1434249/134‑1 (MQ=255)
                                   ttCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                               >  2:1434249/1‑134 (MQ=255)
                                   ttCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgctg                                                                                                                >  2:934327/1‑149 (MQ=255)
                                     ccGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgctgct                                                                                                              >  1:442673/1‑149 (MQ=255)
                                      cGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGta                                                                                                                                                >  1:150777/1‑114 (MQ=255)
                                      cGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGta                                                                                                                                                <  4:1411652/114‑1 (MQ=255)
                                      cGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGta                                                                                                                                                <  2:150777/114‑1 (MQ=255)
                                      cGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGta                                                                                                                                                >  3:1411652/1‑114 (MQ=255)
                                       ggTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGt                                                                                                            >  3:518566/1‑149 (MQ=255)
                                        gTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTGCTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGtt                                                                                                           >  4:1193680/1‑149 (MQ=255)
                                        gTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGta                                                                                                                                                <  3:1173820/112‑1 (MQ=255)
                                        gTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGta                                                                                                                                                >  4:1173820/1‑112 (MQ=255)
                                          tCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCa                                                                                                         <  2:1100734/149‑1 (MQ=255)
                                              gCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGccc                                                                                                     >  1:1015018/1‑149 (MQ=255)
                                                     tGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagcag                                                                                              <  3:1412395/149‑1 (MQ=255)
                                                      ggAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagcagc                                                                                             <  2:442673/149‑1 (MQ=255)
                                                       gagaACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   >  4:482266/1‑94 (MQ=255)
                                                       gagaACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGAtggtgg                                                                                                                                                   <  3:482266/94‑1 (MQ=255)
                                                            ccATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgctcgctc                                                                                       >  2:713296/1‑149 (MQ=255)
                                                             cATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCaccatc                                                                                                                                     >  1:808212/1‑102 (MQ=38)
                                                             cATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCaccatc                                                                                                                                     <  2:808212/102‑1 (MQ=39)
                                                             cATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCaccatc                                                                                                                                     <  2:464696/102‑1 (MQ=39)
                                                             cATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCaccatc                                                                                                                                     >  1:464696/1‑102 (MQ=39)
                                                               tACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGc                                                                                    >  3:118952/1‑149 (MQ=255)
                                                                    ggCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGc                                                                               <  1:934327/149‑1 (MQ=255)
                                                                      ctctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCtt                                                                             <  2:1015018/149‑1 (MQ=255)
                                                                       tctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgct                                                                                                                 >  2:1101847/1‑112 (MQ=255)
                                                                       tctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgct                                                                                                                 <  1:1101847/112‑1 (MQ=255)
                                                                          aCTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                         <  1:1103938/149‑1 (MQ=255)
                                                                            tttGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGc                                                                                                                                   <  2:810962/89‑1 (MQ=255)
                                                                            tttGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGc                                                                                                                                   >  1:810962/1‑89 (MQ=255)
                                                                                 cGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCAc                                                                    <  2:1087579/147‑1 (MQ=255)
                                                                                 cGTTTTTACTGCCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCAc                                                                    >  1:1087579/1‑147 (MQ=255)
                                                                                  gTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                         >  2:66880/1‑141 (MQ=255)
                                                                                  gTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                         <  1:66880/141‑1 (MQ=255)
                                                                                   tttttACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCa                                                                >  3:17099/1‑149 (MQ=255)
                                                                                             cGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCctct                                                      <  4:118952/149‑1 (MQ=255)
                                                                                              gTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagcagc                                                                                             <  2:1194545/109‑1 (MQ=255)
                                                                                              gTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagcagc                                                                                             >  1:1194545/1‑109 (MQ=255)
                                                                                                 aTCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTg                                                  >  4:525720/1‑149 (MQ=255)
                                                                                                 aTCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTg                                                  >  3:367044/1‑149 (MQ=255)
                                                                                                 aTCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTg                                                  >  3:1437601/1‑149 (MQ=255)
                                                                                                    gTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGc                                               >  1:354991/1‑149 (MQ=255)
                                                                                                       ggACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGt                                            >  3:132778/1‑149 (MQ=255)
                                                                                                        gACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGa                                                                                                  >  3:985299/1‑94 (MQ=255)
                                                                                                        gACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGtt                                           <  4:17099/149‑1 (MQ=255)
                                                                                                          cTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTc                                         >  2:6941/1‑149 (MQ=255)
                                                                                                           tGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTcc                                        <  4:518566/149‑1 (MQ=255)
                                                                                                            gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgct                                                                                                                 <  1:1053670/75‑1 (MQ=255)
                                                                                                            gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGgct                                                                                                                 >  2:1053670/1‑75 (MQ=255)
                                                                                                            gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCg                                       <  1:713296/149‑1 (MQ=255)
                                                                                                             aTCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCAc                                                                    <  4:1072154/119‑1 (MQ=255)
                                                                                                             aTCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCAc                                                                    >  3:1072154/1‑119 (MQ=255)
                                                                                                                   gAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTc                                <  3:1193680/149‑1 (MQ=255)
                                                                                                                       ccGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                     <  1:643721/108‑1 (MQ=255)
                                                                                                                       ccGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                     >  2:643721/1‑108 (MQ=255)
                                                                                                                          tGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCg                         <  4:132778/149‑1 (MQ=255)
                                                                                                                            ccGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATActc                                                                                                                            >  2:1251850/1‑48 (MQ=255)
                                                                                                                            ccGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATActc                                                                                                                            <  1:1251850/48‑1 (MQ=255)
                                                                                                                                gTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTATTCAGCCGGGAGCAGCTCGCTCAGCCAGTCTTTAAATTCGCTGCAAAATGCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATc                   >  3:911036/1‑149 (MQ=255)
                                                                                                                                     gtggtgGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCt              <  4:367044/149‑1 (MQ=255)
                                                                                                                                     gtggtgGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCt              >  1:526093/1‑149 (MQ=255)
                                                                                                                                      tggtggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATc                                                           <  1:1231256/103‑1 (MQ=255)
                                                                                                                                      tggtggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATc                                                           >  2:1231256/1‑103 (MQ=255)
                                                                                                                                      tggtggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGt                                     >  1:585821/1‑125 (MQ=255)
                                                                                                                                      tggtggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGt                                     <  2:585821/125‑1 (MQ=255)
                                                                                                                                       ggtggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTgg                <  4:939332/145‑1 (MQ=255)
                                                                                                                                       ggtggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTgg                >  3:939332/1‑145 (MQ=255)
                                                                                                                                             gatgatGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                         <  4:1228460/82‑1 (MQ=38)
                                                                                                                                             gatgatGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                         >  3:1228460/1‑82 (MQ=38)
                                                                                                                                             gATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGcc                                                                                                      >  2:589984/1‑53 (MQ=255)
                                                                                                                                             gATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGcc                                                                                                      <  1:589984/53‑1 (MQ=255)
                                                                                                                                             gATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg      >  3:271838/1‑149 (MQ=255)
                                                                                                                                                 gtggGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcgccgc  >  1:618964/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGCGCCGCTGGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGC  >  NC_000913/3208056‑3208349

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