Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,169 A→T 8.6% noncoding (172/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,1690AT 8.6% 228.8 / 11.1 70noncoding (172/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base A (27/37);  new base T (3/3);  total (30/40)
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.35e-01

CATCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAG  >  NC_000913/2690895‑2691413
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
cATCTCCAGCGCCGCGTTGGTCATGGCTTCACTGCGCCGGGCAGCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCGAAGCGCATTCAGGCTTTGACAGGCTTTCCATTCCAGCACCTACAGGGGCAGCAGAATCTGCAAAAATGCCAGACTTACCAGGTGTCGTAATGAGCGGGGAAAAAAGGGCCCGCGTTTCAAGGTGTTATCTTCGCCAGACCCGAATAGCCTGATCCCAACCGAGGGGAGGCTCGGATACAACAAAGCCGGGa                                                                                                                                                                                                                                                      >  4:82405/1‑275 (MQ=255)
       agggccgcGTTGGTCATCGCTTCCCTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGAGGAAGTTCAGATACAACAAAGCCGGGAATTAcc                                                                                                                                                                                                                                                <  2:316445/271‑1 (MQ=255)
            cgcgTTGGTCACCGCTTCACTGCGCGGGGCAACAATAAGCGTGGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCtt                                                                                                                                                                                                                                          <  4:110648/275‑1 (MQ=255)
                 tGGTCATCGCTTAACTGGGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGGGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTa                                                                                                                                                                                                                                     <  2:300417/275‑1 (MQ=255)
                 gggTCATCGCGTCACTGCGCCGGGCATCAATAAGCGTAGTTCGGTTAACCAGCGCCGCCGGATCTCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTa                                                                                                                                                                                                                                     <  1:47754/274‑1 (MQ=255)
                                     cGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTGGTTATGGAATAAGGCGGTGCCTa                                                                                                                                                                                                                  >  2:432061/1‑274 (MQ=255)
                                       ggCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACt                                                                                                                                                                                                               >  2:179193/1‑275 (MQ=255)
                                                         gCGGTTAACCCGCGCCGCCTGATGCCTAAGCGAATTCAGGTCTTGCAAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATgg                                                                                                                                                                                             <  4:179191/275‑1 (MQ=255)
                                                          cGGTTAACCGGCGCCGCACGGTCGCTAAGCGCAGTAAAGCTGTGCCAGGGTTGCACTGCCAGAACAAACAGGGGTAGCAGAATCGGCAAAAATGCCAGCATTACAAGGTGTCTGAATGAGCGGGGAAAACGGGGCCAGAGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAAGCGTGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                                                                                                                            <  4:131874/275‑1 (MQ=255)
                                                                                          aTTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTGTGTTATGGAATAAGGCGGTGCATAACTCGACGTTTCGCCCGATGGTTGAACTAGCTACTGTGATAGCACAAGTTTGaa                                                                                                                                                            >  4:229751/1‑274 (MQ=255)
                                                                                                 cTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTCACCGAGGGGAAGTTCAGATACAACCAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACCAGACGTTTCGCCCGATGGTTGATATAGCTACGTTCATATCGGAAGTTGGACggaaggc                                                                                                                                                     >  4:387599/1‑275 (MQ=255)
                                                                                                                             aCAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAACCGGACCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGCAGTTCAGATACACCAAAGCCGGGACTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACACGACGTTTCGCCCGATGGTTGATATAGCTACGCTCATATCAGAAGTTGGACGGCAGCCACCTTGTTGTGCCTCATTCGTCTTTTAt                                                                                                                         >  4:457070/1‑275 (MQ=255)
                                                                                                                               aGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGaa                                                                                                                                                                                                                                                     <  6:474955/149‑1 (MQ=255)
                                                                                                                                   gcagcaGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTAc                                                                                                                                                                                                                                                 >  5:456851/1‑149 (MQ=255)
                                                                                                                                         gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCt                                                                                                                                                                                                                                           <  8:79332/149‑1 (MQ=255)
                                                                                                                                          aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCtt                                                                                                                                                                                                                                          <  8:115929/149‑1 (MQ=255)
                                                                                                                                           aTCAGCAAAGATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAGAAGCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGGCGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACTACAAAGCCGGGTAATTCCCGGCttt                                                                                                                                                                                                                                         <  7:203230/149‑1 (MQ=255)
                                                                                                                                           aTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCttt                                                                                                                                                                                                                                         <  7:454169/149‑1 (MQ=255)
                                                                                                                                                 aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATg                                                                                                                                                                                                                                   <  6:305232/149‑1 (MQ=255)
                                                                                                                                                    aaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCCTTGTTATGGaa                                                                                                                                                                                                                                <  7:307850/149‑1 (MQ=255)
                                                                                                                                                             aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                                                                                                                       <  7:270278/149‑1 (MQ=255)
                                                                                                                                                             aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                                                                                                                       <  6:139937/149‑1 (MQ=255)
                                                                                                                                                                aCCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGcc                                                                                                                                                                                                                    <  8:196550/149‑1 (MQ=255)
                                                                                                                                                                 ccAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                                                                                                                   >  5:62008/1‑149 (MQ=255)
                                                                                                                                                                  cAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTa                                                                                                                                                                                                                  >  5:6401/1‑149 (MQ=255)
                                                                                                                                                                      tGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTc                                                                                                                                                                                                              <  8:273963/149‑1 (MQ=255)
                                                                                                                                                                         cGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGAc                                                                                                                                                                                                           <  6:456851/149‑1 (MQ=255)
                                                                                                                                                                          gTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGAAGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAGGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCa                                                                             <  2:403135/274‑1 (MQ=255)
                                                                                                                                                                           tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCa                                                                             <  2:229755/273‑1 (MQ=255)
                                                                                                                                                                              tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                                                                                                                      <  7:78735/149‑1 (MQ=255)
                                                                                                                                                                              tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                                                                                                                      >  8:78735/1‑149 (MQ=255)
                                                                                                                                                                              ggAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGt                                                                        <  2:457075/274‑1 (MQ=255)
                                                                                                                                                                                  cGGGGAACACCGGGCCAGCGTTTTAAGGTGTTACTCTCGTCACACGCGAATAGCCTGACGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTCGTTATGGATTCAGGCGGGTCCTAACTCGACGTTTCCCCCGATGCTTGACATAGCCCCGCTTATCGCAGAAGTTTGCCGGTTGGGCCACAGTTGTGCGTCCTTTGTTTTTTCTGTAGCACTTCCCGAAGGGGCTGACATAAGCCGGTCAATgt                                                                                 >  4:403130/1‑261 (MQ=255)
                                                                                                                                                                                         aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                           <  5:37022/149‑1 (MQ=255)
                                                                                                                                                                                             gggCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGata                                                                                                                                                                                       >  7:395061/1‑149 (MQ=255)
                                                                                                                                                                                                      tttCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGc                                                                                                                                                                              <  6:62008/149‑1 (MQ=255)
                                                                                                                                                                                                             gtgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATc                                                                                                                                                                       <  8:395061/149‑1 (MQ=255)
                                                                                                                                                                                                                 tACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCCTTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCAc                                                                            >  2:435736/1‑236 (MQ=255)
                                                                                                                                                                                                                        ggcAGACGCGAATAGCCTGATTCTGAACAAGGGGAAGTTCAGATATAACAAAGACGTGACTTAACCGGCTGTGTTATGGAAAAAGGCGGGTGCTAAATGTGCGTTTCGCCGGAAGGTTGATATGGCTACGATGATATCAGAAGTTGGGCGCCAGGCACCTTGGTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGTTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTgg                              <  4:427818/273‑1 (MQ=255)
                                                                                                                                                                                                                            gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                                                                                                                        <  7:344130/149‑1 (MQ=255)
                                                                                                                                                                                                                            gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                                                                                                                        <  5:52325/149‑1 (MQ=255)
                                                                                                                                                                                                                             aCGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATg                                                                                                                        >  3:60128/1‑180 (MQ=255)
                                                                                                                                                                                                                               gcgAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                                                                                                                   >  7:442727/1‑87 (MQ=255)
                                                                                                                                                                                                                               gcgAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                                                                                                                   <  8:442727/87‑1 (MQ=255)
                                                                                                                                                                                                                               gcgAAGAGCCTGATGCAAACCGCGGGGAAGGTCACATACAACAAAGCCGTGAATTACCCGCCTTTGTTATGGATTAAGTCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGCCGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTa                       <  2:226873/275‑1 (MQ=255)
                                                                                                                                                                                                                                 gAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGc                                                                                                                    <  1:60128/180‑1 (MQ=255)
                                                                                                                                                                                                                                  aaTAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGt                                                                                                                       >  4:130199/1‑176 (MQ=255)
                                                                                                                                                                                                                                      gCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCa                                                                                                                   <  2:130200/176‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                                                                                                                             <  6:148293/149‑1 (MQ=255)
                                                                                                                                                                                                                                         tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGa                                                                                                                                                             >  4:410176/1‑131 (MQ=255)
                                                                                                                                                                                                                                         tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                                                                                                                           >  7:186139/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                                                                                                                           >  7:296272/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTc                                                                                                                                      >  2:261999/1‑152 (MQ=255)
                                                                                                                                                                                                                                             gCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggc                                                                                                                                                         <  2:410181/131‑1 (MQ=255)
                                                                                                                                                                                                                                               tAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGGTGTGCGTCATTc                                                                                                                                  <  4:261995/152‑1 (MQ=255)
                                                                                                                                                                                                                                                     aGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCAc                                                                            >  2:462421/1‑200 (MQ=255)
                                                                                                                                                                                                                                                      ggggAAGGTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAg  <  2:82406/273‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAttt                                                                                                                            <  5:437164/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGc                                                                                                                    <  8:186139/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCa                                                                                                                   >  6:73561/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aTACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCAc                                                                                                                  >  5:456588/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    acaacaAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGt                                                                                                                >  5:98362/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       acaaAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCg                                                                                                                                 >  2:317439/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                         aaaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                                                                                                              <  6:422686/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                                                                                                              >  5:422686/1‑146 (MQ=255)
                                                                                                                                                                                                                                                                           aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                         >  8:454277/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                         >  7:125258/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                         >  7:102426/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAtt                                                                                                                             <  4:317435/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                         >  8:375266/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             ccGGGACTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTTCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGATGACTTAAGTCGGTGAATGAGCCACTGGTTACTCTTATCCAGTACCTGTCCCAAAAAAGAAAATAGTTTGGTAACGTACTTGTTCTACGTTGcct      >  4:489385/1‑244 (MQ=255)
                                                                                                                                                                                                                                                                                 gAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAg  <  2:489390/246‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aaTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGa                                                                                                  <  6:456588/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
CATCTCCAGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAG  >  NC_000913/2690895‑2691413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: