Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,024 A→G E177E (GAA→GAG slmA → nucleoid occlusion factor SlmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,0240AG98.3% 189.0 / NA 58E177E (GAA→GAGslmAnucleoid occlusion factor SlmA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (33/24);  minor base T (0/1);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGCCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCT  >  NC_000913/3814782‑3815266
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
ggTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGAGCTAATGCTTGAACAGGATCGCCTGCTAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCATGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCCCCGATGAAACCATGCTGGCAAGCCAGATCCTGGCCTTCTGTGACGGTATGCTGTCAAGGCTTGTCGGCAGCGAGTTTAAAt                                                                                                                                                                                                                                               >  2:349945/1‑250 (MQ=255)
ggTTTTAGCGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGCAACGTTTTGTCCGCAGCGAGTTTAAAt                                                                                                                                                                                                                                               >  2:89646/1‑250 (MQ=255)
      ggttaGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGcc                                                                                                                                                                                                                                         <  2:73202/246‑1 (MQ=255)
                      gCCTGACCCGCCTCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGAc                                                                                                                                                                                                                        <  2:210157/251‑1 (MQ=255)
                      gCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGAc                                                                                                                                                                                                                        >  2:136684/1‑251 (MQ=255)
                         tGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGTTCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGcc                                                                                                                                                                                                                     >  2:37471/1‑251 (MQ=255)
                          gACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGccc                                                                                                                                                                                                                    >  1:60905/1‑251 (MQ=255)
                             ccGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCt                                                                                                                                                                                                                 <  1:131724/251‑1 (MQ=255)
                                    atCACTGGGCATGCTCTCATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACAAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGTGTGAAGGTTACACCACCCATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATGTTGACGCTCGCTgggctct                                                                                                                                                                                                          <  2:8480/250‑6 (MQ=255)
                                              aTGCGCTAAGGTGTGAACAGGAGCGCCTGCAAGGGCGCACCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGca                                                                                                                                                                                                <  1:192907/251‑1 (MQ=255)
                                                gcgcTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAACCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGcat                                                                                                                                                                                               <  2:129774/250‑2 (MQ=255)
                                                         tttGAACAGGATAGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGCAAAGAGACTGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGt                                                                                                                                                                                    >  2:438994/1‑250 (MQ=255)
                                                              aCAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGGGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATg                                                                                                                                                                              <  2:141238/251‑1 (MQ=255)
                                                                aGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGAc                                                                                                                                                                            <  1:237028/251‑1 (MQ=255)
                                                                  gATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGACACCCTGCTGGCAAGCCAGATCATGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCAGACGGATGATTTTGACGCCCGCTGGCCGCTACTTGCGGCACAGTTGCAGTAATATGACGc                                                                                                                                                                          >  1:167443/1‑251 (MQ=255)
                                                                   aTCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGcc                                                                                                                                                                         <  1:136684/251‑1 (MQ=255)
                                                                     cGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGGAGGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGAGTTTGTCGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCg                                                                                                                                                                        <  1:438583/250‑1 (MQ=255)
                                                                         tgcaaggcgacattaaacAGCTGGTCGTGCGAATTGAAGCGCAGCAGCGTCAGCGACTGCGTGCTAAGAGAATGCGTGAGGGAGAAGGTTAAACCACCGATGAAAACCTGCTGGCAAGCCAGATCATGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATCCCGCCCGACGGGTGATTTTAACGCCAGCGGGCCGCTAATTGCGGCACAGTTGCAGGAATCTGACGACGGATg                                                                                                                                                                    >  2:258773/18‑250 (MQ=255)
                                                                            aaGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACtt                                                                                                                                                                >  2:342714/1‑251 (MQ=255)
                                                                                gcgcATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTc                                                                                                                                                             >  1:324606/1‑250 (MQ=255)
                                                                                  gcATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTCAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGCTGCAGTAATATGACGCCGGATGACTTTTCATc                                                                                                                                                          >  2:173224/1‑251 (MQ=255)
                                                                                      cAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGc                                                                                                                                                      >  2:343719/1‑251 (MQ=255)
                                                                                           aGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGtt                                                                                                                                                 <  1:199840/251‑1 (MQ=255)
                                                                                            gCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGttt                                                                                                                                                <  1:28531/251‑1 (MQ=255)
                                                                                               gTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCtt                                                                                                                                             >  1:104468/1‑251 (MQ=255)
                                                                                                     gCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGTGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACg                                                                                                                                       <  2:110170/251‑1 (MQ=255)
                                                                                                      cGTATTGAAGCACAGCTGCGTCAGCAATTGCGTAAAAAGAGAATGCGTGAGGGCGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGTCATTCTGTGAAGGTCTGCTGTCCCGTTTTGTCCGCAGCGAGTTTAAATACCGCACGACGGATGATTTTGACGCCGGCTGGCCGCTAATTCAGCAACAGTTACGGTATTATGACGCCGGCTGACTTTACATCCGGCGAGATTCTTTAAACGc                                                                                                                                      >  2:309218/1‑251 (MQ=255)
                                                                                                          ttGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCaa                                                                                                                                   >  1:149585/1‑250 (MQ=255)
                                                                                                            gAAGCGCAGCTGCGTCAGGTATTGCGTGAGAAGAGTATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGAGCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTAGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCAGGCGAGTTTCTTTAACCGCCAAAcc                                                                                                                                >  2:359182/1‑250 (MQ=255)
                                                                                                             aaGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAActc                                                                                                                               >  1:117114/1‑251 (MQ=255)
                                                                                                             aaGCGCAGCAGCGTCAGGTATTGCGTGAAAAGAGCATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCAGGCGCGTTTCTTTAAACGCCAAActc                                                                                                                               >  1:82002/1‑251 (MQ=255)
                                                                                                               gcgcAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCtt                                                                                                                             <  2:93076/251‑1 (MQ=255)
                                                                                                                 gcAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTcg                                                                                                                           >  1:341093/1‑251 (MQ=255)
                                                                                                                         gtcAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTa                                                                                                                                                                                                                                                  >  2:119399/3‑126 (MQ=255)
                                                                                                                         gtcAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTa                                                                                                                                                                                                                                                  <  1:119399/124‑1 (MQ=255)
                                                                                                                          taaGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGCGTTTCTTTAAACGCCAAACTCTTCGCGATAGGcc                                                                                                                  >  2:44011/3‑251 (MQ=255)
                                                                                                                                aTTGCGTGTAACGAGCATGAGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGACGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCCCAGTTGCAGTAATATGACGCCGGATGACTTTTCATCAGGCGAGTTTCTTTAACCGCCAAACTCTTCGCGATAGGCCTTAAcc                                                                                                            >  2:361218/1‑251 (MQ=255)
                                                                                                                                         aaaGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGa                                                                                                   <  1:148000/251‑1 (MQ=255)
                                                                                                                                            gagaATGCGTGAGGGTGAAGGTTACACCACCAATGAAACCCTGCAGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGcc                                                                                                                                                                                                             <  1:436601/144‑1 (MQ=255)
                                                                                                                                            gagaATGCGTGAGGGTGAAGGGTACACCACCGATGAAACACTGCTGGAAAGCCAGAGCCTGGCCTTCTGTGCAGGTATGCTGTCACGTTTTGACCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGcc                                                                                                                                                                                                             >  2:436601/1‑144 (MQ=255)
                                                                                                                                                aTGCGCGAGGGAGAAGGTTACACAACCGATGAAACCCTACTGGCAAGCCAGATCCTGGCCTTCTGTGACGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGCTTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTACTATGACGCCGGATGACTTTTCATCCGGCGGGTTTCTTTAAACGCCACACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCg                                                                                            >  2:278269/1‑251 (MQ=255)
                                                                                                                                                         ggTGAAGGTTACACCACCGATGAAACCCTGCTGGCACGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCg                                                                                   >  2:107006/1‑251 (MQ=255)
                                                                                                                                                            gAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCt                                                                                >  2:329448/1‑251 (MQ=255)
                                                                                                                                                                   gccccacgGATGAAACCCGTCTGGCAAGCCAGATTTTGGCCTGCTGTGAAGGTATGCTGAAACGTTTTGTCCGCAGCGAGTTTAAAGACCGCCCGACGGATGTTTTGGACGCCCGCTGGCCGTGAATTGCGGAACAGTTGCAGTAATATGACGCCGGATGACTTTTCAGCCGGCGAGTTTCTTTAAACGCCAATCATTTCGCGATGGGCCTTAACCGCCGCCACATGTTCCGCCATTTCCGGCTTCTCtt                                                                          <  1:243065/247‑1 (MQ=255)
                                                                                                                                                                   cctccaccGATGAAACCCTGCTGGCAAGTCATATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCACCGATTTTAAATACCGCCCGACGGATGATTTTGACGCCCGTTGGTCGCTAATTGCGGCACAGTTGCAGTAATATTACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGTGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTc                                                                         <  1:390118/248‑1 (MQ=255)
                                                                                                                                                                                     tgctgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCAt                                                                                                                                                           >  2:58841/1‑151 (MQ=255)
                                                                                                                                                                                     tgctgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCAt                                                                                                                                                           <  1:58841/151‑1 (MQ=255)
                                                                                                                                                                                                 aGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGACAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGtgat                                           >  1:217460/1‑251 (MQ=255)
                                                                                                                                                                                                  gATCCTTGCCTCCTGTGGAGTTATGGTGGCACGTTGTGTCCGCAGCGAGTGTAAATCCCGCCCGACGGATGATTTTGACCCCCGCTAGCCCCTAATTCGGGCACCAGTGCAGTAAAATAACGCCGGATGCCTTTTCACCCGGCGAGTCTCGTTAAGCTCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTGCCGTCATTTCCGGCTTCTCCACCAGCTAAGCAACCAGGTGTTTCAGGGtgatg                                          <  2:25448/251‑1 (MQ=255)
                                                                                                                                                                                                        ggtcTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATCCCGCCCGACGGATGATTTTGCCGCCCGCTGGCCGCTAATTGCGGCACAGTTGCATTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATagag                                    <  2:241204/248‑1 (MQ=255)
                                                                                                                                                                                                            ttCTGTGAAGGTATGATGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTCCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATc                                 >  2:36393/1‑250 (MQ=255)
                                                                                                                                                                                                                  gAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGc                          >  1:223301/1‑251 (MQ=255)
                                                                                                                                                                                                                        atgcgggctcGTTTTGCCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAAGTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACGCTTCTCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGt                    <  1:438994/242‑1 (MQ=255)
                                                                                                                                                                                                                        aTGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGAGGGATGATTTTGACGCCCGCTGGCCGCTGATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAACCTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGt                    <  2:256655/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tttGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCa         >  2:136077/1‑250 (MQ=255)
                                                                                                                                                                                                                                      tGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGCGCTGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACt      >  2:9660/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCCGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACtt     >  1:408637/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCt  <  1:343719/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGCCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCT  >  NC_000913/3814782‑3815266

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