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,4540TA97.8% 150.2 / NA 45Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (21/23);  minor base C (0/1);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,867,4550GA100.0% 153.8 / NA 45Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/24);  total (21/24)

ACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867213‑2867697
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  
aCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctacta                                                                                                                                                                                                                                            <  1:192161/251‑1 (MQ=255)
                      gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                                                                                                                                                                      <  1:189171/251‑1 (MQ=255)
                              tACCACCAGACGCAAGTTACTCTCGATCATCCGGCTGCGAGAGGCAACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                                                                                                                                                                                                                <  2:25861/249‑1 (MQ=255)
                              tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                                                                                                                                                                                                                >  1:25861/1‑249 (MQ=255)
                                    caGACGCAAGGTCCTCTCGATCACCCGGCGGCGAAAGTGTACATCTCCCAGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTGAACAGTGGTGAATACCCAATCTCACCAAGGTAAAGCTGAGTCGCTTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTTGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat                                                                                                                                                                                                        <  2:49580/251‑1 (MQ=255)
                                                 ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATc                                                                                                                                                                                           >  1:160320/1‑251 (MQ=255)
                                                                cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                            >  1:397897/1‑251 (MQ=255)
                                                                 ggcgAGAGGCGGCATCTCCACGCAGTGCGCGACGCGCAAAATCACCTTCTTCTTCGGCCGTTAACAGTGGAGAATACCCAATCTCACCAAGGTAAAGCTGCGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCACAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                            >  2:52934/1‑250 (MQ=255)
                                                                      gagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAACTCACCAAGGTAAAGCGAAGGCGCGTCCAACACTCGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTGCGGCCTAAGCGGAATCACTGGGTGCCAATTCAACTAAGGCCTTTTCGTCAAATACCTCAACGCCGTTCCCATCAAATTCCGCATCTTCATTTAAATCATGAACttt                                                                                                                                                                      <  1:498479/251‑1 (MQ=255)
                                                                            aCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                                                                                                                <  1:266826/251‑1 (MQ=255)
                                                                                             gcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGt                                                                                                                                               <  2:235297/251‑1 (MQ=255)
                                                                                              cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg                                                                                                                                              <  2:172416/251‑1 (MQ=255)
                                                                                                   cgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCCCCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCCTAAGGTGGCTcc                                                                                                                                         >  2:341080/1‑251 (MQ=255)
                                                                                                   cgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTcc                                                                                                                                         <  1:97426/251‑1 (MQ=255)
                                                                                                     cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATCACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGCTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                                                                                                                        >  2:373954/1‑250 (MQ=255)
                                                                                                                  ttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTACAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTATTAAAATTCCGCATCTTCATTTAAATCATGCACTTTTAGCGTATTCTGACTCATTAGGTGGCTCCTACCCGTGATCCCtt                                                                                                                          >  2:184974/1‑251 (MQ=255)
                                                                                                                      tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGATGAGTCGCGTCCAACACATGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTCCTAAGGCCTTTTCGTCAAAATCCTCAACGCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg                                                                                                                                              <  1:313422/227‑1 (MQ=255)
                                                                                                                      tcGGCCGGTAACAGGGGTGAATAACCAATCTCACCACGGAAAAGCTGAGTCGCGTCCAACCCACGCTGTGGGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCAATTCTACTAAGACCTTTTCGTCAAGCACCTCAAATCCGTTCTCATCAAATTCCTCATCTTCATTTAAATCATGAACTTTCCGCGTATTCTGACTCATAAGGTg                                                                                                                                              >  2:313422/1‑227 (MQ=255)
                                                                                                                       cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                            >  1:393782/1‑180 (MQ=255)
                                                                                                                       cGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAt                                                                                                                                                                                            <  2:393782/180‑1 (MQ=255)
                                                                                                                        ggCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                                 >  2:163012/1‑238 (MQ=255)
                                                                                                                        ggCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                                 <  1:163012/238‑1 (MQ=255)
                                                                                                                                aCAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCa                                                                                                             >  1:227548/1‑250 (MQ=255)
                                                                                                                                      gtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCGTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCaaaa                                                                                                      <  1:55019/251‑1 (MQ=255)
                                                                                                                                          aTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTCTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAATAGCCt                                                                                                  <  2:188870/251‑1 (MQ=255)
                                                                                                                                          aTAACCAACCTCACCACGGCAAAGCTGAGTCGCGTCCACCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTATGACTCATAAGGTGGCTCCTACCCGTGATCCCTGGACGGAACATTCAGGCACAAGCCt                                                                                                  >  2:355208/1‑251 (MQ=255)
                                                                                                                                             aCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGt                                                                                               >  1:24923/1‑251 (MQ=255)
                                                                                                                                               cAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATCCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTc                                                                                             <  1:360107/251‑1 (MQ=255)
                                                                                                                                                          aGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTAt                                                                                  >  1:450125/1‑251 (MQ=255)
                                                                                                                                                            gTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTaaa                                                                                                                                                                                 <  2:192102/154‑1 (MQ=255)
                                                                                                                                                            gTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTaaa                                                                                                                                                                                 >  1:192102/1‑154 (MQ=255)
                                                                                                                                                                   tGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGc                                                                         >  1:407957/1‑251 (MQ=255)
                                                                                                                                                                         gcgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTTGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATaa                                                                   <  1:172695/251‑1 (MQ=255)
                                                                                                                                                                                acacacGCTGTGTGTCTCCCTGCGATAACAGTTCCTCTTCGGCCAACTCGTTATCGCTGGGATCCAATTAAACCAACGCCTTTTGGACCAAAACCTCAACTCCGTGCTCGACAAATTCCGCATATTCATTTAAATCATGAACTTTCAACGTATTCGGACTCCTAAGCAGGCTCCTACCCGTGATCCCTTGACGCAACATTCACGCAAAAACCTGCTGCCGCCGATTTATCGCTCGGGCCAATAACGccccg                                                            >  2:60337/1‑247 (MQ=255)
                                                                                                                                                                                acacacGCTCTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCTTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAGCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGTTTCCGCCGATTGATCGCTGCGGCAAATAACGCAGCg                                                            <  2:42036/251‑1 (MQ=255)
                                                                                                                                                                                     cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGttt                                                       <  1:369536/251‑1 (MQ=255)
                                                                                                                                                                                             gtcTCCCTGCGATAACAGTTCCCCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGGCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCATCTTCATTTAAATCATGAACTTGCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGTTTCCACCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAttt                                               <  1:520902/249‑1 (MQ=255)
                                                                                                                                                                                             ggCTCCCTGCGATAACAGTTACTCTAAGGACAAATCGTTATCACTGGGTTCCAATTCTACTAAGGACTTTTCGTCAAAAACCTCAACTCCGTTATCATCAACTTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGAGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGTTTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAttt                                               >  1:175687/1‑251 (MQ=255)
                                                                                                                                                                                                 cccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTcccc                                           >  2:69717/1‑251 (MQ=255)
                                                                                                                                                                                                     gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGt                                       >  2:250164/1‑251 (MQ=255)
                                                                                                                                                                                                       gATAACAGTGCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTAGTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAAc                                    <  1:472366/251‑1 (MQ=255)
                                                                                                                                                                                                        aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTGTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCGGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAAc                                    <  1:373954/251‑1 (MQ=255)
                                                                                                                                                                                                        aTAACAGTTCCTCTTCGGCCAAATCGTTATCACGGGGTTCCAATTCTAATAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGTGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTGGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCggggggacggagggaaaagtgtaac                                    <  1:495305/251‑23 (MQ=255)
                                                                                                                                                                                                           aCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGa                                  >  1:380414/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gggTTCCCATTCTACTAAGCCCTTTTCGTCAAAAACCTCAACTCCGTTCGCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCCCCTACCCGTGATCCCTTGACGGAACAGTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCATCGGGTTTACGCATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTTAACTGGTTcc  <  2:175687/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                  
ACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867213‑2867697

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