Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 32.4% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG32.4% 432.2 / 272.4 241intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (69/94);  new base G (73/5);  total (142/99)
Fisher's exact test for biased strand distribution p-value = 9.70e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

Only 100 of 248 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                               
ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                                 <  1:2139904/100‑1 (MQ=255)
 tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                                >  1:1022670/1‑100 (MQ=255)
  gctgcCGGATGACGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                               <  1:1093019/100‑1 (MQ=255)
   ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                              >  2:3712696/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                            >  2:2066277/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                            >  2:1998442/1‑100 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                          >  1:803558/1‑100 (MQ=255)
        ggATGGCGCTACACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTCTACGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                         <  1:2461596/100‑1 (MQ=255)
          gtGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                       >  2:2910299/2‑100 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                       <  2:382467/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                       <  1:3000303/100‑3 (MQ=25)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                       <  2:1302352/97‑3 (MQ=25)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑aa                                                                                                 <  2:1804473/100‑1 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATc                                                                                             >  2:1201301/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                                            <  2:505950/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                                            >  1:774888/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                                            <  2:2025669/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTg                                                                                          >  1:2574668/1‑100 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGc                                                                                         <  2:2958851/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGc                                                                                         <  1:2846191/100‑1 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGc                                                                                         <  1:1991987/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCT‑AC‑AAAATCGTGCa                                                                                        >  1:3459258/1‑100 (MQ=255)
                           ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAAtt                                                                                    <  1:442051/100‑1 (MQ=255)
                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCaa                                                                                 <  2:3577927/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAAt                                                                                <  2:401821/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAAt                                                                                <  1:2910280/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATAc                                                                              <  1:1746054/100‑1 (MQ=255)
                                   gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAt                                                                            >  1:1187763/1‑100 (MQ=255)
                                   gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAt                                                                            <  2:2264695/100‑1 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                           >  1:112978/1‑100 (MQ=255)
                                        cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                                    <  1:998186/100‑1 (MQ=255)
                                            tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t                                                              >  2:3153657/1‑100 (MQ=255)
                                            tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t                                                              >  1:1888014/1‑100 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t‑‑‑t                                                          <  2:138914/100‑1 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t‑‑‑t                                                          <  1:3570604/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑g                                                     <  2:646015/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑Gt                                                    >  2:3621596/1‑100 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTa                                                   <  1:3621568/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTa                                                   <  1:1220134/100‑1 (MQ=255)
                                                   aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAg                                                  <  1:2373010/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTg                                             >  2:282792/1‑100 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCAAATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTg                                             <  1:3108082/99‑1 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATa                                          <  2:1156422/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGc                                       <  1:1661067/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGc                                       >  2:2984976/1‑100 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     >  1:3081429/1‑100 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     >  2:2369793/1‑100 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                   <  2:2601910/100‑1 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  1:1773930/100‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  1:372536/1‑100 (MQ=255)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  2:1291226/1‑100 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑ACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  2:251765/7‑100 (MQ=18)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             <  1:1275269/100‑1 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             >  2:2665227/1‑100 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             <  2:2556385/100‑1 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             >  2:1916605/1‑100 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            <  1:1201295/100‑1 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            >  2:500551/1‑100 (MQ=255)
                                                                         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:2562065/6‑100 (MQ=18)
                                                                         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  1:961340/6‑100 (MQ=17)
                                                                          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:1171033/100‑1 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:1528971/5‑100 (MQ=17)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑AC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:494504/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           >  1:2095719/5‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:3339126/4‑100 (MQ=21)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑TA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3550985/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:3307256/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3500583/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑A‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                        >  1:1938436/3‑99 (MQ=21)
                                                                            aTGCCGGATGCGACGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  2:2447144/1‑100 (MQ=255)
                                                                             ggtcggataaggcgTTCACGCCGCATCCGGCCC‑AC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2762599/11‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  2:3188461/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1510802/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GA‑GA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2698/2‑100 (MQ=17)
                                                                             gaccggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:3456266/11‑100 (MQ=18)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       <  2:3267846/100‑1 (MQ=255)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                        >  2:2232253/1‑97 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:1096657/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:1675406/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GA‑GA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:2602598/1‑99 (MQ=17)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:746243/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  1:2325546/9‑98 (MQ=25)
                                                                               ccGGATAAGGCGTTTACGCCGCATCTGGCCT‑AC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:1059579/1‑98 (MQ=17)
                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:87263/1‑98 (MQ=18)
                                                                                cGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                     >  2:237122/1‑100 (MQ=255)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1769186/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtat                    >  1:2271363/7‑96 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑AC‑AAAAGCATGAAAATTCAATATATTGC‑‑A‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1724518/7‑100 (MQ=18)
                                                                                 ggATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    <  1:1389415/100‑1 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:774892/100‑1 (MQ=255)
                                                                                     gcggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                <  1:2951913/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:3567587/1‑100 (MQ=255)
                                                                                        gcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt             <  2:1649900/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  2:1370108/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAGTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:2169642/100‑1 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  1:2798452/100‑1 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccagatcggaa       >  2:1726952/1‑91 (MQ=255)
                                                                                               cGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      <  1:1441139/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:3389178/100‑1 (MQ=255)
                                                                                                 ccTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  2:1884261/1‑100 (MQ=255)
                                                                                                   |                                                                                                               
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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