Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,697,170 C→T A238V (GCG→GTG)  bcsE → c‑di‑GMP‑binding protein BcsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,697,1700CT91.7% 133.1 / ‑5.1 49A238V (GCG→GTG) bcsEc‑di‑GMP‑binding protein BcsE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/18);  minor base G (4/0);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

AACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGCGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTGGCCTGCGGTGCA  >  NC_000913/3696924‑3697411
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
aaCAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGtgccg                                                                                                                                                                                                                                               >  1:223432/1‑251 (MQ=255)
   aaTTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGANGAATGTTGCTGTACTGGAAGGTGtgccgccg                                                                                                                                                                                                                                            >  2:213914/1‑251 (MQ=255)
         cATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCg                                                                                                                                                                                                                                      >  1:17630/1‑251 (MQ=255)
          aTTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAGACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGGGCCGCCGCTATCgg                                                                                                                                                                                                                                     >  2:278242/1‑251 (MQ=255)
           ttgttgCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGGGTGCCGCCGCTATCgg                                                                                                                                                                                                                                     >  1:38844/1‑250 (MQ=255)
            tgttgCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGCCCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGCAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGGGCCGCCGCTATCGGa                                                                                                                                                                                                                                    >  2:257758/1‑250 (MQ=255)
            tgttgCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGCAAAAGGGGTCAGCGCCCGTCAGCAGATTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAACCGCATTCTGAGTAATGTTGCTGTACTGGAATGTGTGCCGCCGCTATCGGa                                                                                                                                                                                                                                    >  1:412981/1‑250 (MQ=255)
            tgttgCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGGGTGCCGCCGCTATCGGa                                                                                                                                                                                                                                    >  2:104900/1‑250 (MQ=255)
                  ttGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTgg                                                                                                                                                                                                                             >  1:98274/1‑251 (MQ=255)
                           aCCGTTCCCTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAGACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTc                                                                                                                                                                                                                    <  1:278242/251‑1 (MQ=255)
                                                ccAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTc                                                                                                                                                                                               <  1:185159/251‑1 (MQ=255)
                                                                                                 gtgcagctaaaaaGGGGTCACGCTCCGTCATCAACTTATCGTTCTGCAACAAAATGTTATCTGGACAGTAATTCAAAGCGAAGAGCCAGAGATCCAACCCCGCAGCGACGAAAAACGCGTTCTTATTAATGTTGCTGTACTGGAAGTTGTGCCGCCGCTATCGTAACACGGGCAACTGTTCAACAATAACGAAGTCCTTTTCAATGAAGCCCGTACCGCTCAGGCGGCGACTGTGGTCTTGTCTTTACatc                                                                                                                                              <  2:310892/243‑3 (MQ=255)
                                                                                                      aCGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTTGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTTAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAAt                                                                                                                                         <  1:130512/251‑1 (MQ=255)
                                                                                                      aCGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAAt                                                                                                                                         >  1:235195/1‑251 (MQ=255)
                                                                                                              ggggTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAAGGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAAt                                                                                                                                 <  2:410554/251‑1 (MQ=255)
                                                                                                                   cAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGc                                                                                                                            >  2:141365/1‑251 (MQ=255)
                                                                                                                    aGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGcc                                                                                                                           >  1:325515/1‑251 (MQ=255)
                                                                                                                            gTCAGCAGCTTAGCGATCAGCAACAAAATGGAATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGAGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGACGCCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCCCAAATCGAGCCCCTGGcc                                                                                                                    >  2:54700/1‑251 (MQ=255)
                                                                                                                              cagcagCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGCGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCg                                                                                                                  >  1:224848/1‑251 (MQ=255)
                                                                                                                                   gcTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGCCGAAAAACGCATTCTGAGCAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGATCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCANAATGCGCAAATCGAGCCACTGGCCCgcagc                                                                                                              >  2:216042/1‑250 (MQ=255)
                                                                                                                                             cAGCAACAAAATGGTATCTGGACATTAGTTCAAAGGGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATAccc                                                                                                    <  2:349158/250‑1 (MQ=255)
                                                                                                                                                  acaaAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCCTTCATACCCTGCgcc                                                                                              >  2:13152/1‑249 (MQ=255)
                                                                                                                                                  aaaaaaTGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTc                                                                                              <  2:98274/249‑1 (MQ=255)
                                                                                                                                                  aaaaaaTGGTACTGGGACATTAGTTCAAAGGGAAGAGGAGGAGACCCACCCACGCAGCGTCGAAAACCTCATTCTGATTAATGTTGCTGTACTGGATGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCTAGCCACTGGCCCGCAGCATTCATACCCTGCGTc                                                                                              <  2:311530/249‑1 (MQ=255)
                                                                                                                                                   caaAATGGAATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGGGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCg                                                                                             >  1:195296/1‑251 (MQ=255)
                                                                                                                                                     aaagggTATCTGGACATTAGTTAAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCGGTACTGGAAGGTGTGCCGCCGTTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGcc                                                                                           <  1:54700/247‑1 (MQ=255)
                                                                                                                                                     aaaTGGTATCTGGACATTAGTTCAAAGCGAANGGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGcc                                                                                           <  2:235195/251‑1 (MQ=255)
                                                                                                                                                     aaaTGGTATCTGGACAGTAGTTCAAAGCGAAGAGGCGGAGATCCACCCACGCAGGGACGAAAAACGCATTCTGATTAATGTTGCTGAACTGGAAGGGGTGCCGCCGCTATCGGAACACTGGCAACTGTCCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGcc                                                                                           <  2:130661/251‑1 (MQ=255)
                                                                                                                                                         ggTAACTGGACATTAGTTCAAAGCGAATAGGCGGATATCCAACCACGCAGCGACTAAAAACGCTTTCGGATTAATGTTGCTGTACGGGAAGGTGTGCCGCCGCTATCGGACCACTGGCAACTTTTCAACAATAACGAAGTCCTGTTCAATGAGGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTGCTTTACAGCAAATTGCGCAAATCGAGCCACTGGCCCGCAGCATTCATTCCCTGCGTCGCCAgc                                                                                        <  2:257636/250‑1 (MQ=255)
                                                                                                                                                                         gTTCAAAGCGAAGAGGCGGAGATCCAAACACGCAGCGACGAAAAAAGCATTCTGAGTAATGTTGCTGTACAGGAAGGTGAGCCGCCGCTATCGGAACACTGGCAACTGTTCCACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTAAAGCAAAATGGGCAAATCGAACCACTGGCTCGCAGCATTCATACCCTGCGTCTCCAGCGCGGAAGAGCGATGa                                                                         >  2:48246/1‑250 (MQ=255)
                                                                                                                                                                            cAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAAt                                                                     >  1:28222/1‑251 (MQ=255)
                                                                                                                                                                                gCGAAGAGGCGGAGATCCAACCCCGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTc                                                                 <  1:13152/251‑1 (MQ=255)
                                                                                                                                                                                 cGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCg                                                                >  1:409352/1‑251 (MQ=255)
                                                                                                                                                                                       ggCGGAGACCCAACCACGCAGCGACGAAAAACGCATTCGGAGTAATGTGGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCggg                                                          <  2:17630/251‑1 (MQ=255)
                                                                                                                                                                                          ggAGATCCAACAACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCCGTTCAATGAAGCACGTAACGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGAGGGAGTGCGATGAAAATCCTCTTGCGGGaa                                                        >  2:39368/1‑250 (MQ=255)
                                                                                                                                                                                             gaTCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAAt                                                     >  1:459827/1‑250 (MQ=255)
                                                                                                                                                                                                       cGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCt                                          >  1:453421/1‑251 (MQ=255)
                                                                                                                                                                                                          aGCGACGAAAAAAGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAACTACCGCTAGCCTgcg                                       >  2:31355/1‑251 (MQ=255)
                                                                                                                                                                                                             gacgaAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTgcgcgc                                    >  1:438981/1‑251 (MQ=255)
                                                                                                                                                                                                              acgaAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGcc                                   >  1:128032/1‑251 (MQ=255)
                                                                                                                                                                                                              acgaAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGcc                                   >  2:158538/1‑251 (MQ=255)
                                                                                                                                                                                                                              gAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGttt                     <  1:141365/249‑1 (MQ=255)
                                                                                                                                                                                                                              gAGTAATGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGt                   >  1:111867/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaTGTTGCTGTACGGGACGTTGTGCCGCCGCNATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTAtt               <  2:223432/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGTTGCTGTACTGAAAGTTGGGCAGTCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGCCGACGACGGAGGTATTTTCTTTACAGCCAAATGCGCAATTCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGTTGAAAACCCTCTTGCGGGGAACTACCTCTGTCCCGCGCCCCGCCGCTGAACATTTGTTAtgc              >  2:275564/1‑249 (MQ=255)
                                                                                                                                                                                                                                    tGTTGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGCAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAACTACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTgg             >  1:227175/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tGCTGTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGCGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGCAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTGGCCt          >  2:332488/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cTGTACTGGAAGGAGTGCCGCCGCTATCGGAACACTGGAAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTCGTGCGATGAAAATCCTCGTGCGGGCAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTGGCCTGc        >  1:131328/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gTACTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTGGCCTGCgg      <  1:213914/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aCTGGAAGGTGTGCCGCCGCTATCGGAACACTGGCAACGGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTGGCCTGCGGt     <  1:85044/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ggAAGGTGTGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTGGCCTGCGGTGCa  >  1:89019/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
AACAATTTTCATTGTTGCTTGAGGAATACCGTTCACTTTTTGGTCTTGCCAGTTTGCGTTTTCAGGGTGACCAACATTTGCTGGATATTGCCTTCTGGTGCAACGAAAAAGGGGTCAGCGCCCGTCAGCAGCTTAGCGTTCAGCAACAAAATGGTATCTGGACATTAGTTCAAAGCGAAGAGGCGGAGATCCAACCACGCAGCGACGAAAAACGCATTCTGAGTAATGTTGCTGTACTGGAAGGTGCGCCGCCGCTATCGGAACACTGGCAACTGTTCAACAATAACGAAGTCCTGTTCAATGAAGCCCGTACCGCTCAGGCGGCGACGGTGGTCTTTTCTTTACAGCAAAATGCGCAAATCGAGCCACTGGCCCGCAGCATTCATACCCTGCGTCGCCAGCGCGGTAGTGCGATGAAAATCCTCGTGCGGGAAAATACCGCTAGCCTGCGCGCCACCGATGAACGTTTGTTATTGGCCTGCGGTGCA  >  NC_000913/3696924‑3697411

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