Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,454 2 bp→AA coding (97‑98/993 nt) 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,4540TA100.0% 253.2 / NA 69Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (29/40);  total (29/40)
*NC_0009132,867,4550GA98.6% 243.8 / NA 69Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/40);  minor base C (1/0);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

AGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACT  >  NC_000913/2867211‑2867691
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
aGACCACGATTGCCATACCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                                                                                                                                                                                                         <  2:148301/251‑1 (MQ=255)
aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                                                                                                                                                                                                         <  1:354529/251‑1 (MQ=255)
        aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCt                                                                                                                                                                                                                                 <  1:55661/251‑1 (MQ=255)
            ccATAACGGCGGGCAATTTTTACCGCCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTc                                                                                                                                                                                                                             <  2:324475/251‑1 (MQ=255)
                cccggcggGCAATTTTTACCACCAGACGCCGGTTTCTCGCGGTCATCCTGCGGCGAGAGGCGTCATCTCCACGCAGTGCGCGACGCGCAAAATAGACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACTCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTGTCACTGGGGTCCAATTCTACTAAGGCTTTTTCGTCa                                                                                                                                                                                                                         <  2:121061/249‑1 (MQ=255)
                      gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAcc                                                                                                                                                                                                                   <  1:499256/251‑1 (MQ=255)
                      gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAAcc                                                                                                                                                                                                                   >  2:37946/1‑251 (MQ=255)
                        gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                                                                                                                                                                 <  1:167023/251‑1 (MQ=255)
                                 accaccAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                                                                                                                                                        <  1:364878/251‑1 (MQ=255)
                                           gCAAGTTCCTCTCGATCAGCCGGCGGCGAGAGGCGACATCTCCAAGCAGTGCGCGACGCGGAAAATAAACTTCTTCGTCGTCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCTGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCATAAACCTCAACTCCGTTCTCATCAAAt                                                                                                                                                                                              <  2:335230/251‑1 (MQ=255)
                                           gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAt                                                                                                                                                                                              <  1:186506/251‑1 (MQ=255)
                                              aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTcc                                                                                                                                                                                           >  2:108990/1‑251 (MQ=255)
                                                      tcgctcatcCGGCGGCGAGAGGCAACATCTCCCCGCAGTACGCGACCCACAAAATAAACTTCTTCTTCGGACGTTAACAGAGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCCACACACGCAGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCAGTGGGTTCCAATTATACTACGGCCTTTTCGACAAAAAACTCAACTCAGTTCTCATCAAATTCCGCATCTTc                                                                                                                                                                                   >  2:483762/5‑251 (MQ=255)
                                                         atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCATCTTCAtt                                                                                                                                                                                <  1:48531/251‑1 (MQ=255)
                                                                ggcggcGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                          >  1:121810/1‑250 (MQ=255)
                                                                   ggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                       >  1:448922/1‑250 (MQ=255)
                                                                          ggCGACATCGCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATGCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTGTCa                                                                                                                                                               <  2:230988/251‑1 (MQ=255)
                                                                            cGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGc                                                                                                                                                             >  2:377546/1‑251 (MQ=255)
                                                                                atcaccACGCAGTGCGCGACGCGCAAAAGAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAt                                                                                                                                                         <  2:65262/247‑1 (MQ=255)
                                                                                   tcCACGCAGTGCGCGACGCGCAAAATAACCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGATGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCCCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAAATACGTTCTCCTCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCt                                                                                                                                                      >  2:282642/1‑251 (MQ=255)
                                                                                    ccACGCAGTGCGCGACGCGCAAAATAAACTTCTCCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCt                                                                                                                                                      <  1:484866/250‑1 (MQ=255)
                                                                                     cACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                     >  1:54671/1‑250 (MQ=255)
                                                                                          aGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAt                                                                                                                                               >  1:495241/1‑251 (MQ=255)
                                                                                                      gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑Ct                                                                                                                                  <  1:480936/251‑1 (MQ=255)
                                                                                                        aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTAc                                                                                                                                >  2:33585/1‑251 (MQ=255)
                                                                                                         aaaTAAACTGCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTAcc                                                                                                                               <  2:261906/251‑1 (MQ=255)
                                                                                                             aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTTAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTg                                                                                                                           <  1:363207/251‑1 (MQ=255)
                                                                                                                    ttcttcGGCCGTTAACAGCGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGGATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCtt                                                                                                                    >  1:73694/1‑251 (MQ=255)
                                                                                                                       ttcGGCCGTTAACAGTGGTGCATAACCAATCCCACCAAGGAAATGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAACTCATGAACTTTCAGCGTATTCTGACTCATAAGGTAGCTC‑CTACCCGCGATCCCTTGAc                                                                                                                 >  2:420400/1‑251 (MQ=255)
                                                                                                                        tcGGCCGTTAACAGTGGTGAATAACCCATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGGGGCTC‑CTACTCGTGATCCCTTGAc                                                                                                                 <  2:423060/250‑1 (MQ=255)
                                                                                                                          ggCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACgg                                                                                                               >  1:70471/1‑250 (MQ=255)
                                                                                                                           gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCCGCGTATTCTGACTCATAAGGTGGCTA‑CTACCCGTGATCCCTTGACGGa                                                                                                              >  2:225355/1‑250 (MQ=255)
                                                                                                                             cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAAc                                                                                                            >  1:425113/1‑250 (MQ=255)
                                                                                                                               ttAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACAtt                                                                                                         >  2:94881/1‑251 (MQ=255)
                                                                                                                                 aaCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTc                                                                                                        <  1:364209/250‑1 (MQ=255)
                                                                                                                                  aCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCAAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCa                                                                                                       >  1:152043/1‑250 (MQ=255)
                                                                                                                                       ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCaaa                                                                                                 <  1:413776/251‑1 (MQ=255)
                                                                                                                                             tAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                                                                           >  1:171243/1‑251 (MQ=255)
                                                                                                                                              aaCCAATCTCCCCAAGGTAAAGCTGCGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCCCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCCTCAACTTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCACGCCACAGCCTgg                                                                                          >  2:157217/1‑251 (MQ=255)
                                                                                                                                              aaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTgg                                                                                          >  1:66770/1‑251 (MQ=255)
                                                                                                                                                   aTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccg                                                                                     >  1:240122/1‑251 (MQ=255)
                                                                                                                                                       cACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCACTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgccg                                                                                  >  1:96885/1‑250 (MQ=255)
                                                                                                                                                          cAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                                                                              <  1:150445/251‑1 (MQ=255)
                                                                                                                                                             ggTAAAGCTGTGCCGCGTCCAACACACGCGGT‑GGGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTGCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCGTGACAGAACATTCAGGCAAAAGCCTGGTGCCGCCGATTTATc                                                                           <  2:426804/251‑1 (MQ=255)
                                                                                                                                                              gtaatgcggAGTCTCGTCCAACACACG‑TGTTGTGGCTCCCTTCGATAACAGGTCCTCGTCGGCCAAATCTTTATCCCTGGGTTCCAATTCTACGAAGGCGTTTTCTTCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCATCTTCATTTAAATCATGAACCTTCAGCGTTTACTTACTCATAAGGTGG‑TCTCTGTCCTTGCTCCCTTGACGGAACATTCAATCTAAATCCTTGTTCCGCCCATGTATcc                                                                          <  2:333249/243‑2 (MQ=255)
                                                                                                                                                                aaaGCTGAGTCGCGTCCAACACACGCTGT‑GGGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCACTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCt                                                                        >  1:154666/1‑251 (MQ=255)
                                                                                                                                                                   gCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCg                                                                     >  2:354035/1‑251 (MQ=255)
                                                                                                                                                                       aGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCaa                                                                 <  1:37946/251‑1 (MQ=255)
                                                                                                                                                                           gcggccAACACCCGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTGCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGCATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGTTTCCGCCGATTTATCGCTGCGGCAAATaa                                                             <  2:121810/247‑1 (MQ=255)
                                                                                                                                                                                 aaCACTCGCTGT‑GGTGCTGCCTGTGATAACAGTACCTCTTCTGCCAAATCTTTATCACTGGGTTCCAATTCTACTACGGCACTTTCTTCAAAAACCTCAACCCCGTTCTCATCAAATTCTGCAGCTTCATTTAAATCAGGAACTTTGAGCGTATTCTAACTCATAAGGTGGTTC‑CTACCCTTGATCCCTTAACGGAACATTCACGCAAAAGCCTGGTTCCTCCGATTTATCGCTGCGGCAAATAACGCAGc                                                       <  2:472437/251‑1 (MQ=255)
                                                                                                                                                                                       cGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGttt                                                 <  1:377546/251‑1 (MQ=255)
                                                                                                                                                                                        cgttgtctgtcTCCCTGCGATAACAGTTACTCTTCGGCCACATCGTTATCACTTGTTTCCAATTCTTGTAAGGCCTTTTAGTTAAAAACCTCACCTCCGTTCTCATCAAATTCTGCATCTACATTTAAATCATTAACTTTCCGCGTATTCTTACTCATAAGGTGGCTC‑CTACCCGTTATCCCCTGACGGAACATGCAAGCAAAAGCTTGTCTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCTGGttt                                                 <  2:104719/241‑1 (MQ=255)
                                                                                                                                                                                           gt‑gtgGCTCCCTGCTATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTTCTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTTGCTC‑CTACCCGTGATCCCTGGACTGAAGATTTAAGCAAAAGCCGGGTTCCGCCGATTTATCGCTGCGGCAACTAACGCATCGGGTTTAc                                               <  2:383266/249‑1 (MQ=255)
                                                                                                                                                                                              gtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAt                                           <  2:171243/251‑1 (MQ=255)
                                                                                                                                                                                                 gCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTc                                        <  1:325163/251‑1 (MQ=255)
                                                                                                                                                                                                   tCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTccc                                      >  2:438712/1‑251 (MQ=255)
                                                                                                                                                                                                    cccTGTGATAACAGTTCCTCTTCGGCAAAATCGTTATCACTGGGTTCCAATTCTACTAAGTCCTTTTCGTCAAAAACCTCAACCCCGTTCTCATCAAATTCCGCAGCTTCATTGAACTCCTGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAAAATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTcccc                                     <  2:54671/251‑1 (MQ=255)
                                                                                                                                                                                                    cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTcccc                                     >  2:43174/1‑251 (MQ=255)
                                                                                                                                                                                                    cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTccc                                      <  1:446177/250‑1 (MQ=255)
                                                                                                                                                                                                        gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGt                                 <  1:215063/251‑1 (MQ=255)
                                                                                                                                                                                                                attgccTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAAGCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGTGC‑CTTCCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTGATCTCTGCGGCAGAGAGCGCAGCGGGTAGACGGATGGCCCCTTGTAACGAAtt                         <  2:152043/247‑1 (MQ=255)
                                                                                                                                                                                                                aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCt                                                                        <  1:140861/204‑1 (MQ=255)
                                                                                                                                                                                                                aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCt                                                                        >  2:140861/1‑204 (MQ=255)
                                                                                                                                                                                                                aGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAtt                         >  1:206049/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                      <  2:23114/173‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                      >  1:23114/1‑173 (MQ=255)
                                                                                                                                                                                                                                   cGTTATCACTGGGTTCCAATGCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCAGCAAAGTCCGCATCTTCATTTAAATCATGAACGTGCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAGTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGtt       <  2:70471/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACt  <  1:274604/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACt  <  1:94881/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
AGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGT‑GTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTC‑CTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACT  >  NC_000913/2867211‑2867691

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