Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,173 A→T 5.3% noncoding (168/184 nt) glmY ← sRNA activator of glmS mRNA, glmZ processing antagonist

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,691,1730AT 5.3% 380.7 / 7.6 113noncoding (168/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base A (52/50);  new base T (3/3);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 113 total aligned reads displayed.

AGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACA  >  NC_000913/2691028‑2691319
                                                                                                                                                  |                                                                                                                                                  
agcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAccc                                                                                                                                                 >  4:135969/1‑149 (MQ=255)
 gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCg                                                                                                                                                <  2:876379/149‑1 (MQ=255)
 gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCg                                                                                                                                                >  2:712194/1‑149 (MQ=255)
 gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCg                                                                                                                                                >  3:538795/1‑149 (MQ=255)
    gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCt                                                                                                                                             <  1:126754/149‑1 (MQ=255)
    gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCt                                                                                                                                             <  1:436085/149‑1 (MQ=255)
     aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCtt                                                                                                                                            <  3:685024/149‑1 (MQ=255)
     aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt                                                                                                                                            <  3:50657/149‑1 (MQ=255)
     aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt                                                                                                                                            <  1:624625/149‑1 (MQ=255)
      aTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCttt                                                                                                                                           <  4:108453/149‑1 (MQ=255)
      aTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCttt                                                                                                                                           <  3:825842/149‑1 (MQ=255)
       tCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTg                                                                                                                                          >  4:600174/1‑149 (MQ=255)
       tCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTg                                                                                                                                          >  1:44566/1‑149 (MQ=255)
         aGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCTTTGtt                                                                                                                                        <  1:563078/149‑1 (MQ=255)
          gCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTa                                                                                                                                       <  4:795151/149‑1 (MQ=255)
          gCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTa                                                                                                                                       >  1:302937/1‑149 (MQ=255)
            aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGTGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                     <  4:598134/149‑1 (MQ=255)
            aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGGGCTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                     >  2:820707/1‑149 (MQ=255)
            aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                     >  3:455674/1‑149 (MQ=255)
              aaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGa                                                                                                                                   >  4:845760/1‑149 (MQ=255)
               aaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGaa                                                                                                                                  <  4:313870/149‑1 (MQ=255)
                        aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                         <  2:387580/149‑1 (MQ=255)
                          tACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGc                                                                                                                       <  4:131501/149‑1 (MQ=255)
                            ccAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                     <  4:456956/149‑1 (MQ=255)
                            ccAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                     <  2:882478/149‑1 (MQ=255)
                            ccAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                     >  2:85032/1‑149 (MQ=255)
                               gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAc                                                                                                                  >  3:793728/1‑149 (MQ=255)
                                    cGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGAc                                                                                                             <  1:706389/149‑1 (MQ=255)
                                     gTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACACGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACg                                                                                                            >  2:367344/1‑149 (MQ=255)
                                      tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGt                                                                                                           >  3:560548/1‑149 (MQ=255)
                                            gCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGcc                                                                                                     <  3:206407/149‑1 (MQ=255)
                                                 gAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATg                                                                                                >  1:152169/1‑149 (MQ=255)
                                                   aaaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                              <  2:721931/149‑1 (MQ=255)
                                                    aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                             <  3:295030/149‑1 (MQ=255)
                                                    aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAATCGACGTTTCGCCCGATGGtt                                                                                             >  4:218597/1‑149 (MQ=255)
                                                          gCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata                                                                                       >  1:344764/1‑149 (MQ=255)
                                                          gCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata                                                                                       >  1:779908/1‑149 (MQ=255)
                                                          gCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata                                                                                       >  1:732221/1‑149 (MQ=255)
                                                              gCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTa                                                                                   <  3:845760/149‑1 (MQ=255)
                                                                gTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                 >  3:395536/1‑149 (MQ=255)
                                                                 tttCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGc                                                                                <  4:795838/149‑1 (MQ=255)
                                                                 tttCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAACTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGc                                                                                <  1:156356/149‑1 (MQ=255)
                                                                     aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                              >  4:480418/1‑131 (MQ=255)
                                                                     aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                              <  3:480418/131‑1 (MQ=255)
                                                                     aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGat                                                                            <  1:33363/149‑1 (MQ=255)
                                                                      aGGTGTTACTCTCGTCAGACACGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGata                                                                           <  1:367344/149‑1 (MQ=255)
                                                                       ggTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata                                                                                         >  1:431400/1‑134 (MQ=255)
                                                                       ggTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata                                                                                         <  2:431400/134‑1 (MQ=255)
                                                                       ggTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTAc                                                                                  >  3:439150/1‑141 (MQ=255)
                                                                        gtgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTg                                                                                            >  3:729322/1‑130 (MQ=255)
                                                                        gtgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTg                                                                                            <  4:729322/130‑1 (MQ=255)
                                                                         tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACCAAGCCGGC‑AATTACCCGGCTTTGTTATGGAATAATGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATGTCa                                                                        <  4:330468/149‑1 (MQ=255)
                                                                         tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                        <  1:85032/149‑1 (MQ=255)
                                                                         tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                        <  4:438033/149‑1 (MQ=255)
                                                                              ctctCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                        <  3:102017/144‑1 (MQ=255)
                                                                              ctctCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                        >  4:102017/1‑144 (MQ=255)
                                                                                 tcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa                                                                     >  4:663345/1‑144 (MQ=255)
                                                                                 tcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa                                                                     <  3:663345/144‑1 (MQ=255)
                                                                                 tcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTgg                                                                >  2:729486/1‑149 (MQ=255)
                                                                                 tcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCGACGCTGATATCAGAAGTTgg                                                                >  3:725118/1‑149 (MQ=255)
                                                                                  cGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGa                                                               <  1:820532/149‑1 (MQ=255)
                                                                                  cGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGa                                                               <  3:225469/149‑1 (MQ=255)
                                                                                   gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                              <  1:75080/149‑1 (MQ=255)
                                                                                          gcgAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggcaggc                                                       >  4:302238/1‑149 (MQ=255)
                                                                                               tAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCtt                                                  >  4:472117/1‑149 (MQ=255)
                                                                                                    tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGGCGGCAGGCACCTtgttgt                                             >  3:612487/1‑149 (MQ=255)
                                                                                                    tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                             <  2:152169/149‑1 (MQ=255)
                                                                                                      aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGTCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                           <  4:80428/149‑1 (MQ=255)
                                                                                                      aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTTGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                           <  1:186383/149‑1 (MQ=255)
                                                                                                      aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                           <  4:697306/149‑1 (MQ=255)
                                                                                                             ccGAGGGGAAGTTCAGATACAGCAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATGAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                    <  2:732221/149‑1 (MQ=255)
                                                                                                             ccGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                    >  2:738197/1‑149 (MQ=255)
                                                                                                             ccGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                    <  3:302238/149‑1 (MQ=255)
                                                                                                             ccGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                    >  4:120111/1‑149 (MQ=255)
                                                                                                                aGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTa                                 <  3:356771/149‑1 (MQ=255)
                                                                                                                aGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTa                                 <  1:718781/149‑1 (MQ=255)
                                                                                                                     aaGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTa                            <  3:478908/149‑1 (MQ=255)
                                                                                                                     aaGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTa                            >  1:331534/1‑149 (MQ=255)
                                                                                                                       gTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATg                          <  3:453729/149‑1 (MQ=255)
                                                                                                                        ttCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGt                         >  1:597267/1‑149 (MQ=255)
                                                                                                                           aGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGc                      >  2:619546/1‑149 (MQ=255)
                                                                                                                           aGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGc                      >  3:245385/1‑149 (MQ=255)
                                                                                                                            gATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCa                     <  4:560548/149‑1 (MQ=255)
                                                                                                                                 aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                >  2:384572/1‑149 (MQ=255)
                                                                                                                                    aaaGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa             >  1:161318/1‑149 (MQ=255)
                                                                                                                                    aaaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa             >  2:521705/1‑149 (MQ=255)
                                                                                                                                     aaGCCGGGTAATT‑CCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa            >  2:549724/1‑149 (MQ=255)
                                                                                                                                     aaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTGCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa            >  1:781352/1‑149 (MQ=255)
                                                                                                                                      aGCCGGGTAATT‑CCCGGCTCTGTTGTGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg           >  4:387852/1‑149 (MQ=255)
                                                                                                                                      aGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg           >  3:649001/1‑149 (MQ=255)
                                                                                                                                      aGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg           >  3:440404/1‑149 (MQ=255)
                                                                                                                                       gCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAgg          <  4:421235/149‑1 (MQ=255)
                                                                                                                                        ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAt                                >  4:375213/1‑126 (MQ=255)
                                                                                                                                        ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAt                                <  3:375213/126‑1 (MQ=255)
                                                                                                                                        ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAggg         <  3:288930/149‑1 (MQ=255)
                                                                                                                                         cGGG‑AATTACCGGGCTTTGTTATGGAATAACGCGGTGCCTAACTCGACGTTTTGCCCGATGATTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTCTATGTAGCACGTCCCGCAggcg        >  4:856197/1‑147 (MQ=255)
                                                                                                                                          ggg‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGc       >  4:477989/1‑149 (MQ=255)
                                                                                                                                              aaTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGa    >  3:336487/1‑149 (MQ=255)
                                                                                                                                              aaTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGa    >  2:283330/1‑149 (MQ=255)
                                                                                                                                                ttACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACa  <  1:738197/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACA  >  NC_000913/2691028‑2691319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑