Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 29.1% intergenic (‑181/+127) 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,9710CA29.1% 506.1 / 212.7 237intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (69/99);  new base A (62/7);  total (131/106)
Fisher's exact test for biased strand distribution p-value = 8.48e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 253 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070
                                                                                                   |                                                                                                           
aGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                                                             <  2:3824800/100‑1 (MQ=255)
    aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                         <  1:2273616/100‑1 (MQ=255)
    aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                         <  1:2311326/100‑1 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                      >  1:2698222/1‑100 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGAATGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                       <  2:2632717/99‑1 (MQ=255)
          gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                   <  1:821901/100‑1 (MQ=255)
          gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                   <  1:1494949/100‑1 (MQ=255)
          gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                   <  2:3586213/100‑1 (MQ=255)
             ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                >  2:1812275/1‑100 (MQ=255)
             ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                >  1:60143/1‑100 (MQ=255)
               gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                              <  1:1774437/100‑1 (MQ=255)
                ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                             >  1:2797547/1‑100 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                           >  2:1896680/1‑100 (MQ=255)
                   caggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                          <  1:2549909/98‑1 (MQ=255)
                   caggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                          <  1:3440299/98‑1 (MQ=255)
                    cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                           >  2:276027/1‑98 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                        <  1:1981352/100‑1 (MQ=37)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTGTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                       <  2:295458/100‑1 (MQ=255)
                        tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑a                                                                                    >  1:2932017/1‑100 (MQ=255)
                             cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑aaaa                                                                              >  2:2187834/1‑100 (MQ=255)
                              tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAAt                                                                             <  1:3629351/100‑1 (MQ=255)
                              tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAAt                                                                             >  1:3256633/1‑100 (MQ=255)
                               gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATc                                                                            >  1:1916829/1‑100 (MQ=255)
                                gaCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTACGCTATGTTTATGCCGAATGCGTCTTTAACGCCTTATCCGTCCT‑Acgca                                                                                >  2:1485222/2‑93 (MQ=25)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                           <  1:3510278/100‑1 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                           <  2:2286037/100‑1 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                           <  1:3017425/100‑1 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                           <  2:3814777/100‑1 (MQ=255)
                                 aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGt                                                                          <  2:834424/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTg                                                                         <  1:3787153/100‑1 (MQ=255)
                                        ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAAtt                                                                   <  1:3346184/100‑1 (MQ=255)
                                         gCACAGCGCTACCGGGTATTCAGATGGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTc                                                                  <  2:1924443/100‑1 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCa                                                                 >  1:1399735/1‑100 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAAt                                                               <  1:1250809/100‑1 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATa                                                              <  2:3482615/100‑1 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATa                                                              <  2:128994/100‑1 (MQ=255)
                                                  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTg                                                         >  2:3650752/1‑100 (MQ=255)
                                                     cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑g                                                    <  2:1039507/100‑1 (MQ=255)
                                                        gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑At                                                >  2:146694/1‑100 (MQ=255)
                                                          ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAACTTCAATACATTG‑CA‑GA‑At‑tt                                             <  2:2411721/100‑1 (MQ=255)
                                                          ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑At‑tt                                             <  2:433215/100‑1 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑At‑ttt                                            >  1:1104764/1‑100 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑At‑ttt                                            >  2:1644592/1‑100 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑At‑ttt                                            >  2:1259000/1‑100 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAATGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑g                                        >  2:1517042/1‑100 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑g                                        <  2:1370813/100‑1 (MQ=255)
                                                                 ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAgg                                    <  1:2593722/100‑1 (MQ=255)
                                                                    ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCt                                 >  1:1123435/1‑100 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTg                                <  1:183132/100‑1 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTg                                <  1:3828759/100‑1 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                               >  2:2556248/1‑100 (MQ=255)
                                                                       aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGAt                              >  1:1326639/1‑100 (MQ=255)
                                                                          aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAg                           <  1:1216033/100‑1 (MQ=255)
                                                                          aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAg                           >  1:302455/1‑100 (MQ=255)
                                                                           aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                          <  1:723624/100‑1 (MQ=255)
                                                                             tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                        <  1:239360/100‑1 (MQ=255)
                                                                                gcTTTGTTTATGCCGGATGCGTCGTTAATGCCTTTTTTTGTTT‑AC‑AAAATCGCGCCAACTCAAAACATTG‑CA‑GA‑AG‑TTTC‑‑ATAGGCCTGATAAGCGCAg                      <  2:1717930/99‑1 (MQ=12)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGAtgcgctccgct                    >  2:2467479/1‑90 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                    >  2:1503135/1‑100 (MQ=255)
                                                                                    tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATTGTGTAAATTCAATAAATTG‑AT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  2:2006015/8‑100 (MQ=21)
                                                                                    tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAC‑AAAAACCTGCAAATTCAATATATTG‑CA‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCa                  >  1:3159563/8‑100 (MQ=12)
                                                                                     ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa               >  2:3139033/7‑100 (MQ=18)
                                                                                     gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  2:3528445/100‑1 (MQ=255)
                                                                                     gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  1:875368/1‑100 (MQ=255)
                                                                                     gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:2122604/1‑100 (MQ=255)
                                                                                      tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:3694331/100‑1 (MQ=255)
                                                                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:2980829/100‑1 (MQ=255)
                                                                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  2:1140647/100‑1 (MQ=255)
                                                                                       gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCAAATTCAACATATTG‑CA‑AT‑TC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:2104335/5‑100 (MQ=18)
                                                                                       gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:2324274/5‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1637322/4‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑AC‑AAAAGCACGCAAATTCAATATATTG‑CA‑GA‑GA‑TTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:813508/4‑100 (MQ=14)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:890024/4‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATTGTGTAAATTCAATAAATTG‑AT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             <  1:2006025/97‑1 (MQ=21)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACCTGCAAATTCAATATATTG‑CA‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:3612369/4‑100 (MQ=18)
                                                                                        tATGCCGGCTGCGGCGTTAACGCCTTCTCCGGCCT‑AC‑TAAATCGTGCAAATTCAATACATTG‑CA‑TA‑AT‑TTTC‑‑GTATTGCTGATAAGCGGAGCGCATCAgg             <  2:2712831/100‑1 (MQ=17)
                                                                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCGTGCAAATTCAATACATTG‑CATGG‑GC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:2734088/3‑100 (MQ=17)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc            >  1:2884127/3‑100 (MQ=21)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:2163579/3‑100 (MQ=21)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:70706/3‑100 (MQ=17)
                                                                                          tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:3904781/1‑100 (MQ=255)
                                                                                          ggCCGGATAAGGGGCTTACGCCGCATCCGGCAT‑AC‑AAGAACGTGCAAATTCAACAAATTGTGATGC‑AT‑CCT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:746507/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:45136/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:135189/2‑100 (MQ=17)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1747833/2‑100 (MQ=17)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:2326511/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1535773/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCAAATTCAACATATTG‑CA‑AT‑TC‑TCTT‑‑GTAGGTCTGATAAGCGTAGCGCATCAGGctg          >  1:563498/1‑98 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa          >  1:1590760/1‑99 (MQ=21)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:592406/1‑99 (MQ=21)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:2994279/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACCTTGCTAATTCAATATATTG‑CA‑GGGAC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag          >  1:2198249/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  2:202301/1‑100 (MQ=255)
                                                                                            ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGccta         <  2:1588541/100‑4 (MQ=255)
                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAACGCATCAGGc            >  2:1849155/1‑97 (MQ=17)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:3235920/100‑1 (MQ=255)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:302455/100‑1 (MQ=255)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:3804108/100‑1 (MQ=255)
                                                                                                  gcggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  1:985615/100‑1 (MQ=255)
                                                                                                   cggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttttttct  >  2:1180819/1‑96 (MQ=255)
                                                                                                   |                                                                                                           
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070

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: ‑