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,3150CT93.5% 169.8 / ‑0.9 63P402S (CCC→TCC) tamBtranslocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (28/30);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGG  >  NC_000913/4445074‑4445561
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
gcgcGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACt                                                                                                                                                                                                                                                <  1:37013/250‑1 (MQ=255)
 ggcggaCAGCAGGAGATGGGCGTTAATCATTCCGGTCCGGTGTATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGGGAACAGCAAACATCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACGTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTGCGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTAt                                                                                                                                                                                                                                              <  2:148803/246‑1 (MQ=255)
       cAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGATCTCTATGCGTACGGCTGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCt                                                                                                                                                                                                                                        >  1:72731/1‑251 (MQ=255)
           tGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGAc                                                                                                                                                                                                                                    >  1:222278/1‑251 (MQ=255)
              agcttGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGCCGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTCCTGGCCGTTCACTGGTTAGAAGCAGTAACAGGCGGATGATCTGAAACTGCAACTTACCAGCAAAATGACAGATTACACGCTCTCTGTGCGTACGGCAGTGGAGGGACAGGGGATCCCCTCCGCGACCATTACCCTCGTCGcc                                                                                                                                                                                                                                 >  2:130402/4‑251 (MQ=255)
                   ggCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAg                                                                                                                                                                                                                             <  1:173610/250‑1 (MQ=255)
                       ttAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTAATGGCCGTTCACTGGTGAGAAACAGTATCAGGCGGATGATCTGAAACTGAAACTTAACGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAAt                                                                                                                                                                                                                        >  1:294913/1‑251 (MQ=255)
                       gtAATCTTACCGGTCCGTTGGATATGGATTTACGCGCCCAGACGCGAATGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTTAACAGCAAACAGCTTTACTGGCCGTTCACTGTTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACCTACCGGCAAAATGACCGATTACTCGCCCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGCCTATTTCCCTCGACGCCAAAGGTAAt                                                                                                                                                                                                                        <  2:74446/250‑1 (MQ=255)
                          aTCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGa                                                                                                                                                                                                                      >  1:433588/1‑250 (MQ=255)
                                   gtccggtGGATATGGATTTACGCGTCCAGACGCGACTGGCGGAAGCCGGATTGCCTCTCAACGTGGAAGTGAACAGCAACCAGCTTTTCTGGCCGTTCACTGGTGAGAAGCAGTCACAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGACCCCGTCCGCGACTATTACCCTCGACGCCAAAGTTAATGAACAGCAGGTc                                                                                                                                                                                                            <  2:257574/251‑1 (MQ=255)
                                     ccgggGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCa                                                                                                                                                                                                           <  2:122556/250‑1 (MQ=255)
                                            atatGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGAc                                                                                                                                                                                                   <  1:334609/251‑1 (MQ=255)
                                                ggATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACaaa                                                                                                                                                                                                <  1:43681/250‑1 (MQ=255)
                                                     tccgctcccAGACGCGACTGGCGGAAGCCGGATTGCCGCTAAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACTCGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCAc                                                                                                                                                                                           <  2:406708/244‑1 (MQ=255)
                                                            ccAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAACCCGCTTTACTGGCCGTTCCCTGGTGAGAAGCAGTATCAGGCGGATGATCTGCAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATCCGTACGGAAGTGAAGGGACAGGAGATACCGTCCGCCACTCTTACCCTCGACGCCAAAGGTAATGAACAGCCGGCCCATCTCAACAAAATCCACGTCgcgg                                                                                                                                                                                   >  1:483048/1‑251 (MQ=255)
                                                              aGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCgcggcg                                                                                                                                                                                 <  2:229352/251‑1 (MQ=255)
                                                                     aCTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGTGTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAg                                                                                                                                                                          <  2:201202/251‑1 (MQ=255)
                                                                      cTGGCAGAAGCCGGATTGCCGCTCAACGTGGCAGTGACCAGCAAACAGCTTTACTGGCCGTTCACTGGTGTGAAGCAGTATCAGGCGGATGATATGAAACTGGAACTTAACGGCAAAATGACCGATTACACGCTCTCTATGCGTCCGGCCGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGGCGCCAAGGGTACTGCACAGCAGGTCAATCTCGACCAACTCACCGTCGCGGCGGTGGCAgg                                                                                                                                                                         >  2:207799/1‑251 (MQ=255)
                                                                       tGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAggg                                                                                                                                                                        >  1:393204/1‑251 (MQ=255)
                                                                        ggCGGAAGCCGGATTGCCGCGCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGTATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCTTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGGAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGa                                                                                                                                                                       <  1:39772/251‑1 (MQ=255)
                                                                                 cGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTTCCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACCCACCGTCGCGGCGCTGGAAGGGAAAACTGAAc                                                                                                                                                              <  1:167573/251‑1 (MQ=255)
                                                                                      tGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCaa                                                                                                                                                          <  2:222278/250‑1 (MQ=255)
                                                                                               aCGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTc                                                                                                                                                <  1:483541/251‑1 (MQ=255)
                                                                                                cGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCCGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGCAAACTGAACTCAAAGCGTTGCTCg                                                                                                                                               >  2:128689/1‑251 (MQ=255)
                                                                                                  gggAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCTCCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGAc                                                                                                                                             <  1:81437/250‑1 (MQ=255)
                                                                                                   ggAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACt                                                                                                                                            <  1:391884/251‑1 (MQ=255)
                                                                                                           aCAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAAc                                                                                                                                      >  1:424923/1‑249 (MQ=255)
                                                                                                              gCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGcc                                                                                                                                 <  1:173549/251‑1 (MQ=255)
                                                                                                                   cAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAg                                                                                                                            <  2:424923/251‑1 (MQ=255)
                                                                                                                        ttACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGgc                                                                                                                       >  1:28798/1‑251 (MQ=255)
                                                                                                                            tGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCgg                                                                                                                   >  2:409851/1‑251 (MQ=255)
                                                                                                                               ccGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCACCAGGCCATTAGTTGGCGCGGTg                                                                                                                 >  1:478238/1‑250 (MQ=255)
                                                                                                                                 gTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGACCATGAGTTGGGTCGGTGAGc                                                                                                              >  2:169504/1‑251 (MQ=255)
                                                                                                                                 gTTCACTGGTGAGAAGCAGTATCAGGCAGATGAGTTGAAACTGAAACTGCCCGGCCAAATGACCTATTACACTCTCTCTATGCGTACGGCAGTGAATGGACAGGAGATCCCCGCCGCGACTATTCCCCTCGACGCCAAAGGGAATGAACAGCAGGCCAATCTCGACAAACTCACCGTCGCGGCGCTTGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATGAGGTGGCGCGGTGAGc                                                                                                              <  2:38335/251‑1 (MQ=255)
                                                                                                                                            agaAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAAc                                                                                                   <  2:35722/251‑1 (MQ=255)
                                                                                                                                              aaGCAGTATCAGGCGGACTATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGATATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCGCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACgg                                                                                                 <  1:210912/251‑1 (MQ=255)
                                                                                                                                               aGCAGTATCAGTCGGTGGATCTGAAACCGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGc                                                                                                <  2:379130/251‑1 (MQ=255)
                                                                                                                                               aGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGc                                                                                                <  1:375350/251‑1 (MQ=255)
                                                                                                                                                 cAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCAt                                                                                              <  1:250232/251‑1 (MQ=255)
                                                                                                                                                       cAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTCCACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAccac                                                                                        <  2:273489/251‑1 (MQ=255)
                                                                                                                                                       cAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAaca                                                                                         <  1:119610/250‑1 (MQ=255)
                                                                                                                                                                gatCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCCACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCaa                                                                                  >  1:215878/1‑248 (MQ=255)
                                                                                                                                                                 atCTGAACCTGAAACTTACCGACAAAATGACCGATTACCCTCTCTATATGCGTACTGCAGTGAAGGGTCAGTAGATCCCGTCCGCGACTATTACCCCCGACGGCAAAGGTAAGGACCAGCAGGGCAATCTGGGCAATCTCCCGGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGCCGTTGCTTGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAaga                                                                               <  1:61273/250‑1 (MQ=255)
                                                                                                                                                                         ctgaaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCATCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGACAATCTCGCCAAACTCACCGTCGCGGAGATGGAAGGGAACACTGAACGCAAGGCGTTGCTCGACTGGCCACGGGCAATTACTTGGCGGGGGGAGCTACCGCTGAACGGCCTTAACACCGCCAAACAGGTCCCgg                                                                       >  2:56039/1‑250 (MQ=255)
                                                                                                                                                                           gaaactTACCGGCAAAATGACGTATTAACCGCTCCCTATGCGAACGGTGGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACCCACCGTCGCGGCGCTGGAAGGGAAAACTTAACTCATGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTCCTCGGAc                                                                     <  2:41770/250‑1 (MQ=255)
                                                                                                                                                                             aactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTgg                                                                  >  1:167791/1‑251 (MQ=255)
                                                                                                                                                                                        aaaaTGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCg                                                                                                                                                      >  1:39953/1‑156 (MQ=255)
                                                                                                                                                                                        aaaaTGACCGATGACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCg                                                                                                                                                      <  2:39953/156‑1 (MQ=255)
                                                                                                                                                                                          aaTGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGACGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTcc                                                     >  2:93361/1‑250 (MQ=255)
                                                                                                                                                                                          aaTGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCAGTCCGCGACTATTAACCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACCGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCCTTAACCACGCCAAAGAGATCCCGGACTGGCCGTCGAAACTc                                                      >  2:76508/1‑250 (MQ=255)
                                                                                                                                                                                                  aTTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGCAATGAACAGCAGGTCAATCTCGACACACTCCCCGTCGCGTCGCTGGAAGGGAAAACTGACCTCAAGTCGTTGCTCGACTGGCAACAGGCTATTAGTTGGCTCTGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTg                                             <  1:207799/251‑1 (MQ=255)
                                                                                                                                                                                                      cacGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTa                                         <  1:128689/251‑1 (MQ=255)
                                                                                                                                                                                                        cGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTaa                                        >  1:208986/1‑250 (MQ=255)
                                                                                                                                                                                                         gctcgttaggcgcaCGGCAGGGAAGGGACAGGAGAACCCGTCCGCGAGGCTTGCCCACGATGCTAAAGGTACTCAACAGCAGGTCAATCGCGGCAAACTCCCCGTCGCGTCGCTGGAAGGGAAACCTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGTAGTAACCCCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCGTGATTaaaa                                      <  2:115669/238‑1 (MQ=255)
                                                                                                                                                                                                           tctctATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAAcc                                    >  1:229039/1‑251 (MQ=255)
                                                                                                                                                                                                                   cGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGCCTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGCTTAAAACCCGCGGt                              >  2:207804/1‑249 (MQ=255)
                                                                                                                                                                                                                   cGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGt                              >  1:78914/1‑249 (MQ=255)
                                                                                                                                                                                                                    gTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCCAAGGTAATGAACAACAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGACGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCCTTCACACCGCCAAAGAGTTCCCGGACTGTCCGGCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGc                           >  2:505097/1‑251 (MQ=255)
                                                                                                                                                                                                                     gccgccggTGAAGGGACAGGAGATCCCTCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCCCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGACCTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGTTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGcc                          <  2:363362/244‑1 (MQ=255)
                                                                                                                                                                                                                       cGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAAAACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTaaaaacccaggaagcctg                        >  2:498278/1‑237 (MQ=255)
                                                                                                                                                                                                                        ggCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGt                       >  1:411591/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gggACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACggcggc               >  2:118622/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGATTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGACTGGACAGAGGCACCTg          >  2:405030/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cccgcccgAGACTATTACCCGCGCCGCCTAAGGTAATGAACAGCGGGTCAATCTCGACAAACTCACCGGCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTGGCTCGACGGGCAACAGGCCATTAGTGGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATgg  <  2:473057/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
GCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATGG  >  NC_000913/4445074‑4445561

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