Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,171 T→A 8.1% noncoding (170/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,1710TA 8.1% 235.4 / 5.8 71noncoding (170/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base T (27/38);  new base A (3/3);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAG  >  NC_000913/2690902‑2691413
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
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)
                                                                                                                            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                                                                                                                                                                                                                       <  6:139937/149‑1 (MQ=255)
                                                                                                                                                      aTTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGt                                                                                                                                                                                                                       <  7:270278/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                                                                                                              >  5:422686/1‑146 (MQ=255)
                                                                                                                                                                                                                                                                  aaaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                                                                                                              <  6:422686/146‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                         >  7:102426/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                         >  7:125258/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                         >  8:454277/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)
                                                                                                                                                                                                                                                                             ttACCCGGCTTTGTTATGGAATAAGGAGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATAACAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACa                                                                                                >  8:408154/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
AGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAG  >  NC_000913/2690902‑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: