Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,720,828 C→T intergenic (‑4/‑197) slyA ← / → ydhI global transcriptional regulator/DUF1656 family putative inner membrane efflux pump associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,8280CT97.9% 335.5 / ‑2.3 96intergenic (‑4/‑197)slyA/ydhIglobal transcriptional regulator/DUF1656 family putative inner membrane efflux pump associated protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (47/47);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

TTTGTGTTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCATTTGGTT  >  NC_000913/1720731‑1720923
                                                                                                 |                                                                                               
tttGTGTTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTAt                                                                                               <  1:279475/100‑1 (MQ=255)
 ttGTGTTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTAtc                                                                                              <  1:1350396/100‑1 (MQ=255)
  tgtgTTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTAtct                                                                                             <  1:192742/100‑1 (MQ=255)
   gtgtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTAtctc                                                                                            >  2:907086/1‑100 (MQ=255)
   gtgtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTAtctc                                                                                            >  2:282701/1‑100 (MQ=255)
   gtgtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTAtctc                                                                                            >  1:189098/1‑100 (MQ=255)
    tgtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTcc                                                                                           >  1:51499/1‑100 (MQ=255)
    tgtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTcc                                                                                           >  1:1668683/1‑100 (MQ=255)
    tgtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTcc                                                                                           <  1:1809274/100‑1 (MQ=255)
     gtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCt                                                                                          >  1:1034869/1‑100 (MQ=255)
     gtTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCt                                                                                          >  2:1182756/1‑100 (MQ=255)
        aaCTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTtat                                                                                       <  2:1984120/100‑1 (MQ=255)
         aCTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTtata                                                                                      <  2:2006257/100‑1 (MQ=255)
           tCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAAt                                                                                    <  2:1278025/100‑1 (MQ=255)
               gCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGc                                                                                <  2:256418/100‑1 (MQ=255)
                   tttCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGc                                                                            >  2:1770819/1‑100 (MQ=255)
                   tttCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGc                                                                            <  1:1492880/100‑1 (MQ=255)
                    ttCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCGGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCt                                                                           >  1:1779201/1‑100 (MQ=255)
                     tCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTa                                                                          <  2:1668670/100‑1 (MQ=255)
                       aGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAg                                                                        >  2:2000273/1‑100 (MQ=255)
                        ggCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGa                                                                       >  2:72345/1‑100 (MQ=255)
                        ggCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGa                                                                       >  2:1317837/1‑100 (MQ=255)
                        ggCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGa                                                                       >  2:336286/1‑100 (MQ=255)
                         gCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGat                                                                      <  1:1502459/100‑1 (MQ=255)
                             tGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGAtatta                                                                  <  2:232304/100‑1 (MQ=255)
                             tGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGAtatta                                                                  <  2:232497/100‑1 (MQ=255)
                             tGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGAtatta                                                                  >  2:1747626/1‑100 (MQ=255)
                              ggTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGAtattat                                                                 >  1:300843/1‑100 (MQ=255)
                               gTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAACTAGCTTGCTAAGAtattacg                                                                <  2:582987/100‑3 (MQ=255)
                                    tCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCt                                                           >  1:1027923/1‑100 (MQ=255)
                                     cAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCtt                                                          <  1:72345/100‑1 (MQ=255)
                                     cAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCtt                                                          >  2:274921/1‑100 (MQ=255)
                                      agagCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCttt                                                         >  1:1764989/1‑100 (MQ=255)
                                      agagCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCttt                                                         >  1:1387169/1‑100 (MQ=255)
                                       gagCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCtttt                                                        <  1:2033274/100‑1 (MQ=255)
                                        agCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTa                                                       <  1:907090/100‑1 (MQ=255)
                                         gCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAg                                                      <  2:1240917/100‑1 (MQ=255)
                                         gCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAg                                                      <  1:1140081/100‑1 (MQ=255)
                                         gCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAg                                                      <  2:51499/100‑1 (MQ=255)
                                           aCGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAg                                                      <  2:758284/98‑1 (MQ=255)
                                           aCGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGaa                                                    >  2:1305660/1‑100 (MQ=255)
                                            cGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAAt                                                   <  2:201698/100‑1 (MQ=255)
                                            cGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAAt                                                   >  2:377664/1‑100 (MQ=255)
                                              ccATATGCGCACCAACTGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAg                                                 <  2:1199806/100‑1 (MQ=255)
                                              ccATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAg                                                 >  2:1417707/1‑100 (MQ=255)
                                              ccATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAg                                                 >  2:574170/1‑100 (MQ=255)
                                                atatGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAgtg                                               <  1:498324/100‑1 (MQ=255)
                                                  atGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGtgtg                                             >  1:1898083/1‑100 (MQ=255)
                                                   tGCGCACCCACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgc                                            <  1:1452914/100‑1 (MQ=255)
                                                   tGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgc                                            <  2:211645/100‑1 (MQ=255)
                                                   tGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgc                                            <  1:24690/100‑1 (MQ=255)
                                                   tGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgc                                            >  1:662441/1‑100 (MQ=255)
                                                    gcgcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgca                                           >  2:2043992/1‑100 (MQ=255)
                                                    gcgcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgca                                           >  2:445463/1‑100 (MQ=255)
                                                    gcgcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgca                                           <  1:1625285/100‑1 (MQ=255)
                                                     ngcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag                                          >  2:1209790/2‑100 (MQ=255)
                                                     cgcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag                                          >  1:1630235/1‑100 (MQ=255)
                                                     cgcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag                                          <  1:336287/100‑1 (MQ=255)
                                                     cgcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag                                          >  1:1535224/1‑100 (MQ=255)
                                                      gcACCAACCGTGCCAGATCGGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagc                                         <  2:1779186/100‑1 (MQ=255)
                                                      gcACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcag                                          <  2:1647054/99‑1 (MQ=255)
                                                       cACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagca                                        >  1:1353574/1‑100 (MQ=255)
                                                       cACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTgcagca                                        >  1:1901754/1‑100 (MQ=255)
                                                        aCCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCaa                                       >  1:64854/1‑100 (MQ=255)
                                                           aaCCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTg                                    <  1:445464/100‑1 (MQ=255)
                                                           aaCCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTg                                    <  2:189097/100‑1 (MQ=255)
                                                            aCCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGt                                   >  1:444012/1‑100 (MQ=255)
                                                              cGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTAt                                 >  2:1448434/1‑100 (MQ=255)
                                                                tGCCAGATTAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTg                               <  2:1898065/100‑1 (MQ=255)
                                                                tGCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTg                               >  1:455244/1‑100 (MQ=255)
                                                                 gCCAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGc                              <  1:1358537/100‑1 (MQ=255)
                                                                  ccAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCt                             <  2:1480616/100‑1 (MQ=255)
                                                                  ccAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCt                             <  1:270837/100‑1 (MQ=255)
                                                                  ccAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGAGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCt                             >  2:1438353/1‑100 (MQ=255)
                                                                   cAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTa                            >  1:1569379/1‑100 (MQ=255)
                                                                   cAGATCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTa                            <  1:262541/100‑1 (MQ=255)
                                                                      aTCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTaaaa                         >  1:265077/1‑100 (MQ=255)
                                                                      aTCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTaaaa                         <  1:18220/100‑1 (MQ=255)
                                                                      aTCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTaaaa                         >  1:468506/1‑100 (MQ=255)
                                                                      aTCAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTaaaa                         >  2:1504586/1‑100 (MQ=255)
                                                                        cAGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACa                       <  1:1770833/100‑1 (MQ=255)
                                                                         aGAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACaa                      <  2:1366783/100‑1 (MQ=255)
                                                                          gAACCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACaaa                     <  2:1357625/100‑1 (MQ=255)
                                                                            aCCTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATg                   <  2:150703/100‑1 (MQ=255)
                                                                              cTAGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATg                   <  2:1027918/98‑1 (MQ=255)
                                                                                aGTGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTAtt               >  2:62717/1‑100 (MQ=255)
                                                                                  tGGCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATtgc             <  1:1747640/100‑1 (MQ=255)
                                                                                   ggCGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATtgct            >  2:1555498/1‑100 (MQ=255)
                                                                                     cGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGATTATTGTGCAGCAATTTTATTGCTAAAACAAATGTATttcttc          >  2:1971620/1‑98 (MQ=255)
                                                                                     cGATTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATtgctgc          <  1:274922/100‑1 (MQ=255)
                                                                                       aTTCCAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCAt        <  2:1694347/100‑1 (MQ=255)
                                                                                          ccAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCATTTg     >  1:1660178/1‑100 (MQ=255)
                                                                                          ccAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCATTTg     >  1:1804247/1‑100 (MQ=255)
                                                                                           cAATTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCATTTgg    >  2:1964287/1‑100 (MQ=255)
                                                                                            aaTTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCATTTGGt   <  2:1034864/100‑1 (MQ=255)
                                                                                             aTTTTATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCATTTGGtt  >  1:250628/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TTTGTGTTAACTCCAGCGGTTTCAGGCGATGGTCTATCAGAGCACGCCATATGCGCACCAACCGTGCCAGATCAGAACCTAGTGGCGATTCCAATTTCATCTCCTTATAATTAGCTTGCTAAGATATTATGCGGCTTTTAGAATAGTGTGCAGCAATTGTATTGCTAAAACAAATGTATTGCTGCATTTGGTT  >  NC_000913/1720731‑1720923

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: