Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,354 G→C 6.8% A152G (GCC→GGC)  rsd ← stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,3540GC 6.8% 418.5 / 15.7 130A152G (GCC→GGC) rsdstationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base G (64/57);  new base C (4/5);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 136 total aligned reads displayed.

GGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATC  >  NC_000913/4196207‑4196503
                                                                                                                                                   |                                                                                                                                                     
ggTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcg                                                                                                                                                     >  2:111930/1‑150 (MQ=255)
   aTTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGgcggcg                                                                                                                                                  >  2:570229/1‑150 (MQ=255)
         cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGc                                                                                                                                            <  2:517479/150‑1 (MQ=255)
         cTGATCAGGGTCCGCGGATCCCGATTAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCTCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCCGCCCGCGTCaaac                                                                                                                                            <  2:1329409/150‑3 (MQ=255)
          tGATCAGGTTCCGCGGATACAGATTAAACGGTCTCAGCCAGTTAAGGCTCTCCGACAAGAAATTAGGCCCGCAAATGGGGCAGTGAATGTAAATAACGCGGTAACAGCGCTGAACTAAAGCACGATGTATCACGCGGGCGGCGTCAAGcc                                                                                                                                           >  2:411341/1‑149 (MQ=255)
             tCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGCTTGACGCCGCCCGCGTCAAACatcc                                                                                                                                        <  2:828745/150‑4 (MQ=255)
              cAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag                                                                                                                                       >  1:992969/1‑150 (MQ=255)
               aGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCATTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGCTTGACGCCGCCCGCGTCAAACatcctg                                                                                                                                      <  2:239117/150‑6 (MQ=255)
               aGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGAAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACTCGTGCGGCGTCAAGCACcagc                                                                                                                                      >  1:243179/1‑150 (MQ=255)
               aGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGACCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagc                                                                                                                                      >  1:68758/1‑150 (MQ=255)
                      gcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCa                                                                                                                               >  1:420991/1‑150 (MQ=255)
                         gATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAACTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCt                                                                                                                            <  1:1675336/150‑1 (MQ=255)
                          aTCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCtt                                                                                                                           <  2:1010406/150‑1 (MQ=255)
                           tCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTa                                                                                                                          <  2:750031/150‑1 (MQ=255)
                            cccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTAt                                                                                                                         >  1:1514550/1‑150 (MQ=255)
                             ccGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATc                                                                                                                        >  1:1196252/1‑150 (MQ=255)
                                 aTAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTcc                                                                                                                    <  1:1020811/150‑1 (MQ=255)
                                 aTAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTcc                                                                                                                    <  1:1173359/150‑1 (MQ=255)
                                   aaaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAg                                                                                                                  >  2:390295/1‑150 (MQ=255)
                                    aaCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCACGCACCAGCAGAATCAGCTTATCTTCCAGc                                                                                                                 >  2:1246824/1‑150 (MQ=255)
                                            cAGCCAGTTAATGCTCTCCGACAATAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGGTGGATGTTTGACGCGGGCGGTGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                                                                                                         <  2:868138/150‑1 (MQ=255)
                                            cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACACAgcg                                                                                                         >  1:1525027/1‑150 (MQ=255)
                                            cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                                                                                                         >  2:773779/1‑150 (MQ=255)
                                            cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                                                                                                         >  1:159772/1‑150 (MQ=255)
                                            cAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcg                                                                                                         >  2:1379301/1‑150 (MQ=255)
                                             aGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgc                                                                                                        <  2:475054/150‑1 (MQ=255)
                                             aGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAgcgc                                                                                                        >  1:794544/1‑150 (MQ=255)
                                                cAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCt                                                                                                     <  1:556397/150‑1 (MQ=255)
                                                  gTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAAGCGGGCGGAGTCAAGCAACAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCCCTTc                                                                                                   >  1:1237605/1‑150 (MQ=255)
                                                   ttAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTcc                                                                                                  <  1:1111398/150‑1 (MQ=255)
                                                   ttAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTcc                                                                                                  >  2:940195/1‑150 (MQ=255)
                                                   ttAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTcc                                                                                                  <  1:810513/150‑1 (MQ=255)
                                                    tAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTccc                                                                                                 >  1:1120976/1‑149 (MQ=255)
                                                         cACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGcc                                                                                            <  1:1292406/150‑1 (MQ=255)
                                                              cGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGcc                                                                                       <  1:360315/150‑1 (MQ=255)
                                                                aCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGa                                                                                     >  2:1441073/1‑150 (MQ=255)
                                                                 cAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGAt                                                                                    >  2:295871/1‑150 (MQ=255)
                                                                  aaGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATg                                                                                   <  1:197272/150‑1 (MQ=255)
                                                                    gAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGGTAACAGCGCAGAACTCAAGCAGTATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTc                                                                                 <  2:1268458/150‑1 (MQ=255)
                                                                        ttAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGAt                                                                             <  1:1617461/150‑1 (MQ=255)
                                                                         tAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCCTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACGAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATa                                                                            <  1:1173443/150‑1 (MQ=255)
                                                                          aGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGAtta                                                                           >  2:241251/1‑148 (MQ=255)
                                                                            ccccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATaaaa                                                                         >  2:556516/1‑150 (MQ=255)
                                                                             cccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAc                                                                        >  2:1682219/1‑150 (MQ=255)
                                                                             cccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAc                                                                        >  2:1633345/1‑150 (MQ=255)
                                                                              ccGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAAcc                                                                       <  1:701011/150‑1 (MQ=255)
                                                                                   aaTGGTGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTg                                                                  <  2:1582339/150‑1 (MQ=255)
                                                                                   aaTGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTg                                                                  <  2:1151943/150‑1 (MQ=255)
                                                                                    aTGGGGCAGTTAATGTAAATTACGCGTTAACCGCGCAGAACGCACGCGGGTTGTTTTACGCTGGAGGCGTCAAGCACCAGCATATTCAGATTTTCTTCCCGCACAAAGCTCGCTTCCATAGACTAGCCGATGTCATATAAAACCTGTTgt                                                                 <  2:915231/150‑2 (MQ=255)
                                                                                        ggCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACt                                                             <  1:307306/150‑1 (MQ=255)
                                                                                           aTTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCaa                                                          >  2:681038/1‑150 (MQ=255)
                                                                                             tGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAgg                                                        <  2:895901/150‑1 (MQ=255)
                                                                                              taaTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                       <  2:1318741/149‑1 (MQ=255)
                                                                                              gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAATCCTGTTGACACTCAAGGt                                                       >  1:683959/1‑150 (MQ=255)
                                                                                              gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                       >  2:801025/1‑150 (MQ=255)
                                                                                              gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                       >  2:463985/1‑150 (MQ=255)
                                                                                              gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                       <  2:882872/150‑1 (MQ=255)
                                                                                              gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                       >  2:147944/1‑150 (MQ=255)
                                                                                              gAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGt                                                       <  2:572642/150‑1 (MQ=255)
                                                                                               aaTGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTa                                                      >  2:1334896/1‑150 (MQ=255)
                                                                                                 tGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAAt                                                    <  2:1366703/150‑1 (MQ=255)
                                                                                                      aTTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCa                                               <  2:243179/150‑1 (MQ=255)
                                                                                                      aTTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCa                                               >  1:336447/1‑150 (MQ=255)
                                                                                                         aCGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGa                                            >  2:1143527/1‑150 (MQ=255)
                                                                                                         aCGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGa                                            >  1:380012/1‑150 (MQ=255)
                                                                                                          cgcgTTAACAGCGCTGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCCGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATgct                                           >  1:1254705/1‑148 (MQ=255)
                                                                                                          cgcgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGAt                                           >  1:629183/1‑150 (MQ=255)
                                                                                                            cgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCg                                         <  1:984109/150‑1 (MQ=255)
                                                                                                            cgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCg                                         >  1:1495475/1‑150 (MQ=255)
                                                                                                            cgTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCg                                         <  1:583198/150‑1 (MQ=255)
                                                                                                                 aCAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAACCTCAAGGTAATTTTCATGATCGATGGc                                    >  2:813285/1‑150 (MQ=255)
                                                                                                                  cAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCTGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCAGTTGAAACTCAAGGTAATTATCAAGATCGATGGCg                                   >  1:1460555/1‑150 (MQ=255)
                                                                                                                    gcgcAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAACGTAATTATCATGATCGATGGCGGt                                 >  1:887450/1‑150 (MQ=255)
                                                                                                                        aGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTcc                             <  2:180222/150‑1 (MQ=255)
                                                                                                                         gAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCa                            >  1:516123/1‑150 (MQ=255)
                                                                                                                           aCTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGa                          <  1:612254/150‑1 (MQ=255)
                                                                                                                             tCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACt                        <  1:173239/150‑1 (MQ=255)
                                                                                                                             tCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGGTTCCAGACt                        <  1:1610943/150‑1 (MQ=255)
                                                                                                                              cAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTg                       >  2:946297/1‑150 (MQ=255)
                                                                                                                              cAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTg                       >  2:398546/1‑150 (MQ=255)
                                                                                                                               aaGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTgg                      <  2:1726098/150‑1 (MQ=255)
                                                                                                                               aaGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAAATATCATGATCGATGGCGGTTTCCAGACTgg                      <  2:840935/150‑1 (MQ=255)
                                                                                                                                 gCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGaa                    >  2:30376/1‑150 (MQ=255)
                                                                                                                                  cAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAAt                   <  2:159772/150‑1 (MQ=255)
                                                                                                                                    ggATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCg                 >  1:1479096/1‑150 (MQ=255)
                                                                                                                                     ggtgcttGACGCCGCCCGCGTCAACCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgt                >  2:808304/6‑150 (MQ=255)
                                                                                                                                      gtgcttGACGCCGCCCGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCTTGATCGATGGCGGTTTCCAGACTGGAATCgta               >  1:1118178/5‑150 (MQ=255)
                                                                                                                                      aTGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCATCCAGCACAAAGCGCGCTTCCAGCGCCGCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgta               <  2:1280341/150‑1 (MQ=255)
                                                                                                                                       tgcttGACGCCGCCCGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtag              >  1:382444/4‑150 (MQ=255)
                                                                                                                                       tGTTTGCCGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGAGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtag              <  1:111930/150‑1 (MQ=255)
                                                                                                                                       tGTTTGACGCGGGCGGTGTCAAGCACCAGCAGAATCAGCTTATCATCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGAAATTATCATGATCGATGGAGGTTTCCAGACTGGAATCgtag              >  1:1567905/1‑150 (MQ=255)
                                                                                                                                       tGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtag              <  2:1400136/150‑1 (MQ=255)
                                                                                                                                       tGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtag              <  1:1032155/150‑1 (MQ=255)
                                                                                                                                        gTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagt             >  1:930058/1‑150 (MQ=255)
                                                                                                                                        gTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagt             >  1:1597462/1‑150 (MQ=255)
                                                                                                                                        gTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagt             <  2:68758/150‑1 (MQ=255)
                                                                                                                                        gTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagt             >  2:665451/1‑150 (MQ=255)
                                                                                                                                        gTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCgtagt             <  1:1380540/150‑1 (MQ=255)
                                                                                                                                              cgcgGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCa       <  1:773779/150‑1 (MQ=255)
                                                                                                                                                 gggcggCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATaa    <  1:1379301/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCAGCACAAAGCGCGCTTCCAGCGCCTCGCCGATGTCAGATAAAACCTGTTGAAACTCAAGGTAATTATCATGATCGATGGCGGTTTCCAGACTGGAATCGTAGTAATCCATAATC  >  NC_000913/4196207‑4196503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: