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.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGT  >  NC_000913/3705918‑3706114
                                                                                                    |                                                                                                         
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                                                       >  2:465919/1‑100 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                                                       >  2:3006490/1‑100 (MQ=255)
     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                                                      <  1:1359974/100‑1 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGc                                                                                                      <  2:1137404/98‑1 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCAg                                                                                                   <  1:661341/100‑1 (MQ=255)
        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                                                                  <  2:2580960/100‑1 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                                                                >  1:792665/1‑100 (MQ=255)
          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt                                                                                               <  2:816131/100‑1 (MQ=255)
            ntaCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGCAGAAttt                                                                                            >  2:1656929/2‑100 (MQ=255)
            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAttt                                                                                            >  1:582416/1‑100 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAttt‑‑‑t                                                                                        <  1:2824942/100‑1 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑Tc                                                                                       <  2:476873/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAgg                                                                                <  1:2576750/100‑1 (MQ=255)
                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGc                                                                               >  2:1267205/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGa                                                                           >  2:1433147/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGa                                                                           >  1:23376/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGa                                                                           <  1:575699/100‑1 (MQ=255)
                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGAt                                                                          >  1:715506/1‑100 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATa                                                                         <  2:358458/100‑1 (MQ=255)
                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATaa                                                                        >  1:919992/1‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGc                                                                      >  2:1145561/1‑100 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCg                                                                     <  2:1483982/100‑1 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCg                                                                     >  2:1636103/1‑100 (MQ=255)
                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGa                                                                   >  1:769878/1‑100 (MQ=255)
                                  ncTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                 >  1:2018889/2‑100 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                 >  1:804200/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                <  2:311916/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                <  2:578073/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                <  1:1122282/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                >  1:1461731/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                >  1:2441741/1‑100 (MQ=255)
                                    tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAgcgc                                                               >  1:1794836/10‑100 (MQ=18)
                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTCCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcgc                                                               <  1:1510908/100‑1 (MQ=255)
                                     ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCa                                                              >  1:602934/9‑100 (MQ=18)
                                      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                             >  1:2353187/8‑100 (MQ=17)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                          <  2:1622741/100‑1 (MQ=255)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                         >  2:2779604/4‑100 (MQ=21)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATG‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:1399973/4‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:1546451/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  2:1255125/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:856168/3‑100 (MQ=18)
                                            ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       <  1:525882/90‑1 (MQ=21)
                                            ggctggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGCACAGCTA‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:189451/11‑100 (MQ=25)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2344358/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:2546412/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:366627/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                        >  2:494779/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                       >  1:944018/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCC‑‑‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2946719/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:2050703/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAAC‑‑‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                       >  1:493374/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                      >  1:1178203/2‑100 (MQ=21)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                      <  2:2670219/100‑1 (MQ=255)
                                             gCCGGATGCGGCGATAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTC‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGACa                                                       >  2:667640/1‑99 (MQ=255)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                      >  1:3034754/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                      >  2:2927051/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                     >  2:1281867/1‑100 (MQ=21)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCT‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                     >  1:2206618/9‑98 (MQ=18)
                                              ccggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                     >  1:2484979/9‑100 (MQ=18)
                                              caggataaggcgTTCACGCCTTATCCTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATG‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                     >  1:2759267/9‑100 (MQ=14)
                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCC‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                    >  1:2384252/8‑97 (MQ=21)
                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                                                   >  2:770288/8‑99 (MQ=21)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:393530/7‑100 (MQ=17)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:843012/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCT‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:552389/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                   >  1:1650576/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                    >  2:2913916/7‑96 (MQ=17)
                                                ggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                   >  1:1867700/7‑100 (MQ=18)
                                                 gataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttt                                                  >  1:352443/6‑100 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTAT‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatctt                                                  >  2:1867264/6‑95 (MQ=25)
                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt                                                  >  2:1340591/6‑95 (MQ=18)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGTCATTtttg                                                 <  2:603550/100‑1 (MQ=255)
                                                   taaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCC‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                                                >  1:1130765/4‑100 (MQ=18)
                                                   taaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTA‑‑‑T‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctttac                                                >  1:1623825/4‑93 (MQ=21)
                                                      ggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTA‑‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAATtttgcgtt                                             >  1:2113271/1‑97 (MQ=17)
                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                             >  1:983734/1‑100 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                        >  1:2740386/1‑100 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                        <  2:226733/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtac                                     <  2:1448560/100‑3 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                     >  1:2969477/1‑100 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTttttcttttttttttag                                    >  2:1062233/1‑91 (MQ=18)
                                                                      ccGGCCTATAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  1:2836218/1‑100 (MQ=255)
                                                                      ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  1:1701931/1‑100 (MQ=255)
                                                                        ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                           <  2:2181151/100‑1 (MQ=255)
                                                                           cTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                        >  2:2090293/1‑100 (MQ=255)
                                                                           cTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                        >  2:43908/1‑100 (MQ=255)
                                                                             aCAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                      >  1:15807/1‑100 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                     <  1:1325530/100‑1 (MQ=255)
                                                                                gaaaTCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                    <  2:2437148/99‑1 (MQ=255)
                                                                                aaaaTCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                    <  2:191995/100‑1 (MQ=255)
                                                                                aaaaTCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                    <  1:1814020/100‑1 (MQ=255)
                                                                                 aaaTCGTGCAAATTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                   <  2:2396108/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGCAGAATCT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTGTa                    <  2:855359/98‑1 (MQ=255)
                                                                                          aaaTTCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt          <  1:2852497/100‑1 (MQ=255)
                                                                                             ttCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga       <  2:251794/100‑1 (MQ=255)
                                                                                              tCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag      <  2:1753750/100‑1 (MQ=255)
                                                                                              tCAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag      <  1:629218/100‑1 (MQ=255)
                                                                                               cAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga     >  1:1987479/1‑100 (MQ=255)
                                                                                               cAATACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga     <  1:1023507/100‑1 (MQ=255)
                                                                                                  tACATTGCAGAATTT‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGt  >  1:686296/1‑100 (MQ=255)
                                                                                                    |                                                                                                         
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTT‑‑‑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: