Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,455 G→A Q33* (CAG→TAG)  rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,4550GA98.2% 177.7 / NA 55Q33* (CAG→TAG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/29);  minor base C (1/0);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGT  >  NC_000913/2867221‑2867694
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
tGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCtt                                                                                                                                                                                                                                  <  1:352985/250‑1 (MQ=255)
    aaaactgcggGCAATTTTTTCCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCTCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAGCACTGGGTTCCTATTCTACTAAGGCCTTGTCGt                                                                                                                                                                                                                             <  2:116333/245‑1 (MQ=255)
         ggcggGCAATTTTTACCACCAGACGAAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGAGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATATCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGAGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGTGTTCCTATTCTACTAAGGCATTTTCGTCaaa                                                                                                                                                                                                                         >  2:414861/1‑250 (MQ=255)
          gcggGCAATTTTTACCACCAGACGCAAGTCACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCaaaaa                                                                                                                                                                                                                       <  2:399221/251‑1 (MQ=255)
            gggCAATCTTTACCACCAAACGCAAGTTACTCACGCTCACCCGCCGGCGACAGGCGACATCTCCACGCAGTGCGCGACGCGCAAACCACACTTCTTCTTCGACCGTTAACAGTGGCGAATAACCAATCTCACCAAGGAAAAGCTGCGCCGCGTACAACACACGCTGTGAGGCTCCCGCCGATAACAGCTCACCTCCAGCCAAATCGTTCTCCCTGGGCTCCTATTCTACTAACGCCTATTCCTCaccaccc                                                                                                                                                                                                                     >  2:513545/1‑245 (MQ=255)
             ggCAATTTTTCCCACCAGACGCCAGGTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCAAGCATTGCGCGACGCGCAAAATAAACTTCTGCTTCGGCCGTTAACAGTGGTTAATAACCAGTCTCACCAAGTTAAAGCTGAGGCGCGGCCAACACACGCTTTGTGGCGCCCGGCGATAACAGTTCCTCTTCGGCCATATCCTTAGCACTGGGTTCCTTTTCTATGGAGTCCTTTTCGTCATAAACCt                                                                                                                                                                                                                    <  2:223360/251‑1 (MQ=255)
                 aTTTGTACCCCCAGACGCAAGTTACTCTCGATCATCCTGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAAc                                                                                                                                                                                                                <  2:128930/251‑1 (MQ=255)
                   ttttACCACCAGACGCACGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTc                                                                                                                                                                                                              >  1:191645/1‑251 (MQ=255)
                     ttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                                                                                                                                                                                                             <  1:139441/250‑1 (MQ=255)
                      tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCGTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCg                                                                                                                                                                                                            <  1:90024/250‑1 (MQ=255)
                       accaccAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGt                                                                                                                                                                                                           <  1:18415/250‑1 (MQ=255)
                             aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCTCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                    <  2:298004/251‑1 (MQ=255)
                              gACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                    >  1:113033/1‑250 (MQ=255)
                                  cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCGACTCCGTTCTCATCAAAtt                                                                                                                                                                                               <  1:408378/251‑1 (MQ=255)
                                        cctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                         <  2:175476/250‑1 (MQ=255)
                                             cgagcatcCGGCCGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAAGAAACGTCTTCTGCGGACGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACATTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCACCTCCGTTCTCATCATATTCCGCATCt                                                                                                                                                                                       <  2:110036/244‑1 (MQ=255)
                                                    ccGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAAt                                                                                                                                                                             <  2:191645/251‑1 (MQ=255)
                                                      ggcggcGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                            >  1:50336/1‑250 (MQ=255)
                                                      ggcggagcgagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATCAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGATGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGCTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                            >  2:529853/9‑250 (MQ=255)
                                                       gcggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCAt                                                                                                                                                                          <  1:300872/251‑1 (MQ=255)
                                                        cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTGCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCGTCGGCCAAAGCGTTAGCACTGGGTTCCTATTCTACTAAGGCCTTTCCGTCATAAACCTCAACTCCGTTCTCATCAAATTCCGCATCATCAATTAAATCATg                                                                                                                                                                         <  2:370225/251‑1 (MQ=255)
                                                                         cccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCt                                                                                                                                                        <  1:137717/250‑1 (MQ=255)
                                                                          ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                       <  1:318250/251‑1 (MQ=255)
                                                                              gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTc                                                                                                                                                   <  1:534877/251‑1 (MQ=255)
                                                                               cAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTc                                                                                                                                                   >  2:376334/1‑250 (MQ=255)
                                                                                    cgcgACGCGCAAAATAACCTTCTTCTTCGGCCGTTAACAGGGGTGGATAACCAAACTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTATGGCTCCTTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTACCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTc                                                                                                                                                         <  1:551726/239‑1 (MQ=255)
                                                                                    cgcgACGCGCAAAATAAACTTCTTCTTCAGCCGTTAACAGTGGTGAATAACCAATCTCACCCAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAtta                                                                                                                                                         >  2:551726/1‑238 (MQ=255)
                                                                                            gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGATGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                                                                                                                     <  1:155527/251‑1 (MQ=255)
                                                                                                 aTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCg                                                                                                                                <  1:194420/251‑1 (MQ=255)
                                                                                                    aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                              >  1:103293/1‑250 (MQ=255)
                                                                                                         cttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCt                                                                                                                        >  1:251783/1‑251 (MQ=255)
                                                                                                               cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                                  <  1:202552/251‑1 (MQ=255)
                                                                                                                ggCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGCCgg                                                                                                                  >  2:305951/1‑250 (MQ=255)
                                                                                                                ggCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGGAAAGCTGAGTCGCGTCCAACACACGCTGGGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGACTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATCAGGTGGCTCCTACCCGTGATCCCTTGACgg                                                                                                                  >  2:419477/1‑250 (MQ=255)
                                                                                                                 gCCGTTAACAGTTGTTAATAACCAATCTCGCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGGTCCTCTTCGTCCAAATCGTTATCACTGGGTTCCTAGTCTACTAAGTCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCATCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGaa                                                                                                                <  1:244832/251‑1 (MQ=255)
                                                                                                                 gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGa                                                                                                                 >  2:413499/1‑250 (MQ=255)
                                                                                                                           gtggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGCTCCCTTGACGGAACATTCAAGc                                                                                                       >  1:217716/1‑250 (MQ=255)
                                                                                                                               tgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCaaa                                                                                                    >  2:4885/1‑249 (MQ=255)
                                                                                                                                      ccAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCCATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAACAGCCTGGtt                                                                                           >  2:425906/1‑251 (MQ=255)
                                                                                                                                                 aaGGGAAAGCTGAGTAGCGTCAAACACACGCTGTGTGGCTCCCTGCGGTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTCCTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                                                                                 >  2:88545/1‑250 (MQ=255)
                                                                                                                                                    gTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCg                                                                             >  1:251479/1‑251 (MQ=255)
                                                                                                                                                      aaaGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCt                                                                           >  2:310033/1‑251 (MQ=255)
                                                                                                                                                         gCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTGTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCg                                                                        <  2:144739/251‑1 (MQ=255)
                                                                                                                                                               gcgcgTCCAACACACGCTGTGTGGCTCCCGGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAGTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAAt                                                                  <  2:217716/250‑1 (MQ=255)
                                                                                                                                                                    tCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGc                                                             >  2:139706/1‑251 (MQ=255)
                                                                                                                                                                       aacacgagCTGTGTGGCCCCCTGCGGTAACAGTTCCTATACGGCCAAATCGTTATCACTGGGTTCCTATTATACTCAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGCCTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCACGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGc                                                          >  2:241709/8‑251 (MQ=255)
                                                                                                                                                                        acacacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCaa                                                                                                                                                                                                  >  1:497455/1‑114 (MQ=255)
                                                                                                                                                                        acacacGCTGTGTGGCTCCCTGCGATAACAGGTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTGCGTCAAAAACCTCAACGCCGTTCTCATCaa                                                                                                                                                                                                  <  2:497455/114‑1 (MQ=255)
                                                                                                                                                                          acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCAt                                                                                                                                                                          <  2:441368/136‑1 (MQ=255)
                                                                                                                                                                          acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCAt                                                                                                                                                                          >  1:441368/1‑136 (MQ=255)
                                                                                                                                                                            ccGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGtt                                                     <  1:302716/250‑1 (MQ=255)
                                                                                                                                                                                 gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTAc                                                  <  1:256501/249‑1 (MQ=255)
                                                                                                                                                                                                aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAAc                                 <  1:184010/251‑1 (MQ=255)
                                                                                                                                                                                                                           tATCACCGGGTTCCTCTTCTACTAAGCCCTTTTCGTCAACAACCTCCACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAAc      >  2:79396/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAACACTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGCTTACGGATTTCCCCTTGTAACGACTTTCAAAATGCAAGCGTGTTGAACTGGt  >  2:102649/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
TGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGT  >  NC_000913/2867221‑2867694

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: ‑