Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,953,879 C→T H146H (CAC→CAT cmoB → tRNA U34 carboxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,953,8790CT96.7% 186.2 / NA 60H146H (CAC→CATcmoBtRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/26);  minor base G (0/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

AGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCAACTGGTCAATGAAGGTGAACTGGTGCTGGAAACGCTGG  >  NC_000913/1953640‑1954123
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
aGCGAAGAGCCACTTAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATTCCGAGTTGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGAGTTAGCCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTAGCACATGGGGCGCATGATTTGCGCTGGCGCGCAGCTCGCGGTggg                                                                                                                                                                                                                                            <  2:413973/251‑1 (MQ=255)
aGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGAGCGCATCTCGCGGTggg                                                                                                                                                                                                                                            >  2:542892/1‑251 (MQ=255)
aGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGATGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGATCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTggg                                                                                                                                                                                                                                            >  1:609277/1‑251 (MQ=255)
    aaGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGTAGGCGCGCATCTCGCGGTGGGGATc                                                                                                                                                                                                                                        <  2:309064/251‑1 (MQ=255)
      gagCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGa                                                                                                                                                                                                                                      <  1:199017/251‑1 (MQ=255)
          cACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCt                                                                                                                                                                                                                                  >  2:91568/1‑251 (MQ=255)
                gcgcAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCa                                                                                                                                                                                                                             >  1:409633/1‑250 (MQ=255)
                         aaaTTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCtg                                                                                                                                                                                                                   <  1:123774/251‑1 (MQ=255)
                               aaCNCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAg                                                                                                                                                                                                             <  2:596242/251‑1 (MQ=255)
                                  gCATGGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGttt                                                                                                                                                                                                          <  2:244786/251‑1 (MQ=255)
                                   cATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTg                                                                                                                                                                                                         >  2:240361/1‑251 (MQ=255)
                                    aTTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTg                                                                                                                                                                                                         >  2:341113/1‑250 (MQ=255)
                                       gAAACGCTGATGCGCAACCTGGTGCCGTGGCGCAGAGGGCCGTTCTCACTGTATGGCGTCTACATCGATACCGAATGGCGTTCCGACTGGAAAGGGGATCGTGTTCTGCCCCATCTGTCTTATTTAACCGGACGCACCATTCTGTTTGTCTGCTGCTGCAGCGGTTATGTCATGTGGCTCATGATTGGCGCAGGCGCGGATCTCGGGGTGGGTGTCGACCCTACGCAACTGTTTCTCTGCCAGTTTGAAGc                                                                                                                                                                                                     <  2:634790/251‑1 (MQ=255)
                                        aaaCGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCg                                                                                                                                                                                                    >  1:614938/1‑251 (MQ=255)
                                              tGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGt                                                                                                                                                                                              >  2:286522/1‑250 (MQ=255)
                                              tGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCATCAACATCGATACCGATTGGCGTTCCGACTGGAAATGGGATCGTGTTATGCCCCCTCTTTCTGATTTAACCGGAAGCACCATTCTGTATGTCGAATGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTTGCGGTGGAGATCGACCCTACGCAACTGTTTCTGTGCCCGTTTGAAgggttgagt                                                                                                                                                                                              >  2:161189/1‑243 (MQ=255)
                                                  gcgcAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAac                                                                                                                                                                                          >  1:197626/1‑251 (MQ=255)
                                                    gcAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCTCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTGTGAAGCGGTGCGTAAacta                                                                                                                                                                                        <  1:53974/251‑1 (MQ=255)
                                                      aaCCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAactact                                                                                                                                                                                      <  1:84538/251‑1 (MQ=255)
                                                        ccTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTgg                                                                                                                                                                                    >  1:247530/1‑251 (MQ=255)
                                                                cGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGAt                                                                                                                                                                           >  1:572899/1‑251 (MQ=255)
                                                                   ggCGCAAAGGGCCGTTCTCACTGTCTGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCCACTGTTTCTGTGCCAGTTTGAAGCGGCGCGTAAAATCTGGGGTAATGATc                                                                                                                                                                          >  2:626970/1‑249 (MQ=255)
                                                                      ccAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGTGATCGTGTTCTGCCCCATCTTGCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCGCATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACGGGGTAATGATCAgcgc                                                                                                                                                                     <  1:573693/250‑1 (MQ=255)
                                                                         aaGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAgcgcgcg                                                                                                                                                                  <  1:331272/251‑1 (MQ=255)
                                                                             gCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCAtt                                                                                                                                                              >  2:177519/1‑251 (MQ=255)
                                                                               cGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGAGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTg                                                                                                                                                            >  2:24817/1‑251 (MQ=255)
                                                                                 ttCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGtt                                                                                                                                                          >  1:376225/1‑251 (MQ=255)
                                                                                               cGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTg                                                                                                                                            <  1:230753/251‑1 (MQ=255)
                                                                                                  cAACATCGATACCGAATTGCGTTCCGACCGCAAATGGGATCGTGTTCTGACCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCCGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGCGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGaac                                                                                                                                         >  2:152119/1‑251 (MQ=255)
                                                                                                   aaCATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGaac                                                                                                                                         >  2:186864/1‑250 (MQ=255)
                                                                                                          aTACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCg                                                                                                                                 <  2:73725/251‑1 (MQ=255)
                                                                                                           tACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACGCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCgg                                                                                                                                >  1:585651/1‑251 (MQ=255)
                                                                                                                 aTGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTg                                                                                                                           >  1:298479/1‑250 (MQ=255)
                                                                                                                 aTGGCGTTACGACTTGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCCGACGAACCATTCTGGAGGTCGGTTGCGGCAGCGATTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCTAGTTTGAAGCGGTGCGTAATCTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACATCTTCCGTCACtta                                                                                                                          <  2:531194/251‑3 (MQ=255)
                                                                                                                       ttCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGACCAACTTCCGGCACTGAAAGCCt                                                                                                                    >  2:42174/1‑251 (MQ=255)
                                                                                                                               ggAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTACCCGGACGCACCATTATGGCTGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATAcc                                                                                                            >  2:259181/1‑251 (MQ=255)
                                                                                                                                     gggATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCttt                                                                                                      <  2:879/251‑1 (MQ=255)
                                                                                                                                           gtgtTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCTGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCAGCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCg                                                                                                                                                                                               <  1:657534/157‑1 (MQ=255)
                                                                                                                                           gtgtTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCg                                                                                                                                                                                               >  2:657534/1‑157 (MQ=255)
                                                                                                                                             gtTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATgg                                                                                               <  1:564710/250‑1 (MQ=255)
                                                                                                                                                  gCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCTTCTTTTCGATGGGCg                                                                                            <  1:566069/248‑1 (MQ=255)
                                                                                                                                                         cTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTAcc                                                                                   >  1:349632/1‑250 (MQ=255)
                                                                                                                                                           ttCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCa                                                                                  <  1:91568/249‑1 (MQ=255)
                                                                                                                                                              tGATTTAACAAGACGCACCATACTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGCGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGCAAGACTTTGATATCGTCTTTTCGATGGGCGTGCTTTACCATCGcc                                                                             >  2:626092/1‑251 (MQ=255)
                                                                                                                                                               gATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCg                                                                               <  2:409633/248‑1 (MQ=255)
                                                                                                                                                                 tttAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGtt                                                                          >  1:582174/1‑251 (MQ=255)
                                                                                                                                                                  ttAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGtt                                                                          <  1:344921/250‑1 (MQ=255)
                                                                                                                                                                   tAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTc                                                                         <  1:288868/250‑1 (MQ=255)
                                                                                                                                                                           cGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTGCCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGa                                                                <  2:376225/251‑1 (MQ=255)
                                                                                                                                                                                 aTTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCt                                                          >  1:499648/1‑251 (MQ=255)
                                                                                                                                                                                          gTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGtt                                                 <  1:32212/251‑1 (MQ=255)
                                                                                                                                                                                              gCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGaaa                                             <  1:42174/251‑1 (MQ=255)
                                                                                                                                                                                                 gCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCTCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGTCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAg                                            <  1:626970/249‑1 (MQ=255)
                                                                                                                                                                                                    gCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCaa                                       <  2:392052/251‑1 (MQ=255)
                                                                                                                                                                                                                        cgcATGATTGGCGCAGGCGCGCATCTCGCGGTGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCAACTGGTCAATGAAGGTGAACTg                  <  1:286522/251‑1 (MQ=255)
                                                                                                                                                                                                                           aTGATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAAATACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCAACTGGTCAATGAAGGTGCACTggg                >  2:16356/1‑250 (MQ=255)
                                                                                                                                                                                                                             gATTGGCGCAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCAACTGGTCAATGAAGGTGAACTGGTGc              >  1:339823/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gcgcAGGCGCGCATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCAACTGGTCAATGAAGGTGAACTGGTGCTGGa          >  1:139918/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gcAGGCGCGCATCTCGCGGGGGGGTTCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTCAACTACTGGGTAATGATCAGCGCGCGCATTTGTTCCCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCATGTGAAAGCTCACCTGGTCAATGACGCTCAACTGGGGCTggcacc       >  2:201219/1‑247 (MQ=255)
                                                                                                                                                                                                                                         cgcgcATCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCAACTGGTCAATGAAGGTGAACTGGTGCTGGAAACGCTgg  >  1:535773/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
AGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCCCCATCTTTCTGATTTAACCGGACGCACCATTCTGGATGTCGGCTGCGGCAGCGGTTATCACATGTGGCGCATGATTGGCGCAGGCGCGCACCTCGCGGTGGGGATCGACCCTACGCAACTGTTTCTGTGCCAGTTTGAAGCGGTGCGTAAACTACTGGGTAATGATCAGCGCGCGCATTTGTTACCGTTAGGTATTGAACAACTTCCGGCACTGAAAGCCTTTGATACCGTCTTTTCGATGGGCGTGCTTTACCATCGCCGTTCACCGCTGGAGCATCTATGGCAGTTGAAAGATCAACTGGTCAATGAAGGTGAACTGGTGCTGGAAACGCTGG  >  NC_000913/1953640‑1954123

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