Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,445,315 C→T P402S (CCC→TCC)  tamB → translocation and assembly module subunit TamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,445,3150CT95.8% 133.0 / ‑2.4 51P402S (CCC→TCC) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (21/25);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

GGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTG  >  NC_000913/4445066‑4445553
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
ggcgCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATcccgcc                                                                                                                                                                                                                                                  <  1:65535/251‑1 (MQ=255)
        gcgcGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACt                                                                                                                                                                                                                                           <  1:73106/250‑1 (MQ=255)
             aaCAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTAcc                                                                                                                                                                                                                                     <  1:323414/251‑1 (MQ=255)
             aaCAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTAcc                                                                                                                                                                                                                                     <  1:99609/251‑1 (MQ=255)
                atgtGGAAATTGGCGTTAATCGTTCACGTCCGGTGGATAGGGATTTACGCGCCAAGGCGCGACTGGCAGAAGCCGGATTGACGCTCAACGTGGAAGTGATGAGCAAACAACATTACTGGCTGGACACTGGTGAGAAGCAGTATCAGGCGGATGATCTAAAGCTGCAGCTTACCGGCAGAATGACCGATTACACGCTCTCTATGCGTACGGAAGTGAAGGGACAGGAGATCTCGGCCGCGACTATTGcccgc                                                                                                                                                                                                                                  <  2:80703/248‑3 (MQ=255)
                 gCTGGAGATTGTCGTTAAGCTTGCCGGTCCGGTGAATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAATGTGGAAGTGAACAGCAAAAAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATTATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCACTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCg                                                                                                                                                                                                                                 <  2:270965/251‑1 (MQ=255)
                 gCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCg                                                                                                                                                                                                                                 >  1:178343/1‑251 (MQ=255)
                      agaTTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGcc                                                                                                                                                                                                                            <  1:87554/251‑1 (MQ=255)
                             cGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTa                                                                                                                                                                                                                     <  1:49733/251‑1 (MQ=255)
                                                gtGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAAAAGCTTTACTGGCCGTTAACTGGAGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTATCTATGCGTACGGCAGTGGAAGGAAAGGAGATCCCGTCAGCGACTATTAACCTCGACGCCAAAGGTAATGAACAGCAGGTCAAtat                                                                                                                                                                                                 >  1:151501/1‑249 (MQ=255)
                                                  ggATATGGATTTACTCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCg                                                                                                                                                                                               <  2:217552/251‑1 (MQ=255)
                                                        ggATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAAc                                                                                                                                                                                         <  1:60161/251‑1 (MQ=255)
                                                         gATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGAAGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACt                                                                                                                                                                                        >  1:368790/1‑251 (MQ=255)
                                                            ttACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCAc                                                                                                                                                                                     >  1:80089/1‑251 (MQ=255)
                                                                   cccAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTcgcg                                                                                                                                                                              <  1:63465/251‑1 (MQ=255)
                                                                    ccAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCCCCGTCgcgg                                                                                                                                                                             >  1:153533/1‑251 (MQ=255)
                                                                     cAGACGCGCCTGGCGGAAGCCGGTTTGCCGCTCAACGTTGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTTGTGAGAAGCAGTATCAGGCTGATTATCTGAAACTGAAACTTACCGTCAAAATGACCGATGACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTTCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCgcggc                                                                                                                                                                            <  1:372726/251‑1 (MQ=255)
                                                                        aCGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCt                                                                                                                                                                         <  1:14514/251‑1 (MQ=255)
                                                                           cgACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGa                                                                                                                                                                      <  1:111889/251‑1 (MQ=255)
                                                                                          ggATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACt                                                                                                                                                       <  1:152624/251‑1 (MQ=255)
                                                                                            aTTGCCGCTCAACGTGGAAGTGAACAGCACACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGCAACTGAAACTTACCGGCAAACTGACCGATTACACGCTCTCTAAGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTAGACGCCTAAGGCAATGAAAAGCAGGGCAATCTCGACAAACTGACCGTCGCGGCGCTGGCAGGGAAAACTGAACTc                                                                                                                                                      >  1:155320/1‑250 (MQ=255)
                                                                                            aTTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTc                                                                                                                                                      >  1:92223/1‑250 (MQ=255)
                                                                                              tGCCGCTCAACGTGGAAGTGAACAGCAAACAGCGTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGAGCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCaa                                                                                                                                                    <  1:165933/250‑1 (MQ=255)
                                                                                              tGCCGCACAACGAGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGGATCAGAAGGATGATCTGAAACTGAAAATTACTGGCCAAATGACAGATTACACGATCACTATGCGGACGGCAGTGAAGTGACAGGAGATCCCGTCCGCGACCAATAACCTCGACGCCAAAGGTAATGAACACACAGTCAATCTCGACAAACTCACCGTAGCGGCGATGGAGGGGcaaagtagactcaa                                                                                                                                                    >  2:9085/1‑240 (MQ=255)
                                                                                               gccgcataaagtagAAGTGAACAGCAAACAGATTGTATGGTCGGTCCCTGGTGAGAAGCAGGATCAGGCGTATGATCCGAAACGGAATCTTACCGGCAAAATGACCGATTACACTCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGAACCCGACCGCGACTATTACCCTTGACGCCAAAGGAAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACGGAACTCAAgg                                                                                                                                                  <  2:303045/238‑1 (MQ=255)
                                                                                                          gggAAGTGTACAGCAACCAGCTTTACGGGGCATTCACTGGTGAGAAGCAGTACCAGGCGGATGATCTGAGACTGAAACGTCCCGGCAAAAAGACCGCTCACCCGCGCTCTATGCGTACGGCAGTGAAGGGACAGGAGAGCCCTTCCGTGACTATTACCCTCGACGGCACAGGTAATGAACAGCAGGTCAATCTCTACAAACTCACCGTCGCGGCGCTGGAAGTGAAATCTGAACTCAAGTCGTTGCTCGAc                                                                                                                                       <  2:19131/250‑1 (MQ=255)
                                                                                                           ggAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAATGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACCGCAGGTCAATATCGAAAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACt                                                                                                                                      >  2:343714/1‑251 (MQ=255)
                                                                                                                       cAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGcc                                                                                                                           >  1:82167/1‑250 (MQ=255)
                                                                                                                         aaCAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGACATCACGGCGTTGCTCGACTGGCAACAGGCCAtt                                                                                                                        >  1:259460/1‑251 (MQ=255)
                                                                                                                                                tGGTGAGAAGCAGTAGCAGGCGGATGATCTGAAACGGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATTCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCa                                                                                                                                                                                                                           <  1:165654/131‑1 (MQ=255)
                                                                                                                                                tGGTGAGAACCAGTATCAGGCGGATGATCTGAAACTGAAACTTCCCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGGCCGCGACTATTACCCTCGCCGCCa                                                                                                                                                                                                                           >  2:165654/1‑131 (MQ=255)
                                                                                                                                                       aaGCGGTATCATGCGGATGAGCGTAAACTGAAACGTACCGGCAAAATGACCGATTACACGCTCTCTAGTCGTACGGCAGTGAAGGGAAAGGAGATCCCGTCCGCGACTATTCCCCGCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTTGCCGTCGCGGCGCTGGAAGTGAAAACTGAACTCAAGGCGGTTCGCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGATTAACgg                                                                                          <  2:178343/251‑1 (MQ=255)
                                                                                                                                                            gTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGAGCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTa                                                                                     <  1:121449/251‑1 (MQ=255)
                                                                                                                                                             tATCAGGCGGATGATCTGAAACTGCAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACCGGCGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTa                                                                                     >  1:259074/1‑250 (MQ=255)
                                                                                                                                                              aTCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTaa                                                                                    <  1:102414/250‑1 (MQ=255)
                                                                                                                                                                    cGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGcc                                                                             <  1:343714/251‑1 (MQ=255)
                                                                                                                                                                               aaactgaaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTccc                                                                 >  1:342433/1‑251 (MQ=255)
                                                                                                                                                                                aactgaaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCg                                                                <  1:358204/251‑1 (MQ=255)
                                                                                                                                                                                aactgaaactTACCGGCAAAATGACCAATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGc                                                                                                        <  1:168396/213‑1 (MQ=255)
                                                                                                                                                                                     aaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTg                                                           >  1:209929/1‑251 (MQ=255)
                                                                                                                                                                                     aaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTg                                                           >  1:70461/1‑251 (MQ=255)
                                                                                                                                                                                        ctTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACgg                                                                                          >  1:285261/1‑218 (MQ=255)
                                                                                                                                                                                                  aaaTGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTc                                              >  1:93314/1‑251 (MQ=255)
                                                                                                                                                                                                    aTGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCTCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGAAAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAacac                                                                                 <  1:10250/215‑1 (MQ=255)
                                                                                                                                                                                                      gACCGATTACCCGCTCTCTAGGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCACGACGCCAAAGGTAATGAACAGCAGGTCAAGCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACGGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATg                                          <  2:70461/251‑1 (MQ=255)
                                                                                                                                                                                                            ttACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTCCCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCaa                                            <  1:228021/243‑1 (MQ=255)
                                                                                                                                                                                                            ttACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCAGCGACTATTACCCTCGACGCCAAAGGTAATGAACCGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGATGGAAGGGAAAAATGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGAGAGCTAACGCTTAACGGCATTAAAACCGCCAAAGAGTGCCCGGAATGGCCGTcagaaatcaaagatcggg                                    >  2:228021/1‑234 (MQ=255)
                                                                                                                                                                                                               cacGCTCTCTATGCGTACGGCAGTGAAGGTACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTACTGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGCAACCTGAACTCAAGACGTTACTCGCCTGGCCACAGGCCCTTAGTTGGCGCGGTGCGCTAACGCTTACCGGCATTAACACCCCACAAGAGTTCCCGGACTGCCCGTCggacctccctgtcttagtg                                  >  2:71452/1‑232 (MQ=255)
                                                                                                                                                                                                                  gctatctATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTAGACAAACTCACCGTCGCGGCGCTGGAAGGGAACACTGAGCTCAAGGCGATGCTCAACGCGCATCAGGCCATTAGTGGGGACGCGCAGCTAAGCGTTAACGTCATGCACACAGCAGAAGAGTTCCCGGACTCGACGGCGAAAGTCAATGGCTTGATTaat                               >  1:250162/5‑249 (MQ=255)
                                                                                                                                                                                                                       ctATGCTTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCCCGACGCCAAAGTTAATGAACAGCAGGTCAATCTCGACAAACTCTCCGTCGCGTCGCTGGAAGTGAAAACTGAACGCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCAGTTTGCTAACGCTTAACGGCATTAACACCGCCAAATAGTTCCCGGACTGGCCGTCGAAACTCAATGGCGTGATTAAAAcc                            <  1:190645/248‑1 (MQ=255)
                                                                                                                                                                                                                            cGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAg                    <  1:315409/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cAGGAGATCCCGTCCGAGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTg  >  1:308097/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
GGCGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTG  >  NC_000913/4445066‑4445553

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