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,4550GA97.2% 115.0 / NA 37Q33* (CAG→TAG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/19);  minor base C (1/0);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCG  >  NC_000913/2867219‑2867698
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
aTTGCCAGAACGGCGGGCAATTTTTTCCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCGTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGGGTCCAACACACGTTTTGTGGCTCCCTGCGATAACAGTTCCGCGTCGGCCAAATCGTTATCACTGGGGTCCTATTCTACTAAGGCCt                                                                                                                                                                                                                                       <  2:150256/251‑1 (MQ=255)
    ccATAACGGCGGGCCATTTTTACACCCAGACGCAAATTACTCTCCATCCTCCGGCGGCGAGGGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGGGGTGAATAAACAATCTCACCAAAGTAAAACTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGGTCCTCTTTGGCCCAATCGTTATCACTGGGTTCCTTTTTTACTCAAGCCTTTTc                                                                                                                                                                                                                                   >  2:210976/1‑251 (MQ=255)
              gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGCGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAAcc                                                                                                                                                                                                                         >  1:172121/1‑251 (MQ=255)
                  aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCa                                                                                                                                                                                                                      >  1:122126/1‑250 (MQ=255)
                     ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTc                                                                                                                                                                                                                  >  1:101118/1‑251 (MQ=255)
                        tACCACCAGACGCAAGTTACTCTCGATCATCCGGTGGCGAGCGGCGCCATCTCCACGAAGTGCACGACGCGCACCATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACAAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACCGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGt                                                                                                                                                                                                               >  2:148362/1‑251 (MQ=255)
                           caccaGACGCAAGTTACGCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCTTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtct                                                                                                                                                                                                            <  1:155308/251‑1 (MQ=255)
                                     aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTc                                                                                                                                                                                                  <  2:33086/251‑1 (MQ=255)
                                      aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTcc                                                                                                                                                                                                 <  1:271129/251‑1 (MQ=255)
                                        ttACTCTCGCTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGATCTCATCAAATTCCGc                                                                                                                                                                                               >  1:281600/1‑251 (MQ=255)
                                        tgaCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCTCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCTATTCTACTAAGGCCGTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg                                                                                                                                                                                                <  2:102227/248‑1 (MQ=255)
                                                          cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACAAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGACAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                             >  1:158136/1‑251 (MQ=255)
                                                                 agGCGACATCTACACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTc                                                                                                                                                                      >  1:50352/1‑251 (MQ=255)
                                                                       cATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                                                                                                                 >  2:263016/1‑250 (MQ=255)
                                                                         cccccACGCAGTGCGTGCCTCGTAAAATAAACGTCTGCTTCGTCCGTTAACAGTTTCGAATAACCAATCTCTCCAAGGTAAAGCGGAGTCGCGTCCAACACACGCTGTGTTGCTCCCTGCGATAACAGTTTCTCTGCGGCCAACTCGTCATCACTGGGGTTCTAATCTACTAAGGCTTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAATTCATGAACCTTCAGCGTAtt                                                                                                                                                              <  2:11479/248‑1 (MQ=255)
                                                                                  aGTGCGCGACGCGCAAAATACACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCACAACCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAt                                                                                                                                                     >  2:17840/1‑251 (MQ=255)
                                                                                         gACGCGCAAAATAAACTTCTTCTTCGGCCGGTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTgg                                                                                                                                              <  2:301086/251‑1 (MQ=255)
                                                                                               cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa                                                                                                                                                                                    <  2:86749/207‑1 (MQ=255)
                                                                                               cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTa                                                                                                                                                                                    >  1:86749/1‑207 (MQ=255)
                                                                                                      aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGa                                                                                                                                 <  1:337811/251‑1 (MQ=255)
                                                                                                       aCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACATTTCCGCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTTCTAAGGCCTTTTCGTCAAAATGCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGAt                                                                                                                                <  1:246755/251‑1 (MQ=255)
                                                                                                                tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGAc                                                                                                                        <  2:292924/250‑1 (MQ=255)
                                                                                                                    ccGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                   <  1:79595/251‑1 (MQ=255)
                                                                                                                                                                   gcgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATaa                                                                    <  2:106936/251‑1 (MQ=255)
                                                                                                                                                                                   gtgtgGCTCCCTGCGATAACATTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACgg                                                    <  1:189785/251‑1 (MQ=255)
                                                                                                                                                                                       ggCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAtt                                                 <  1:263016/250‑1 (MQ=255)
                                                                                                                                                                                           cccTGAGATAAAAGTTCATCTTAGCCCAAATGCTTATCACTGGGTTCCTATACTACTAACGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAGTTCCGAATCTTCATTTAAATCCTAAACTTTCAGCGCACTCTGACTCATAAGGTGGCTCCTACACGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCAGCGGCAACTAACGACGCGGGTTTACGGATTTccac                                            >  2:181975/1‑249 (MQ=255)
                                                                                                                                                                                            ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATaa                                                                                                                                                   >  2:263853/1‑147 (MQ=255)
                                                                                                                                                                                            ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCAGCTTCATTTAAAGCATGAACGTTCAGCGTATTCTGACTCATaa                                                                                                                                                   <  1:263853/147‑1 (MQ=255)
                                                                                                                                                                                                      cAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAt                                 <  2:101118/251‑1 (MQ=255)
                                                                                                                                                                                                                       aaTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtg                >  2:105487/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtgt               >  1:65997/1‑251 (MQ=255)
                                                                                                                                                                                                                                cACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAACAGCCTGGTTCCGCCGATTTATCGCTGCGGCAACTAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTgg       >  2:142872/1‑251 (MQ=255)
                                                                                                                                                                                                                                cACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTg        <  1:143653/250‑1 (MQ=255)
                                                                                                                                                                                                                                  cTGGGTTCCTCTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTgttt     >  1:238175/1‑248 (MQ=255)
                                                                                                                                                                                                                                  cTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCCTTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGCACTGGtt     >  1:223815/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggTTCCTATTCTACTAAGGCCTGTTCGTCAAAAACCTCAACTCCGTTCTTATCAAATTCCGCATCATCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTCCCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTAGCGCTGCGTCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGTTTCCg  <  2:238175/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
ATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCG  >  NC_000913/2867219‑2867698

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: