Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC43.9% 217.7 / inf 278intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (80/76);  new base C (114/8);  total (194/84)
Fisher's exact test for biased strand distribution p-value = 2.24e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-02

Only 100 of 281 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083
                                                                                                  |                                                                                                            
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                       <  2:2753336/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGca                                                                                                      <  2:2484250/100‑2 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                    >  2:2183530/1‑100 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                                    <  1:1215384/100‑1 (MQ=37)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                                   <  1:1010223/100‑2 (MQ=25)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                                 >  2:460815/1‑100 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                                 >  2:2027214/1‑100 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                                <  2:558666/100‑1 (MQ=255)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                              >  2:2927379/1‑100 (MQ=255)
                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                            >  2:2432225/1‑100 (MQ=255)
                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATAGTTGACAAAGTGCGCTTTGGTTAAGCCAGATGAGGGGTAAACTCCTTATCCGGCCTACAAAATc                                                                                            >  1:2186233/1‑100 (MQ=21)
                  aCGGAGAGGCACAGCGCCACCGAGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                            <  1:2413519/98‑1 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                          >  1:1478884/1‑100 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTCTTGACAAAGTGCGATTTGTTTATGCCGGATGCAGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                          <  1:2680874/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                        <  1:708527/100‑1 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                       >  1:2811676/1‑100 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                     <  1:1108216/100‑1 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                    <  2:2584305/100‑1 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACTCCTTATCCGGCCTACAAAAGCGTGCaaa                                                                                     >  2:1026532/1‑98 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                               <  2:3049609/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                              <  1:2702218/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                              >  2:1979468/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                          >  2:417166/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                          <  2:1550325/100‑1 (MQ=255)
                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                        <  2:1082357/100‑1 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t                                                              <  1:3343268/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t‑‑‑tt                                                         <  2:212359/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t‑‑‑tt                                                         >  2:2303659/1‑100 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGGGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑Gt                                                    >  1:1851281/1‑100 (MQ=255)
                                                 ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAg                                                  >  2:2701999/8‑100 (MQ=18)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGc                                                <  2:1707828/100‑1 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCt                                              >  2:2837289/1‑100 (MQ=255)
                                                      tgACAAAGTGCTGTTTGGTTCTGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACAAAATCGTGCAGATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTg                                             <  1:1577931/100‑1 (MQ=17)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTg                                             >  1:1367345/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGa                                            >  2:2217062/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGAt                                           <  1:377126/100‑1 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGc                                       <  2:3334909/100‑1 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     <  2:1283972/100‑1 (MQ=255)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  2:797891/1‑100 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAA‑‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  2:3229828/7‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑C‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  2:2256683/7‑100 (MQ=18)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCTGCCTACAAAATCGTGCAAATTCAATACATGGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGAGAAGCGGAGCGCATCAg                           >  1:2027745/1‑100 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           >  1:2256146/1‑100 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCAT‑‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:2786412/5‑100 (MQ=21)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCT‑‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:3063648/5‑100 (MQ=18)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          >  2:1173356/1‑100 (MQ=255)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:2596096/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3200515/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:1134503/3‑97 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1472879/3‑100 (MQ=14)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAAT‑‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:2502633/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1230055/3‑100 (MQ=21)
                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  2:409009/1‑100 (MQ=255)
                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                         <  2:265218/100‑1 (MQ=255)
                                                                           ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑‑G‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2732408/11‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCT‑‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCa                        >  2:3306163/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:2083269/2‑96 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGAT‑‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:843489/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGACTACAAAACTTTGCAAATTCAATATATTGCAAT‑‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGAGCATCAGGCa                        >  1:2289709/2‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                      >  2:1851877/1‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:3015468/1‑99 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑C‑‑‑CAC‑‑GTAGGCCTGATAATCGTAGCGAATAAgaaaa                       >  2:1368949/1‑96 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                      >  1:102586/1‑100 (MQ=21)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:681183/9‑98 (MQ=17)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  2:628674/9‑98 (MQ=21)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc                      >  1:2781434/9‑98 (MQ=18)
                                                                             ccGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:958944/1‑100 (MQ=18)
                                                                             caggataaggcgTTCACGCCGCATCCTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGA‑‑T‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  2:195146/9‑99 (MQ=14)
                                                                              ctgataaggcgTTTACGCCTTATCAGGCCTACAAGATCGTGCAAATTCAACATATTGCCACT‑‑C‑‑‑ACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  2:2460863/8‑98 (MQ=25)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGAT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                     >  2:747371/8‑97 (MQ=21)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagct                     >  2:3039669/8‑97 (MQ=18)
                                                                               ggataaggcgTTTACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAA‑‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                    >  1:642239/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCGAATTCAATATATTGCAGAGC‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1900459/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:3447080/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1664766/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAA‑‑‑CTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  2:1937670/7‑97 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCGTACGAAATCGTGCTAATTCAATATATTGTGGGCT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:923982/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGA‑‑C‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagct                     >  1:2986906/7‑96 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                    >  2:1387556/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  2:2554132/7‑98 (MQ=21)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    >  2:2951293/1‑100 (MQ=255)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGAT‑‑C‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattta                   >  1:2575390/6‑99 (MQ=17)
                                                                                 ataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGAT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg                  >  1:374889/5‑100 (MQ=21)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  1:2951246/100‑1 (MQ=255)
                                                                                   aaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcg                >  1:2012771/3‑99 (MQ=17)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCTGGCATTtttgcgt               <  2:2359807/100‑3 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGAGAATCAGGCATTtttgc                 >  1:3231777/1‑97 (MQ=255)
                                                                                      gcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAT‑‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAatttggcgtt              >  2:1765652/1‑94 (MQ=25)
                                                                                      gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATGCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt             <  1:3290504/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:36698/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           >  1:1338954/1‑100 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:1653748/1‑100 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  1:211583/1‑100 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  1:1701729/1‑100 (MQ=255)
                                                                                            aCGCCTTATTCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  1:96064/100‑1 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:3412789/100‑1 (MQ=255)
                                                                                                cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   >  2:2888958/1‑100 (MQ=255)
                                                                                                 ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  >  1:496411/1‑100 (MQ=255)
                                                                                                 ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  2:876350/100‑1 (MQ=255)
                                                                                                 ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑ATC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  1:2197750/1‑98 (MQ=255)
                                                                                                  |                                                                                                            
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083

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