Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0170CT31.3% 329.1 / 163.8 213intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (67/67);  new base T (60/1);  total (144/69)
Fisher's exact test for biased strand distribution p-value = 2.95e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 215 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGT  >  NC_000913/3705918‑3706114
                                                                                                    |                                                                                                         
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                       >  2:465919/1‑100 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                       >  2:3006490/1‑100 (MQ=255)
     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                                                      <  1:1359974/100‑1 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                                                      <  2:1137404/98‑1 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGTAAATTCAATACATTGCA‑g                                                                                                   <  1:661341/100‑1 (MQ=255)
        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                                                                  <  2:2580960/100‑1 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                                                                >  1:792665/1‑100 (MQ=255)
          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At                                                                                               <  2:816131/100‑1 (MQ=255)
            ntaCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt                                                                                            >  2:1656929/2‑100 (MQ=255)
            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt                                                                                            >  1:582416/1‑100 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt‑‑‑t                                                                                        <  1:2824942/100‑1 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑Tc                                                                                       <  2:476873/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAgg                                                                                <  1:2576750/100‑1 (MQ=255)
                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGc                                                                               >  2:1267205/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGa                                                                           >  2:1433147/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGa                                                                           >  1:23376/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGa                                                                           <  1:575699/100‑1 (MQ=255)
                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGAt                                                                          >  1:715506/1‑100 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATa                                                                         <  2:358458/100‑1 (MQ=255)
                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATaa                                                                        >  1:919992/1‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGc                                                                      >  2:1145561/1‑100 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCg                                                                     <  2:1483982/100‑1 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCg                                                                     >  2:1636103/1‑100 (MQ=255)
                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGa                                                                   >  1:769878/1‑100 (MQ=255)
                                  ncTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                 >  1:2018889/2‑100 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                 >  1:804200/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                <  2:311916/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                <  2:578073/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                <  1:1122282/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                >  1:1461731/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                >  1:2441741/1‑100 (MQ=255)
                                    tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CA‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAgcgc                                                               >  1:1794836/10‑100 (MQ=18)
                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTCCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcgc                                                               <  1:1510908/100‑1 (MQ=255)
                                     ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCa                                                              >  1:602934/9‑100 (MQ=18)
                                      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑CA‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                             >  1:2353187/8‑100 (MQ=17)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                          <  2:1622741/100‑1 (MQ=255)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                         >  2:2779604/4‑100 (MQ=21)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑TG‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:1399973/4‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:1546451/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  2:1255125/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CA‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:856168/3‑100 (MQ=18)
                                            ggtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑GA‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       <  1:525882/90‑1 (MQ=21)
                                            ggctggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑CA‑GC‑TA‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:189451/11‑100 (MQ=25)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2344358/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:2546412/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:366627/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                        >  2:494779/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑AT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                       >  1:944018/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑CC‑‑‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2946719/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:2050703/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑AC‑‑‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                       >  1:493374/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                      >  1:1178203/2‑100 (MQ=21)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                      <  2:2670219/100‑1 (MQ=255)
                                             gCCGGATGCGGCGATAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TC‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGACa                                                       >  2:667640/1‑99 (MQ=255)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑AT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                      >  1:3034754/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                      >  2:2927051/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                     >  2:1281867/1‑100 (MQ=21)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                     >  1:2206618/9‑98 (MQ=18)
                                              ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                     >  1:2484979/9‑100 (MQ=18)
                                              caggataaggcgTTCACGCCTTATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑TG‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                     >  1:2759267/9‑100 (MQ=14)
                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                    >  1:2384252/8‑97 (MQ=21)
                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                                                   >  2:770288/8‑99 (MQ=21)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:393530/7‑100 (MQ=17)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:843012/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:552389/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑CT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                   >  1:1650576/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑TT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                    >  2:2913916/7‑96 (MQ=17)
                                                ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                   >  1:1867700/7‑100 (MQ=18)
                                                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑TG‑GG‑CC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttt                                                  >  1:352443/6‑100 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑AT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatctt                                                  >  2:1867264/6‑95 (MQ=25)
                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑CT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt                                                  >  2:1340591/6‑95 (MQ=18)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAACACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGTCATTtttg                                                 <  2:603550/100‑1 (MQ=255)
                                                   taaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                                                >  1:1130765/4‑100 (MQ=18)
                                                   taaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑TA‑‑‑T‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctttac                                                >  1:1623825/4‑93 (MQ=21)
                                                      ggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TA‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAATtttgcgtt                                             >  1:2113271/1‑97 (MQ=17)
                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                             >  1:983734/1‑100 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                        >  1:2740386/1‑100 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                        <  2:226733/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtac                                     <  2:1448560/100‑3 (MQ=255)
                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                     >  1:2969477/1‑100 (MQ=255)
                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTttttcttttttttttag                                    >  2:1062233/1‑91 (MQ=18)
                                                                      ccGGCCTAT‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  1:2836218/1‑100 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  1:1701931/1‑100 (MQ=255)
                                                                        ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                           <  2:2181151/100‑1 (MQ=255)
                                                                           cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                        >  2:2090293/1‑100 (MQ=255)
                                                                           cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                        >  2:43908/1‑100 (MQ=255)
                                                                             aC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                      >  1:15807/1‑100 (MQ=255)
                                                                              c‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                     <  1:1325530/100‑1 (MQ=255)
                                                                                gaaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                    <  2:2437148/99‑1 (MQ=255)
                                                                                aaaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                    <  2:191995/100‑1 (MQ=255)
                                                                                aaaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                    <  1:1814020/100‑1 (MQ=255)
                                                                                 aaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                   <  2:2396108/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑CT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTGTa                    <  2:855359/98‑1 (MQ=255)
                                                                                          aaaTTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt          <  1:2852497/100‑1 (MQ=255)
                                                                                             ttCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga       <  2:251794/100‑1 (MQ=255)
                                                                                              tCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag      <  2:1753750/100‑1 (MQ=255)
                                                                                              tCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag      <  1:629218/100‑1 (MQ=255)
                                                                                               cAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga     >  1:1987479/1‑100 (MQ=255)
                                                                                               cAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga     <  1:1023507/100‑1 (MQ=255)
                                                                                                  tACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGt  >  1:686296/1‑100 (MQ=255)
                                                                                                    |                                                                                                         
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGT  >  NC_000913/3705918‑3706114

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