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,3150CT97.2% 224.8 / NA 74P402S (CCC→TCC) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (36/34);  total (37/37)
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

CGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGC  >  NC_000913/4445070‑4445555
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
cgcTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCcgcg                                                                                                                                                                                                                                             <  1:376686/251‑1 (MQ=255)
      gcgAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTAtt                                                                                                                                                                                                                                       <  1:383845/251‑1 (MQ=255)
         aaCAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAAATTACCGGCAAAATGACCGATTACACGCTCTCTTTGCGTACGGCAGTGAAGGGACAGGGGATCCCGTCCGCGGCTATTAcc                                                                                                                                                                                                                                    >  2:376700/1‑251 (MQ=255)
           cAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCt                                                                                                                                                                                                                                  <  1:7238/251‑1 (MQ=255)
                ggAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACg                                                                                                                                                                                                                             <  2:94906/251‑1 (MQ=255)
                 gagaTTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGGTTGCCGCTCAACGTGGAAGTGAACAGAAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGGTGATCTGAAACTGAAACTTACCGGCAAAATAACCGATTACAAGCTCTCTCTGCGTAAGGCAGTGAAAGGGCAGGGGAACCCGTCCGCGAACATTAcccccgcccc                                                                                                                                                                                                                            >  2:264629/1‑244 (MQ=255)
                       ggCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAg                                                                                                                                                                                                                       <  1:369385/250‑1 (MQ=255)
                       ggCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAg                                                                                                                                                                                                                       <  1:244809/250‑1 (MQ=255)
                                  ttCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAAcag                                                                                                                                                                                                            <  2:410480/250‑1 (MQ=255)
                                  ttCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAAca                                                                                                                                                                                                             <  1:244607/249‑1 (MQ=255)
                                         ccggtGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCa                                                                                                                                                                                                     <  1:355698/250‑1 (MQ=255)
                                          cggtGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAAt                                                                                                                                                                                                   <  1:311159/251‑1 (MQ=255)
                                            gtGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGAACGATTAAACGCCCTCTATGCGGACGGCAATGAAGGGACAGGAGATCCCGGACCGGACTATTAACCTAGCAGCCAAAGGGAATGAAcaaccaggcaatac                                                                                                                                                                                                 >  2:133874/1‑239 (MQ=255)
                                            gtGGATATGGAGTTCCGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAAtct                                                                                                                                                                                                 <  2:200183/251‑1 (MQ=255)
                                                          aCGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCg                                                                                                                                                                                   >  1:299914/1‑251 (MQ=255)
                                                              gcCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGCCGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTcgc                                                                                                                                                                               >  1:257665/1‑251 (MQ=255)
                                                                              cGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGCAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGaa                                                                                                                                                                >  2:24259/1‑250 (MQ=255)
                                                                                 aaGCCGGATTGCAGCTAAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGACGGATGATATGAAACTGAAAATTACCGGCAAAATGAACGATTAAACGCTCTCTATGCGGACGGCAGTGAAGGGAAAGGAGGTCCCGTCCGCGACTCTTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGGCACACTTACCGTCGTGGCGTTGGATGGGAAAAc                                                                                                                                                             >  2:216103/1‑250 (MQ=255)
                                                                                  aGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGCAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTg                                                                                                                                                           >  1:467360/1‑251 (MQ=255)
                                                                                     cGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAAc                                                                                                                                                        <  1:303736/251‑1 (MQ=255)
                                                                                          tGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCaa                                                                                                                                                    <  2:282546/250‑1 (MQ=255)
                                                                                           gCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGGATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAg                                                                                                                                                   <  2:23702/250‑1 (MQ=255)
                                                                                             cGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCg                                                                                                                                                >  1:190988/1‑251 (MQ=255)
                                                                                             cGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAACAGGTCAAACACGACAAACCAACCGACGCGGCGCTGGGAGGGAAAACTGAACTCAAgg                                                                                                                                                  >  2:350095/1‑249 (MQ=255)
                                                                                                 cAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGGATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATTAACAGCAGGTCAATCCCGACAAACTCACCGTCGCGGCGCTGGACGGGTAAACTGAACTCAAGGCGTTGc                                                                                                                                            >  2:115657/1‑251 (MQ=255)
                                                                                                  aaCGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAACCTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGAGTACGGCAGTGACGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGCCGTTGCt                                                                                                                                           >  1:62020/1‑251 (MQ=255)
                                                                                                        gAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTg                                                                                                                                     <  1:262401/251‑1 (MQ=255)
                                                                                                                 aGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACCCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGAGTTGCTCGAATGGCAACAGGc                                                                                                                            >  2:3208/1‑251 (MQ=255)
                                                                                                                    aaaCAGCGTTACTGGCCGTTCACTGGTGAGAAGCAGTAGCAGGCGGATGATCTGAACCTGAAACTTACCGGCAAAATGACCGATTGGACGCTCTCTAGGCGTACGGCAGGGAAGGGACAGGAGATCCCGTCCGCTACTATTACCCTCGACGCCAAAGGTAAGGAACAGCAGGGGAATCGCTACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCAt                                                                                                                         <  1:225287/251‑1 (MQ=255)
                                                                                                                        aGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGt                                                                                                                     <  1:444276/251‑1 (MQ=255)
                                                                                                                                 gtccGTTCACTGGTGAGAAGCAGTATAAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACACTCGACGCCAAAGGGAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTc                                                                                                                                          >  1:296160/3‑221 (MQ=255)
                                                                                                                                   ccGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTg                                                                                                           >  1:156980/1‑250 (MQ=255)
                                                                                                                                        cACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAg                                                                                                                                                                                                                       <  2:198894/137‑1 (MQ=255)
                                                                                                                                        cACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAg                                                                                                                                                                                                                       >  1:198894/1‑137 (MQ=255)
                                                                                                                                        cACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTa                                                                                                      >  1:440922/1‑250 (MQ=255)
                                                                                                                                         aCTGGTGAGAAGCAGTATCCGGCTGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAAc                                                                                                    <  2:122028/251‑1 (MQ=255)
                                                                                                                                         aCTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCGCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTTAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAAc                                                                                                    <  1:461841/251‑1 (MQ=255)
                                                                                                                                            ggTGAGAAGCAGTATCAGGCGTAGGATCGTAAACGGAAACTTACCGGCAAAATGCCCGATTACACGTACTCTATGCGTACGGCAGTGAAGGGACAGGAGAACCCGCCCTCGACTATTACCCTCGACGCCAAAGGTAACGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGTAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGTCCATTAGTTGGCGCGGTGAGCTTACGCt                                                                                                 <  2:94150/251‑1 (MQ=255)
                                                                                                                                            ggTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGc                                                                                                  <  1:264629/250‑1 (MQ=255)
                                                                                                                                               gagaAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTTCCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACGCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGTGGTGAGCTAACGCTTaa                                                                                              <  1:428967/251‑1 (MQ=255)
                                                                                                                                               gagaAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGACATTAGTTGGCGCGGTGAGCTAACGCTTa                                                                                               >  1:433328/1‑250 (MQ=255)
                                                                                                                                                   aGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGc                                                                                          <  1:329948/251‑1 (MQ=255)
                                                                                                                                                       gTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTa                                                                                      <  1:102626/251‑1 (MQ=255)
                                                                                                                                                       gTATCAGGCGGATGATCTGAAACTGAAACTTACCAGCAAAATGACCGATTACACGCCTTCTATGCGTATGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTa                                                                                      <  2:289846/251‑1 (MQ=255)
                                                                                                                                                         agcAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCCCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTaa                                                                                     <  1:169100/248‑1 (MQ=255)
                                                                                                                                                              gCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGc                                                                               >  1:375871/1‑251 (MQ=255)
                                                                                                                                                                     atCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGCGAGCTAACGCTTAACGGCATTAACACCGCCAAag                                                                          >  1:216344/1‑249 (MQ=255)
                                                                                                                                                                      tCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACCCAGCCAAagagg                                                                       >  1:246192/1‑250 (MQ=255)
                                                                                                                                                                         gaaacttaaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGGACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACGATTCCCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGTTGAGCTAACGCTTAACGGCATTAACACCGCCAAagag                                                                        <  2:375871/247‑1 (MQ=255)
                                                                                                                                                                         gaaactgaaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTc                                                                     <  1:407262/250‑1 (MQ=255)
                                                                                                                                                                         gaaactgaaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTc                                                                     <  2:440922/250‑1 (MQ=255)
                                                                                                                                                                                aatcttcccgaaaaaaTGACCGATTACAAGCTATCTATGCTTACGGCAGTTAAGGGACAGGAGATCCCGACCGCGATTATTACCCTCGACGCCAAAGGTAAGGAACAGCAGGTCAATCTCGACAAACTCTCCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGTCAACAGGCCATTATTTGGCGCGGTGATCTAACGCTTAACGGCATTAACTCCGCCAAAGAGTTCCCGGACTg                                                             <  2:224197/239‑1 (MQ=255)
                                                                                                                                                                                      aCCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGt                                                        <  1:129301/250‑1 (MQ=255)
                                                                                                                                                                                        cGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCg                                                      >  1:454948/1‑250 (MQ=255)
                                                                                                                                                                                              aaTGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTc                                                >  1:144135/1‑250 (MQ=255)
                                                                                                                                                                                                  aCCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGCGTTCCCGGACTGGCCGTCGAAACTCAATgg                                           >  1:172753/1‑251 (MQ=255)
                                                                                                                                                                                                  aCCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGACATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGCGTTCCCGGACTGGCCGTCGATACTCAATgg                                           >  1:451778/1‑251 (MQ=255)
                                                                                                                                                                                                    cGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCt                                         >  2:132503/1‑251 (MQ=255)
                                                                                                                                                                                                         acacGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCCCCGTCGCGGCGCTGGAAGGGAAAACTGACCTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGAtt                                    >  2:124308/1‑251 (MQ=255)
                                                                                                                                                                                                         acacGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGAtt                                    >  2:125918/1‑251 (MQ=255)
                                                                                                                                                                                                              ctctctATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAAc                               >  1:78096/1‑251 (MQ=255)
                                                                                                                                                                                                               tctctATGCGTACGGCAGTGACGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGCCGCCAAAGGTAATGAACAGCAGGTCAATCTCGACCAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCCACGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAAcc                              >  2:245955/1‑251 (MQ=255)
                                                                                                                                                                                                                    gggctgccGGCAGGGAAGGGACGGGAGACCCCGTCCGAGACTATACCCCTCGACGCCAAAGGTAAAGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCGGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCgg                         <  2:475300/244‑1 (MQ=255)
                                                                                                                                                                                                                      gCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACGCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGt                        <  1:403361/250‑1 (MQ=255)
                                                                                                                                                                                                                        gTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGc                     >  1:104371/1‑251 (MQ=255)
                                                                                                                                                                                                                          aCGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAATGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGGGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGGCCAAACTCAATGGCTTTATTAAAACCCGCgggaacat                   >  2:407573/1‑245 (MQ=255)
                                                                                                                                                                                                                               aGTGAAGGGACAGGAGATCCCGTCCGCGACTATTCCCCGCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACCCAAGGCGTTGCTCGACTGGCAACAGGCAATTAGTGGGCGCGGGGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTATCCTGTAcg              <  2:458147/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAAACCGCGGGAGCCTGTAcggcg           >  2:22145/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACggcggc         >  1:394101/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gACAGTAGATCCCGTCCGCGACTATTACCCCCGACGCCAAAGGTAATGAACAGCAGGTCAACCTCGACAAACCCACCGCCGCGGCGCTGGAAGGGAAAACTGACCACAAGGCGTTGCGCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCAGTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACACAATGGCTTGCTTAAAACCCGCGGTAGCCTGTACGGCTGCAcc      <  2:205643/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTACCGCTTAACGGCATTAACACCGCCAAATAGTTCCCGGTCTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGAAGGCTGTAAggcggccac      >  2:109107/1‑248 (MQ=255)
                                                                                                                                                                                                                                       gACAGGAGAGCCCGTCCGCGACTATTCACCTCGACTCAAAAGGTATGAACAGCAGGTCAAGCACGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGACCGCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATGAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCTGTAGCCTGTACGTCGGCACCt     <  2:334128/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gaTCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGc  >  1:239439/1‑248 (MQ=255)
                                                                                                                                                                                                                                              gaTCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCCGGc  >  2:185371/1‑248 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
CGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGC  >  NC_000913/4445070‑4445555

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