Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,068,922 G→T 7.6% intergenic (‑114/+25) argO ← / ← mscS L‑arginine exporter/small conductance mechanosensitive channel MscS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,068,9220GT 7.6% 270.7 / 8.3 83intergenic (‑114/+25)argO/mscSL‑arginine exporter/small conductance mechanosensitive channel MscS
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base T (3/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCA  >  NC_000913/3068773‑3069071
                                                                                                                                                     |                                                                                                                                                     
cccAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCagag                                                                                                                                                       >  2:1816566/1‑150 (MQ=255)
       gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGTCCAGAGAGtatta                                                                                                                                                <  2:1689592/150‑1 (MQ=255)
        cAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtattat                                                                                                                                               >  2:1035734/1‑150 (MQ=255)
          aGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCTAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                                                                             <  1:1725197/150‑1 (MQ=255)
          aGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                                                                             <  1:1864173/150‑1 (MQ=255)
           gACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTg                                                                                                                                            >  1:1755058/1‑150 (MQ=255)
             ccTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                                                                          >  2:190760/1‑150 (MQ=255)
             ccTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                                                                          >  2:1678334/1‑150 (MQ=255)
             ccTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                                                                          >  2:1675072/1‑150 (MQ=255)
              cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGt                                                                                                                                                    <  2:1055880/139‑1 (MQ=255)
              cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGt                                                                                                                                                    >  1:1055880/1‑139 (MQ=255)
               ttGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                                                                        <  1:1771776/150‑1 (MQ=255)
               ttGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                                                                        >  2:246110/1‑150 (MQ=255)
                tGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTc                                                                                                                                       <  1:1708741/150‑1 (MQ=255)
                tGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTc                                                                                                                                       >  2:1102596/1‑150 (MQ=255)
                  aaaGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTctc                                                                                                                                     <  2:496360/150‑3 (MQ=255)
                      taataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGtt                                                                                                                                 >  1:401405/1‑150 (MQ=255)
                        ataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGa                                                                                                                               <  2:1755058/150‑1 (MQ=255)
                          aaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                                                                             >  2:1435570/1‑150 (MQ=255)
                              aaaCACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGc                                                                                                                         <  2:131190/150‑1 (MQ=255)
                                 cacGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGc                                                                                                                      <  2:144731/150‑1 (MQ=255)
                                   cGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCtt                                                                                                                    <  1:891635/150‑1 (MQ=255)
                                     ttGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGCGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTg                                                                                                                  <  1:1816566/150‑1 (MQ=255)
                                     ttGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTg                                                                                                                  <  1:190760/150‑1 (MQ=255)
                                            ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                           >  1:916201/1‑150 (MQ=255)
                                            ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                           >  2:1684792/1‑150 (MQ=255)
                                            ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                           <  2:1351626/150‑1 (MQ=255)
                                             tCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCttt                                                                                                          <  1:66013/150‑1 (MQ=255)
                                                 aaTTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCAcc                                                                                                      >  1:593500/1‑150 (MQ=255)
                                                     aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCt                                                                                                  <  1:246110/150‑1 (MQ=255)
                                                     aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCt                                                                                                  <  1:1884490/150‑1 (MQ=255)
                                                      aGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCtt                                                                                                 >  2:1349409/1‑150 (MQ=255)
                                                        tATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaa                                                                                               <  1:577905/150‑1 (MQ=255)
                                                        tATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaa                                                                                               <  1:1435570/150‑1 (MQ=255)
                                                             gATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTc                                                                                          <  2:1018576/150‑1 (MQ=255)
                                                             gATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTc                                                                                          <  2:1714599/150‑1 (MQ=255)
                                                               tGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGACGCCAGAGAGTATTATCTGGCCTCAGAGTTGATTACGCAGCTTTGTCTTCTTTCACACGCTTAAAGTTcac                                                                                        >  1:904083/1‑150 (MQ=255)
                                                                gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                                                                       <  1:247896/150‑1 (MQ=255)
                                                                gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                                                                       >  1:28212/1‑150 (MQ=255)
                                                                 cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAAATAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACAt                                                                                      >  2:1692140/1‑150 (MQ=255)
                                                                      aCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCAtt                                                                                 <  1:450430/150‑1 (MQ=255)
                                                                      aCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCAtt                                                                                 <  1:1675072/150‑1 (MQ=255)
                                                                            tCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCggg                                                                           <  1:982740/150‑1 (MQ=255)
                                                                                gggATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACg                                                                       <  1:1102596/150‑1 (MQ=255)
                                                                                 ggATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACgg                                                                      >  2:1839246/1‑150 (MQ=255)
                                                                                 ggATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACgg                                                                      <  2:47450/150‑1 (MQ=255)
                                                                                    ttGTTATTAGCGGAAATTGATAATATTAATTGGCGATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGaa                                                                   <  1:200704/150‑1 (MQ=255)
                                                                                            aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATAc                                                           <  2:546115/150‑1 (MQ=255)
                                                                                                  aaTTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggcgg                                                     <  2:1213725/150‑1 (MQ=255)
                                                                                                    ttGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCa                                                   >  2:1249740/1‑150 (MQ=255)
                                                                                                           tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAtt                                            >  2:176085/1‑150 (MQ=255)
                                                                                                           tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAtt                                            <  1:1684792/150‑1 (MQ=255)
                                                                                                            atTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTc                                           <  1:876212/150‑1 (MQ=255)
                                                                                                            atTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTc                                           >  2:896494/1‑150 (MQ=255)
                                                                                                             ttAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCa                                          >  1:1003658/1‑150 (MQ=255)
                                                                                                              tAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCAc                                         >  2:193697/1‑150 (MQ=255)
                                                                                                                 ttGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGtt                                      >  1:476395/1‑150 (MQ=255)
                                                                                                                  tGGCAATCAGAAAACCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTCCGCAGCTTTGTGTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGttt                                     <  2:1754357/150‑1 (MQ=255)
                                                                                                                   ggCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTa                                    >  1:1064318/1‑150 (MQ=255)
                                                                                                                    gCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTaa                                   >  1:808884/1‑150 (MQ=255)
                                                                                                                      aaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGa                                                              >  1:1490842/1‑121 (MQ=255)
                                                                                                                      aaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGa                                                              <  2:1490842/121‑1 (MQ=255)
                                                                                                                      aaTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATa                                 <  2:916201/150‑1 (MQ=255)
                                                                                                                          aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCATCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCt                             >  2:1934868/1‑150 (MQ=255)
                                                                                                                          aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCt                             >  1:1029581/1‑150 (MQ=255)
                                                                                                                           gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAAt                                  <  1:187145/144‑1 (MQ=255)
                                                                                                                           gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAAt                                  >  2:187145/1‑144 (MQ=255)
                                                                                                                             aaaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCa                          <  2:593500/150‑1 (MQ=255)
                                                                                                                              aaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAg                         >  1:80565/1‑150 (MQ=255)
                                                                                                                                aCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGca                       <  2:904083/150‑1 (MQ=255)
                                                                                                                                     taaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATcc                  >  2:1985242/1‑150 (MQ=255)
                                                                                                                                      aaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATccc                 >  2:1309190/1‑150 (MQ=255)
                                                                                                                                      aaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCAATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATccc                 >  2:64579/1‑150 (MQ=255)
                                                                                                                                       agagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCa                >  1:968708/1‑150 (MQ=255)
                                                                                                                                        gagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAg               >  1:21839/1‑150 (MQ=255)
                                                                                                                                         gctgagGCCAGATAATACGCTCTGGCCTCAGCGGTGATTAAGCAGCTTTGTCTTATTTCACCCGCTGAAAGTTCACATCCATTTGCGGGGACGGGAAGTTGATACCGGCGGCAGCAAATTCACGGGTAATACGCTCCCGCACATCCCAgg              >  2:394816/4‑149 (MQ=255)
                                                                                                                                         agagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGt              >  1:1037130/1‑150 (MQ=255)
                                                                                                                                          gagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGAATACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTATAATACGCTCCAGCACATCCCAGTa             >  2:4136/1‑150 (MQ=255)
                                                                                                                                          ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTa             >  1:885508/3‑150 (MQ=255)
                                                                                                                                            gagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTaca           >  2:1381629/1‑150 (MQ=255)
                                                                                                                                            gagGCCAGAAAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACCTCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCCAATTCACGTTTAATACGCTCCAGCACAGCCCAGTacc           >  1:863645/1‑149 (MQ=255)
                                                                                                                                              ggCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACg         <  1:293387/150‑1 (MQ=255)
                                                                                                                                                ccAGAGAGTATTATCTGGCCTCAGCGTTGATAACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGtt       >  1:315399/1‑150 (MQ=255)
                                                                                                                                                 cAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGttt      <  1:1839246/150‑1 (MQ=255)
                                                                                                                                                  agagagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACTTCCCAGTACACGtttt     <  1:1035734/150‑1 (MQ=255)
                                                                                                                                                     gagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCa  >  1:199583/1‑150 (MQ=255)
                                                                                                                                                     gagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCa  <  1:1249740/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCA  >  NC_000913/3068773‑3069071

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