Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,720,294 T→C 5.9% intergenic (+33/+154) cspA → / ← hokA RNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,720,2940TC 5.9% 462.5 / 12.3 158intergenic (+33/+154)cspA/hokARNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (66/77);  new base C (3/6);  total (73/85)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 161 total aligned reads displayed.

CACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAG  >  NC_000913/3720145‑3720443
                                                                                                                                                     |                                                                                                                                                     
  cTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGc                                                                                                                                                     >  1:1136143/1‑150 (MQ=255)
  cTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGc                                                                                                                                                     <  1:1624682/150‑1 (MQ=255)
  cTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGc                                                                                                                                                     >  2:176430/1‑150 (MQ=255)
     ctctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAt                                                                                                                                                  >  2:1154131/1‑150 (MQ=255)
       ctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAttt                                                                                                                                                <  1:561168/150‑1 (MQ=255)
        tGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTg                                                                                                                                               <  1:1232069/150‑1 (MQ=255)
         gCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTgg                                                                                                                                              >  1:29629/1‑150 (MQ=255)
          atATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAAATCCCCGCCAAATGGc                                                                                                                                             <  1:127904/149‑1 (MQ=255)
           tATCCAGACCGATGGATACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGca                                                                                                                                            <  1:1296414/150‑2 (MQ=255)
            aTCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGGCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATATCTGATTAAAAGCACAGAATCTAAGATCCCTGCCATTTTGCgg                                                                                                                                           >  1:1570566/1‑150 (MQ=255)
              ccAGAACGATGGGTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCGTCACCAGCGACAGCGGCGCTAAAGGCCCGGCAGCGGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                                                                                                                         <  1:336901/150‑1 (MQ=255)
               cAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGGAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                        >  2:925916/1‑150 (MQ=255)
               cAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                        >  1:506147/1‑150 (MQ=255)
               cAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCGGTAATCTCTGCGTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCggggt                                                                                                                                        >  2:501515/1‑149 (MQ=255)
                  aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCGCCAAATGGCAGGGAtctt                                                                                                                                    <  2:898778/150‑1 (MQ=255)
                   aCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtctt                                                                                                                                    <  2:1382442/150‑1 (MQ=255)
                      aTGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttttt                                                                                                                                 <  2:826988/150‑1 (MQ=255)
                      aTGGTTACAAATCTCTGGACGAAGGGCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTGAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttttt                                                                                                                                 <  1:1670257/150‑1 (MQ=255)
                       tGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTa                                                                                                                                <  1:757614/150‑1 (MQ=255)
                           tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTg                                                                                                                            <  1:687463/150‑1 (MQ=255)
                            aCAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGt                                                                                                                           <  1:1584006/150‑1 (MQ=255)
                            aCAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAGGCCCGGCCGCTGGTAACGGAACCCGCCGGGAATTAATGCTTAAAAGCACAAAATAAAACTTGCCTGCCATTTTGGGGGGTTTTTTTTATTTGtt                                                                                                                          >  1:561966/1‑150 (MQ=255)
                               aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtttt                                                                                                                        >  1:1390655/1‑150 (MQ=255)
                                aTCTCTGGACGAAGGTCAGAACGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                       <  2:208585/150‑1 (MQ=255)
                                aTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                       <  2:513836/150‑1 (MQ=255)
                                 tctcTGGAAGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                      <  2:1347928/150‑1 (MQ=255)
                                    ctGGACGAAGGTCAGAAAGTGTCCATCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGa                                                                                                                   <  1:294000/150‑1 (MQ=255)
                                      ggACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                 >  2:958353/1‑150 (MQ=255)
                                      ggACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTAAGGaaa                                                                                                                 >  1:812826/1‑150 (MQ=255)
                                      ggACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGATCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAATAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                 >  2:1134301/1‑150 (MQ=255)
                                           aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                            >  2:561744/1‑150 (MQ=255)
                                           aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                            <  1:1154131/150‑1 (MQ=255)
                                              gTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAaataat                                                                                                         >  2:57902/1‑150 (MQ=255)
                                                cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGAACCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcg                                                                                                       <  1:246123/150‑1 (MQ=255)
                                                cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCAGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcg                                                                                                       >  2:1529057/1‑150 (MQ=255)
                                                    aaGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgatcg                                                                                                   <  2:543871/150‑1 (MQ=255)
                                                    aaGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTATTTCCTGAAAacaaataaaaaaatc                                                                                                   <  1:191950/150‑8 (MQ=255)
                                                           cTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaata                                                                                            >  2:730636/1‑150 (MQ=255)
                                                             tCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATaaa                                                                                          >  2:221991/1‑150 (MQ=255)
                                                              cACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATaaaa                                                                                         >  1:602034/1‑150 (MQ=255)
                                                                 cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCt                                                                                      >  2:1723060/1‑150 (MQ=255)
                                                                    cGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAATATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtatt                                                                                   <  2:68566/150‑1 (MQ=255)
                                                                     gAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCCGCCTGTAATCTCTGCTTAAAAGCACAGAAGCTAAGATCGCTGCCATTTGGCGGGGATTGTTTTATTTGTTTTCAGGATATAAATAATCGCTCGCGTAAGAAAATCtatta                                                                                  >  1:118506/1‑150 (MQ=255)
                                                                      aaaGCGGCGCAAAAGGGCCGTCAGCTGGTAACGTAACCATCCTGTAATCTCTGCTTATAATCAAAGAATATAAGATCCCTGCCATTTGGCGGGGTTTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattat                                                                                 <  2:1314823/150‑1 (MQ=255)
                                                                              gcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtg                                                                         <  1:1289635/150‑1 (MQ=255)
                                                                                tAAAGGCCCGTCAGCTGGTAAAGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgtg                                                                       <  2:874530/150‑1 (MQ=255)
                                                                                tAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCTGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgtg                                                                       <  1:1493946/150‑1 (MQ=255)
                                                                                  aaGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTgaa                                                                     <  1:527950/150‑1 (MQ=255)
                                                                                      cccGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAAt                                                                 >  1:615565/1‑150 (MQ=255)
                                                                                      cccGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAAt                                                                 <  2:29629/150‑1 (MQ=255)
                                                                                      cccGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTAAAAAGCACAGAATCTAACATCCCTGCCAGTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGCTCGCGTAATAAAATCTATTATTATTTTTGTgaagatt                                                                 >  2:870023/1‑148 (MQ=255)
                                                                                           cAGCTGGTAACGTAACCAGCATGTAATCTCTGCTAAAAAGCACAGAATCTAAGATCCCTGCCATTTGGAGGTGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAAtt                                                            <  1:416844/150‑1 (MQ=255)
                                                                                             gCTGGTAACGTAACCAGCTTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                          <  2:1140379/150‑1 (MQ=255)
                                                                                               ttgtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTggg                                                        <  1:925916/148‑1 (MQ=255)
                                                                                                ggTAACGTAACCAGCCTGTAATCACTGCTTAAATGCACAGAATCTAAGATCCCAGCCATTTGGCGGGGATTATATAATTTGTTATCAGGAAATAAATAATCGATCGCGTAATAAAATATATTATTATTTTTGTGAAGAATAAAATTggct                                                       >  1:116447/1‑148 (MQ=255)
                                                                                                   aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                    >  1:370660/1‑150 (MQ=255)
                                                                                                   aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                    >  1:505210/1‑150 (MQ=255)
                                                                                                   aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                    >  2:484998/1‑150 (MQ=255)
                                                                                                    acgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTATTAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                   >  2:46468/1‑150 (MQ=255)
                                                                                                        aacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgag                                               <  2:434390/150‑1 (MQ=255)
                                                                                                         acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                              <  2:1222176/150‑1 (MQ=255)
                                                                                                            aGCCTGTAATCTCTTCTTAAAAGCACAGTATCTTAGATCCCTGGCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGTATAAATTTGGGTGCAATGAGAAtt                                           <  2:1256257/150‑2 (MQ=255)
                                                                                                             gCCTGTAATCTCTGCTTAAAAGCACAGAATCAAAGATCCCTGCCATTTGGCGGGGATTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGAAATGATAATgc                                          >  1:1336039/1‑150 (MQ=255)
                                                                                                                tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGATATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCa                                       >  1:1334549/1‑150 (MQ=255)
                                                                                                                tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCa                                       <  1:484998/150‑1 (MQ=255)
                                                                                                                  tAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAAc                                     >  1:1320689/1‑150 (MQ=255)
                                                                                                                    aTCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGc                                   >  1:1621577/1‑150 (MQ=255)
                                                                                                                       tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGTATTTTTTTATTTTTTTTCAGGTAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgt                                <  2:291659/150‑1 (MQ=255)
                                                                                                                       tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATTCGCAACGCCgt                                >  2:1078165/1‑150 (MQ=255)
                                                                                                                       tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgt                                <  2:1054636/150‑1 (MQ=255)
                                                                                                                        ctGCTTAAAAGCACAGAATCTAAGATCCCAGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATACATTTGGGTGCAATGAGAATGCGCAACGCCgta                               <  1:70411/150‑1 (MQ=255)
                                                                                                                          gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaag                             >  2:1381949/1‑150 (MQ=255)
                                                                                                                          gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaag                             <  2:1621577/150‑1 (MQ=255)
                                                                                                                          gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGAAATGAGAATGCGCAACGCCgtaag                             >  2:1176230/1‑150 (MQ=255)
                                                                                                                           cTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTAGTTTTCAGGAAATAAAAAATCGCTCGCGTAATAAAATCTATTATTATTTTAGTGAAGAATAAATTTGGGTGCAATGATAATTCGCTACGCAgtaagt                            >  1:97990/1‑150 (MQ=255)
                                                                                                                             tAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaagtaa                          >  2:718305/1‑150 (MQ=255)
                                                                                                                              aaaaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGATTAAATTTGGGTGCAATGAGAATGCGCAACGCCGtaagtaag                         >  2:274409/1‑150 (MQ=255)
                                                                                                                              aaaaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGtaagtaag                         >  2:1212681/1‑150 (MQ=255)
                                                                                                                               aaaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTTTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaagtaatg                        >  2:823917/1‑148 (MQ=255)
                                                                                                                               aaaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAgg                        <  2:634527/150‑1 (MQ=255)
                                                                                                                                 aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCg                      <  2:450433/150‑1 (MQ=255)
                                                                                                                                 aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCg                      >  1:782998/1‑150 (MQ=255)
                                                                                                                                   cacaGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCggg                    <  1:208919/150‑1 (MQ=255)
                                                                                                                                   cacaGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAAAGCCGTAAGAAAGGCggg                    >  2:1312214/1‑150 (MQ=255)
                                                                                                                                   cacaGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGAAACGCCGTAAGTAAGGCggg                    >  2:1254928/1‑150 (MQ=255)
                                                                                                                                      aGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaat                 <  1:733027/150‑1 (MQ=255)
                                                                                                                                         aTCTAAGATCCCTGCCATTTGGTGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaataat              >  2:142932/1‑150 (MQ=255)
                                                                                                                                         aTCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaataat              <  1:730636/150‑1 (MQ=255)
                                                                                                                                           cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt            <  2:445060/150‑1 (MQ=255)
                                                                                                                                           aaaaaaTCCCCGCCAATTGGCAGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt            >  1:91560/6‑150 (MQ=255)
                                                                                                                                           aaaaaaTCCCAGCCAAAAGGCAGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATACGCAACGCCGTAAGTAAGGCGGGAATAAttt            >  2:132123/6‑150 (MQ=255)
                                                                                                                                            aaaaaTCCCCGCCAAATGGCAGGGATTTTTTTATGTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATATTTTc           >  1:1048514/5‑150 (MQ=255)
                                                                                                                                             aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTcc          >  2:746982/1‑150 (MQ=255)
                                                                                                                                               gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccg        >  1:1412835/1‑150 (MQ=255)
                                                                                                                                                aTCCCTGCCATTTGGCGTGGATTTTTTTAATTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAAATTGGGTGCAATGAGAATGCGCAACCACGTAAGAAATGCGTGAATAATTTCacgc       >  2:480778/1‑150 (MQ=255)
                                                                                                                                                 tCCCTGCCATTTGGCGGGGATTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccgcc      <  2:880969/150‑1 (MQ=255)
                                                                                                                                                  cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcct     <  1:476845/150‑2 (MQ=255)
                                                                                                                                                   ccTGCCATTTGGTGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcctt    <  2:1136143/150‑3 (MQ=255)
                                                                                                                                                    cTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGaa   >  1:67533/1‑150 (MQ=255)
                                                                                                                                                     tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgccttac  <  1:57902/150‑5 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACTTCTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAG  >  NC_000913/3720145‑3720443

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: