Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,900,610 C→T W537* (TGG→TGA bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,900,6100CT87.2% 108.0 / 7.0 47W537* (TGG→TGAbglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (20/21);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCC  >  NC_000913/3900367‑3900841
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
ccggaaccggaaCACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGAttt                                                                                                                                                                                                                                  <  1:434598/251‑1 (MQ=255)
        ggaacACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGAGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATAGGATTTCAGCATgg                                                                                                                                                                                                                          >  1:492081/1‑251 (MQ=255)
             acccaccACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGCAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCACTTTTTATATTCATCGTCATATTCCTTTCATTCAGTTCCATTACCATCAGATTTCAGCCTGGGCAcc                                                                                                                                                                                                                     >  2:552534/1‑251 (MQ=255)
                caccacTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGAc                                                                                                                                                                                                                  >  1:467852/1‑251 (MQ=255)
                        cAGAGACATTTTTTGCCCCAGGGGAAAAGTAGCGATAGGTCACGCTTGTACCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTTCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAGTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACt                                                                                                                                                                                                          <  2:38504/251‑1 (MQ=255)
                                   tttGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCttt                                                                                                                                                                                                >  1:440506/1‑250 (MQ=255)
                                     tGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTAt                                                                                                                                                                                              <  1:552534/250‑1 (MQ=255)
                                      gCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTg                                                                                                                                                                                            <  1:506408/251‑1 (MQ=255)
                                      gCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTg                                                                                                                                                                                            >  1:225679/1‑251 (MQ=255)
                                                cAAAGTAGCGTAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAgg                                                                                                                                                                                  >  1:427105/1‑251 (MQ=255)
                                                cAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAgg                                                                                                                                                                                  >  1:416198/1‑251 (MQ=255)
                                                 aaaGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGt                                                                                                                                                                                 >  1:250756/1‑251 (MQ=255)
                                                                     tGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGCTTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGtttt                                                                                                                                                             >  1:58792/1‑251 (MQ=255)
                                                                      gTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTc                                                                                                                                                            <  1:423566/251‑1 (MQ=255)
                                                                         aTCGGCCTTGCCCTGAGCCACATACTGATTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATCTCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCCGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAATGGTGAAGCCATCCAGTTCCTTTTcccg                                                                                                                                                         >  2:365847/1‑249 (MQ=255)
                                                                                gttgcctggatcccaTACTGTGTTCCCGGGATAGTGGGTGCCGGGTTGGCTGGCTTATCCGCTGCTATGGTCCATTTAGCAGAATTACCGCTGGTAATAGTATGCGTTAAAGCAGCAATTTTTATATTCATTGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGGCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCttt                                                                                                                                                   <  2:300264/237‑1 (MQ=255)
                                                                                   aCCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGAAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGACTTTATTGTCTTCAAAGGTGAAGCCATCCCGTTCGTTTTCCAGGGCTTTGatt                                                                                                                                               >  2:13632/1‑250 (MQ=255)
                                                                                              aTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCga                                                                                                                                     >  2:130232/1‑249 (MQ=255)
                                                                                                                 ggtggggtccGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAAGTTCCGGatg                                                                                                                  <  2:416198/243‑3 (MQ=255)
                                                                                                                 ggTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGAc                                                                                                                   <  2:250756/249‑1 (MQ=255)
                                                                                                                 agtgggttcctggcctgctggcATCTCTGCTGTTATGGCCCACTGTGCAGAAATACCGCTGGCAATAGCTATCGTTAACGCATCAATTATTATACTCATCGTAATATTCCTTTCATGAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCTAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGTGCTTTGAGATACGTGGCGTAGAAACGAATGTCCGGAcg                                                                                                                  <  2:297346/235‑1 (MQ=255)
                                                                                                                    gggggCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACGGAGCAGAAATACCGCTTGCAATTGCCAGCGTTAAAGCATCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCGTCAAAGGTGAAGCCATCCAGTTCGTTGTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGAcgcg                                                                                                                <  2:372087/246‑1 (MQ=255)
                                                                                                                      gTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAgg                                                                                                             <  1:389352/251‑1 (MQ=255)
                                                                                                                        gCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTc                                                                                                           >  1:225526/1‑251 (MQ=255)
                                                                                                                                 gctggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTg                                                                                                  <  1:13632/251‑1 (MQ=255)
                                                                                                                                                                   aaaTACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGttt                                                                <  1:116231/251‑1 (MQ=255)
                                                                                                                                                                     aTACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCTTCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGtttt                                                               >  1:275113/1‑250 (MQ=255)
                                                                                                                                                                     aTACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGtttt                                                               >  1:14079/1‑250 (MQ=255)
                                                                                                                                                                             ggCAATAGCCAACGTTAAAGCAGCAAATTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCTCCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCa                                                        <  2:14079/249‑1 (MQ=255)
                                                                                                                                                                             ggCAATAGCAAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTGTCATTAAGTGCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCTTCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCa                                                        <  2:275113/249‑1 (MQ=255)
                                                                                                                                                                                aaTAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCttt                                                                                                                                                   <  2:140453/155‑1 (MQ=255)
                                                                                                                                                                                aaTAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCttt                                                                                                                                                   >  1:140453/1‑155 (MQ=255)
                                                                                                                                                                                 aTAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCa                                                  <  2:225526/251‑1 (MQ=255)
                                                                                                                                                                                   aGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCAt                                                 <  2:467852/250‑1 (MQ=255)
                                                                                                                                                                                    gCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTa                                               <  2:197095/251‑1 (MQ=255)
                                                                                                                                                                                         cGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGtt                                                                                                                                                               >  1:32226/1‑134 (MQ=255)
                                                                                                                                                                                         cGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTGCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGtt                                                                                                                                                               <  2:32226/134‑1 (MQ=255)
                                                                                                                                                                                            tAAAGCAGCAATGTTTATATATATCGTAATATTCCTTTCAATAAGTTCCATGACCATCAGATTTCAGCCTGGGCACCGACATCAAACGGTTCGTTTTTATTGGCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGGGGCGTAGAAACGAATTTCCGGACGCGAGTTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACGCATTAAATttat                                        <  2:351864/250‑1 (MQ=255)
                                                                                                                                                                                                  agcaATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTAtt                                 <  1:174361/251‑1 (MQ=255)
                                                                                                                                                                                                      aTTTTTATATTCATCGTAATATCCCTTTCATTAAGTTCAATTACCTTCAGATTTCAGCATGGACACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAACGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAAAATACTGGTATTGACTTTAACGATATGGAAGAGCGGGGTTTTATAACCAGAATCATTAATTTTATTACTATTCGc                              >  2:77931/1‑250 (MQ=255)
                                                                                                                                                                                                             tatTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATttatta                                      >  2:67654/1‑235 (MQ=255)
                                                                                                                                                                                                             tatTCATCGTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATttatta                                      <  1:67654/235‑1 (MQ=255)
                                                                                                                                                                                                                     gTAATATTCCTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTgg              >  1:522783/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGCCTCATTAAATTTATTACTATTCGCATCTTTGtt                     >  2:437100/1‑236 (MQ=255)
                                                                                                                                                                                                                             ccTTTCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGtt                     <  1:437100/236‑1 (MQ=255)
                                                                                                                                                                                                                             catttaaTTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGCATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACAAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGGGAAAt       >  2:154983/7‑250 (MQ=255)
                                                                                                                                                                                                                                 tCATTAAGTTCCATTACCATCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGCTATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGAATCTTTGTTTTGCTGGGTGAAATAGccc  >  1:337780/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
CCGGAACCGGAACACCCACCACTACAGAGACATTTTTTGCCCCAGGGGCAAAGTAGCGAAAGGTCACGCTGTTATCGGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCC  >  NC_000913/3900367‑3900841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: