Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,967,036 T→A 6.3% L207Q (CTG→CAG)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,0360TA 6.3% 362.2 / 8.2 120L207Q (CTG→CAG) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (53/59);  new base A (5/3);  total (58/62)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

AACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACA  >  NC_000913/3966887‑3967185
                                                                                                                                                     |                                                                                                                                                     
aaCTGCTCGCGTTCTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTct                                                                                                                                                       <  2:1436182/150‑1 (MQ=255)
aaCAGCTCGCGTACAGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTct                                                                                                                                                       <  2:1734960/150‑1 (MQ=255)
  cTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTctgc                                                                                                                                                     <  1:681901/150‑1 (MQ=255)
   tGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTctgct                                                                                                                                                    <  1:1962692/150‑1 (MQ=255)
        gcgTACTGTATCTGCCATCACCTTTCGGTCGTGGTCATCGTGGTATGATTGTGGCACCGCAGAAAGCCGTTAAAACCATGCTGCTTCAGAACATTGCGCAGAGCTTTGCTTACAACCACCCGGATTGTGTGCCGATGTTTCTGGTGATcg                                                                                                                                               <  2:780108/150‑1 (MQ=255)
        gcgTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGCGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcg                                                                                                                                               >  1:281140/1‑150 (MQ=255)
          gTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcgac                                                                                                                                             >  1:30371/1‑150 (MQ=255)
          gTACGGGATCTGGCATCACCTATCGGTCGGGGTCAGCGGGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATcgac                                                                                                                                             <  2:842406/150‑1 (MQ=255)
             cTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGaa                                                                                                                                          <  1:942062/150‑1 (MQ=255)
               ggATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCCCCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACg                                                                                                                                        >  2:250321/1‑150 (MQ=255)
                  tCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTcc                                                                                                                                     >  1:2089878/1‑150 (MQ=255)
                  tCTGGCATCACCTATCGGTCGGGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGGcc                                                                                                                                     <  2:1853074/150‑1 (MQ=255)
                       cATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGgaag                                                                                                                                >  1:236287/1‑150 (MQ=255)
                         tCACCTATCGGTCGTGGTAAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGgaagaa                                                                                                                              >  1:1538136/1‑150 (MQ=255)
                          aaCCTATCGGTCGTGGTCAGCTTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGTTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGaagaag                                                                                                                             <  1:1059533/149‑1 (MQ=255)
                          aaCCTATCGGTCGGGGTCAGCGTGGTCTGATTGTGGCACCGCCGAATGCCGGTAAAACCATGCTGGTGCAGAAGATTGCTCAGAGCATTGCTTACTACCACCCTGATTGTTTGCTGATGGTTCTGCGGATCGACGAACGTCCGGaagaag                                                                                                                             <  1:924487/149‑1 (MQ=255)
                           aCCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGt                                                                                                                            <  2:1567305/150‑1 (MQ=255)
                             cTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCTGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCCGAAGAAGTaa                                                                                                                          >  1:1451736/1‑150 (MQ=255)
                              tATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAAc                                                                                                                         >  1:1237830/1‑150 (MQ=255)
                              tATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGGGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAAc                                                                                                                         <  1:1590613/150‑1 (MQ=255)
                              tATCGGTCGTGGGCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCAGCTGATCGACGAACGTCCGGAAGAAGTAAc                                                                                                                         <  2:1941202/150‑1 (MQ=255)
                                  ggTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAACCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCgag                                                                                                                     >  1:432822/1‑150 (MQ=255)
                                  ggTCGTGCGCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCgag                                                                                                                     <  1:1900753/150‑1 (MQ=255)
                                  ggTAGTGGTCAGCGTGGTCTGATTGTGGCACCGACGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTATGCTGATCGACGAAAGTACGGAAGAAGTAAACgat                                                                                                                     >  1:1356324/1‑149 (MQ=255)
                                    gcGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGAt                                                                                                                   <  2:1935444/149‑1 (MQ=255)
                                     ggTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACTTTTCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCAGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATg                                                                                                                  <  1:271767/149‑1 (MQ=255)
                                     cGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATg                                                                                                                  >  2:1492662/1‑150 (MQ=255)
                                      gTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCGCAGAGCATTGCTTACAACCCCCCGGATTGTGTGCTGATGGTGCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGc                                                                                                                 <  2:351530/150‑1 (MQ=255)
                                        ggTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCATAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAg                                                                                                               <  2:1599780/150‑1 (MQ=255)
                                         gTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGc                                                                                                              >  2:1131791/1‑150 (MQ=255)
                                          tCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCg                                                                                                             <  2:1648435/150‑1 (MQ=255)
                                          tCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTACTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCg                                                                                                             <  1:1236826/150‑1 (MQ=255)
                                           cAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCAGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGt                                                                                                            >  2:2020697/1‑150 (MQ=255)
                                             gCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCt                                                                                                          >  1:1007017/1‑150 (MQ=255)
                                             gCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCt                                                                                                          <  2:2089878/150‑1 (MQ=255)
                                                 ggTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGCTGCAGCGTCTGGTa                                                                                                      >  2:391593/1‑150 (MQ=255)
                                                    cTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGATATGCAGCGTCTGGTaaaa                                                                                                   <  1:462756/150‑1 (MQ=255)
                                                        ttGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGATGATGGTTCTGCTGATCGACGCACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTg                                                                                               >  2:1243859/1‑150 (MQ=255)
                                                         tgtgGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATAGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGa                                                                                              <  1:1922757/150‑1 (MQ=255)
                                                             gCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCCCCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgtt                                                                                          >  1:733725/1‑150 (MQ=255)
                                                              cACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgttg                                                                                         >  1:1596259/1‑150 (MQ=255)
                                                              cACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCAGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAgttg                                                                                         >  1:1555366/1‑150 (MQ=255)
                                                                 cgccgAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACCCCCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGGCTGGTAAAAGGTGAAGttgttg                                                                                      >  1:1420641/1‑150 (MQ=255)
                                                                 cgccgAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGttgttg                                                                                      <  1:1883396/150‑1 (MQ=255)
                                                                     gtaaGCCGGTAAGACCTTGCTTCTGCAGAACATTGCTCAGAGCATTGCTTACCACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTTACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTc                                                                                  <  2:274219/148‑1 (MQ=255)
                                                                     gAAAGCCGTTAAAACATTGCTGCTGCAGAAAATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTc                                                                                  <  1:996989/150‑1 (MQ=255)
                                                                      aaaGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGGCTGGTAAAAGGTGAAGTTGTTGCTTCt                                                                                 <  2:1407335/150‑1 (MQ=255)
                                                                       aaGCCGGTAAATCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACATCCACCCGGAATGTGTGCTGATGCTTCTGCTGATCGCCGAACGACCGCAAGAAGTAACCGAGATGCAGCGTCTGTAAAAGGTGAAGTTGTTGCTTCTAc                                                                               >  2:424751/1‑150 (MQ=255)
                                                                         gCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGAGTGAGTGCTGATGGTTCTGCTGATCGACGATCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGATGTAGCATCTAcc                                                                              >  2:1723425/1‑150 (MQ=255)
                                                                            ggTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCttt                                                                           >  1:1818315/1‑150 (MQ=255)
                                                                               aaaaCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGAc                                                                        <  1:223883/150‑1 (MQ=255)
                                                                                aaaCCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACg                                                                       <  1:952215/150‑1 (MQ=255)
                                                                                       gctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGc                                                                <  2:1853772/150‑1 (MQ=255)
                                                                                        ctgctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGGGCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCa                                                               <  2:243571/150‑1 (MQ=255)
                                                                                          gctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCAtc                                                             <  1:1040400/150‑1 (MQ=255)
                                                                                          gctgcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCAtc                                                             >  1:1433605/1‑150 (MQ=255)
                                                                                             gcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCg                                                          >  1:1767018/1‑150 (MQ=255)
                                                                                             gcAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCAGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCg                                                          >  2:1483160/1‑150 (MQ=255)
                                                                                             gcAGAACATTGCTCAGAGCATTGCTTACAACAAACCGGATTGAGTGCTGATGGTTCTGCTGATCGAAGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTATAAGGTGAAGTTGTTGCTTCTACCTTTGACTAACAAGCAtctag                                                          >  2:854325/1‑148 (MQ=255)
                                                                                             gcAGAACATTGCTCAGAGCATTGCGTACAACCACCCGGATTGTGGGCTGATGGTTGTGCTGAGCGTCGAACGTCCGGACGAAGTACCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCg                                                          <  1:1257198/150‑1 (MQ=255)
                                                                                                gAACATTGCTCTGAGCATTGCTTACAACCACCCGGTTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAAGCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCa                                                       <  2:1318363/150‑1 (MQ=255)
                                                                                                gAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCa                                                       >  2:1695822/1‑150 (MQ=255)
                                                                                                   cATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGt                                                    <  2:1243304/150‑1 (MQ=255)
                                                                                                      tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCa                                                 <  1:250321/150‑1 (MQ=255)
                                                                                                       gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAg                                                <  1:740119/150‑1 (MQ=255)
                                                                                                       gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAg                                                >  1:1095287/1‑150 (MQ=255)
                                                                                                        cTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAgg                                               <  2:1917918/150‑1 (MQ=255)
                                                                                                         tCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGt                                              >  1:1638201/1‑150 (MQ=255)
                                                                                                         tCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGGTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGGCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGt                                              <  1:1335875/150‑1 (MQ=255)
                                                                                                          cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGtt                                             >  1:1161679/1‑150 (MQ=255)
                                                                                                           agagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTg                                            <  1:1159689/150‑1 (MQ=255)
                                                                                                            gagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGc                                           >  2:2012374/1‑150 (MQ=255)
                                                                                                              gCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCAGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCgg                                         <  1:2020697/150‑1 (MQ=255)
                                                                                                               cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTTGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGa                                        <  2:2080912/150‑1 (MQ=255)
                                                                                                               cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGGTGCGGa                                        <  1:106636/150‑1 (MQ=255)
                                                                                                                   gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGaaagg                                    >  2:183102/1‑148 (MQ=255)
                                                                                                                   gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATg                                    <  2:264852/150‑1 (MQ=255)
                                                                                                                    cTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATgg                                   >  1:1973725/1‑150 (MQ=255)
                                                                                                                     ttACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGt                                  <  1:1512157/150‑1 (MQ=255)
                                                                                                                      tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTg                                 <  2:1222463/150‑1 (MQ=255)
                                                                                                                       aCAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTACCCGAGATGCAGCGTCTGGTTAAAGGTGAAGTTGTTGCTACTACCTATGACGAACCCGCATCTCGCCACCTTCAGCTTGCGGAAATGGTGa                                >  2:1047503/1‑150 (MQ=255)
                                                                                                                        cAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTACCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTgct                               >  1:657837/1‑148 (MQ=255)
                                                                                                                        cAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGAt                               >  1:1980716/1‑150 (MQ=255)
                                                                                                                              cccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCCCGTTCAGGTTGCGGAAATGGTGATCGAGaa                         >  2:1967437/1‑150 (MQ=255)
                                                                                                                              cccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCCGAAATGGTGATCGAGaa                         >  1:1588417/1‑150 (MQ=255)
                                                                                                                                 ggATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGc                      >  2:1518650/1‑150 (MQ=255)
                                                                                                                                 ggACTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTTGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGc                      <  2:440046/150‑1 (MQ=255)
                                                                                                                                  gATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGcc                     <  2:1007017/150‑1 (MQ=255)
                                                                                                                                   tttttgtgCTGATGTTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCATCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGAATCTCGCCACGTTCAGTTTGCGGAAATGGTGATCGAGAAGGCCa                    <  1:1100856/146‑1 (MQ=255)
                                                                                                                                      gtgtgCTGATGGTTCAGCTGATCGACGAACGTCCGGAAGAAGTACCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAc                 >  1:505646/1‑150 (MQ=255)
                                                                                                                                            tGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTgg           >  1:931690/1‑150 (MQ=255)
                                                                                                                                             gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGt          <  2:658488/150‑1 (MQ=255)
                                                                                                                                               tGGTGCTGCTGAGCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg        <  1:701435/150‑1 (MQ=255)
                                                                                                                                               tGGGTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg        <  1:1254135/150‑1 (MQ=255)
                                                                                                                                                ggTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGa       <  1:1731469/150‑1 (MQ=255)
                                                                                                                                                gcgtCTGCTGATCGACGCACGTCGGGATGAAGTAACCTTGATGCAGCGTCTTGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAATTGGTGATCGAGAAGGCCAAACGCCTGGTTGa       <  1:1723425/147‑1 (MQ=255)
                                                                                                                                                 gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAg      >  1:1331612/1‑150 (MQ=255)
                                                                                                                                                 gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAg      >  2:1519082/1‑150 (MQ=255)
                                                                                                                                                   tCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGca    >  2:1663284/1‑150 (MQ=255)
                                                                                                                                                     tgctgATCGACGAACGTCTGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGCTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca  >  1:1961790/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACA  >  NC_000913/3966887‑3967185

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