Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 155,463 A→T *247K (TAA→AAA)  yadV ← putative fimbrial chaperone YadV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913155,4630AT96.2% 163.3 / ‑6.0 52*247K (TAA→AAA) yadVputative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (24/26);  minor base G (0/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACC  >  NC_000913/155218‑155691
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
aCATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                                                                                                                                                                                  <  2:159348/251‑1 (MQ=255)
aCATAAACACTGACGTCGTAAACACCCGGGATAGCGTGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGAGGCATTCTGCCCCATCAGGAGGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATGCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                                                                                                                                                                                   <  2:304451/250‑1 (MQ=255)
 cATAAACACAGACGAAGTAAACACCCGGTATAGCGGGGTTACATTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAAAACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCt                                                                                                                                                                                                                                 >  2:362846/1‑251 (MQ=255)
      acacTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                                                                                                                                                                                            >  2:570910/1‑251 (MQ=255)
                 gtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                                                                                                                                                                                                 >  1:409377/1‑251 (MQ=255)
                  gaaaCACCCGGTATAGCTGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATGCGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCGATGCTTCCTGCTGAATTGTATCCTTGCCTGAGTTCAGCCTGTCATTACCTTCAATTg                                                                                                                                                                                                                <  2:242370/250‑1 (MQ=255)
                   aaaCACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGCTCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTTTATCCTTGCCTGATTTCAGCCTGGCCTTACCTTCAATTgc                                                                                                                                                                                                               >  2:134207/1‑251 (MQ=255)
                           ggTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa                                                                                                                                                                                                       <  2:180158/251‑1 (MQ=255)
                            gTATAGCTGGATTACCATCACTGTACCGTCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGTAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACTGGTCACATTATGATAATTTTTAGTATATTCTATAGTCACGCTATACTTCCTGCGAAATTATATCCTTGCCTGATTTCAGCCTGGCATTTCCTTCAATTGCGCCTCCaaa                                                                                                                                                                                                      <  2:64586/251‑1 (MQ=255)
                             tataGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa                                                                                                                                                                                                       >  2:423488/1‑249 (MQ=255)
                                 gCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAt                                                                                                                                                                                                 >  1:559009/1‑251 (MQ=255)
                                     ggTTACCTTCACTGTACCGGCTTAGATAAATATTAGAGGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACGGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGAGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATATTCACTCTATGCTTCCTGCGTAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCTCCAAAGTCATTaaa                                                                                                                                                                                             <  2:140717/251‑2 (MQ=255)
                                       ttACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATgg                                                                                                                                                                                           >  1:405603/1‑251 (MQ=255)
                                          ccTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGt                                                                                                                                                                                        <  2:573604/251‑1 (MQ=255)
                                                     ccGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGGCATTAATGGCGTAAAAATGCACt                                                                                                                                                                             <  1:27597/251‑1 (MQ=255)
                                                       ggCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACtttt                                                                                                                                                                          <  1:237280/251‑1 (MQ=255)
                                                         cTGAGATCAATAGTAGATGCATTCTGCTCCATCAGTAAGGGATGGTCATATTAAACGAGTTCAACACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGCTTTCAGCCTGGCATTACCTTCAATTGCGCCACCCAAGTCATTAATGGCGTAAAAATGCACTTTTg                                                                                                                                                                         >  2:509811/1‑251 (MQ=255)
                                                            agaTCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAg                                                                                                                                                                      <  2:559009/251‑1 (MQ=255)
                                                            agaTCAATATTAGAGGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATCATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATAATATCCTTTCCTGATTTCAGCCTGGCATTACCTTAAATTGCGCCACCAACGTCATTAATGGCATAAAAATGCACtttttccg                                                                                                                                                                      >  2:339291/1‑247 (MQ=255)
                                                                  atatTAGATTCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGATCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACACTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCg                                                                                                                                                                >  2:456512/1‑251 (MQ=255)
                                                                       aGATGCGTTCGGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTAGCGATACGGGGCAGATGATGATAATTTTTAGTATAGTCTATATGCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATGGCAGCCTTTCATTGCCTGCAATTGCTCCTCCAAAGTCATTAATGGCGTTAAAATGCGCTTTTGcttatggcgctttg                                                                                                                                                           <  2:480594/251‑14 (MQ=255)
                                                                                cTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAgg                                                                                                                                                  >  1:341306/1‑251 (MQ=255)
                                                                                 tGCCCAATCAGGAAGGTATGGTCATATTCACCGCGTTCAGCACTGAAAGAAGTATTGCAATACAGCAGCGCGCAAAACGTGACGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATATTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAACCTGGCATTACCTTCAATTGCGCAACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGACTTCGCTTTGCCATTaaagc                                                                                                                                                 >  2:532326/1‑248 (MQ=255)
                                                                                        tCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGAc                                                                                                                                          >  1:3182/1‑251 (MQ=255)
                                                                                          aGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACtt                                                                                                                                        <  1:515343/251‑1 (MQ=255)
                                                                                           ggAAGGTATGGTCATATTCAACGAGATCAGCACTGAAAGCCGTATTGCAATACAGCAGCGCGCAAAACGTGGCGAAACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCTCTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACttt                                                                                                                                       >  2:540584/1‑251 (MQ=255)
                                                                                            gAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCTCTTTTGCAGAATTCGTTTTGCCATTAAGGCCATTGACTTTc                                                                                                                                      <  2:78688/251‑1 (MQ=255)
                                                                                            gAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTc                                                                                                                                      <  2:461990/251‑1 (MQ=255)
                                                                                                        atatTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGCTACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTAACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGAGCCACCAAAGTCATTAATGGCGTAAATATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCtcatc                                                                                                                           >  2:442329/1‑250 (MQ=255)
                                                                                                        atatTCAACGAGTTCAGCACTGAAAGCAGTAGTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCtcatca                                                                                                                          <  1:305629/251‑1 (MQ=255)
                                                                                                        atatTCAACGAGTTCAGCAATGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCtcatca                                                                                                                          <  2:341306/251‑1 (MQ=255)
                                                                                                            tCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTa                                                                                                                       <  1:570910/250‑1 (MQ=255)
                                                                                                                cGAGTTCAGGAGTGAAAGCAGGATTGCAATACAGGAGCGGGGAAAACGTGGCGATACGGGTCAGTTGTTGATTATTTTTAGTATATTGTATAGTCACGATATGCTTCCTGCGGAATTATATCCTTGGCTGAGTGCAGCCTGGCATTTCCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACGAAATgg                                                                                                                  <  2:156060/251‑1 (MQ=255)
                                                                                                                 gAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGt                                                                                                                 <  1:27204/251‑1 (MQ=255)
                                                                                                                  aGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTg                                                                                                                <  1:289501/251‑1 (MQ=255)
                                                                                                                       aGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATc                                                                                                           >  2:530300/1‑251 (MQ=255)
                                                                                                                               aaGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTcac                                                                                                   <  1:334315/251‑1 (MQ=255)
                                                                                                                                 gCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACAt                                                                                                 <  1:327379/251‑1 (MQ=255)
                                                                                                                                              aCAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGATTTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCg                                                                                    >  2:89010/1‑251 (MQ=255)
                                                                                                                                              aCAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGGGCAATCATTTTCACATCAATCGGATAGCg                                                                                    <  1:323738/251‑1 (MQ=255)
                                                                                                                                                   agcGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTAc                                                                               >  1:148701/1‑251 (MQ=255)
                                                                                                                                                            aaCGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCtt                                                                      >  1:195477/1‑251 (MQ=255)
                                                                                                                                                                    gATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGCCTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAANCGGATAGCGTTTACCGCTAGCTTCTAAATc                                                               >  2:271278/1‑250 (MQ=255)
                                                                                                                                                                        cGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTACTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATcaccac                                                          >  2:181835/1‑251 (MQ=255)
                                                                                                                                                                           ggcAGATGGTGATAAATTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGCTTTCAGCGTTGCATTACCTTCAATTGCGCAACCAAAGTCATTACTTGCGGCAGCATGCACTTTTGCAGAATTCGCTTGGCCATTAAGGCTATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAACTGGAAANCGTTTACCGCAATCTTCTAAATcaccacagc                                                       >  2:284418/3‑248 (MQ=255)
                                                                                                                                                                               gatgatgatAATTTTTAGTAGATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAACAATGCTGTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGATCTTCATGACCTCATCTCTAAATGGTGCAATCAGTTTCACATCAATCGGATAGGGTTTACCGCTAGCTTCTAAATCACCAATACTaaa                                                   <  2:73716/251‑1 (MQ=255)
                                                                                                                                                                                                         ataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCTCCAAAGTCATTAATGGCGTAAAAATGCCCTTTTGCAGAATTCGCTTTGCCATTTAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTg                         <  1:540584/250‑1 (MQ=255)
                                                                                                                                                                                                              cACGCTATGCTTCCTGCGGAATTATAGCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAGTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCAATTCCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCTCCGCTGCTAAAAAGTACGTAGTAAGTGGTTAGATTGGTCACt                   <  2:237058/251‑1 (MQ=255)
                                                                                                                                                                                                               aCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAANGAGACGTAGTAAGGGGTTGGATTGGTCACTc                  >  2:266768/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttCCTGCGGAATTATATCCCTGCCTGATTTCAGCCTGGCATTAACTTCTATTGCGCCACCAATGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCACTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCGATCGGATAGCGTTTACCGTTAGCTTCTAAATCACCACTGCTAAAAGAGACGAANTAAGGGGTTGGTTTGGTACCTCGTACTGa           >  2:281531/1‑249 (MQ=255)
                                                                                                                                                                                                                          ccTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGc        <  2:148701/250‑1 (MQ=255)
                                                                                                                                                                                                                               ggAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTAcc  <  2:405603/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
ACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACC  >  NC_000913/155218‑155691

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: