Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,720,305 G→A 6.0% intergenic (+44/+143) 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,3050GA 6.0% 419.9 / 9.9 149intergenic (+44/+143)cspA/hokARNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base G (62/78);  new base A (4/5);  total (66/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 154 total aligned reads displayed.

ATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCT  >  NC_000913/3720157‑3720454
                                                                                                                                                    |                                                                                                                                                     
  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)
    aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAATGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAt                                                                                                                                                  <  2:1353573/150‑1 (MQ=255)
      aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCGCCAAATGGCAGGGAtctt                                                                                                                                               <  2:898778/150‑1 (MQ=255)
       aCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtctt                                                                                                                                               <  2:1382442/150‑1 (MQ=255)
          aTGGTTACAAATCTCTGGACGAAGGGCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTGAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttttt                                                                                                                                            <  1:1670257/150‑1 (MQ=255)
               tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTg                                                                                                                                       <  1:687463/150‑1 (MQ=255)
                aCAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGt                                                                                                                                      <  1:1584006/150‑1 (MQ=255)
                  aaaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGttt                                                                                                                                    <  1:1304127/150‑1 (MQ=255)
                   aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtttt                                                                                                                                   >  1:1390655/1‑150 (MQ=255)
                    aTCTCTGGACGAAGGTCAGAACGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                  <  2:208585/150‑1 (MQ=255)
                    aTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                                  <  2:513836/150‑1 (MQ=255)
                     tctcTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                                 <  1:1279917/150‑1 (MQ=255)
                        ctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTGTTTATTTGTTTTCAGGa                                                                                                                              >  2:112163/1‑150 (MQ=255)
                          ggACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                            >  2:958353/1‑150 (MQ=255)
                           gACGAAGGTCAGAAAGAGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATATGGCGGGGATTTTTTTATTTGTTTTCAGGaaat                                                                                                                           <  1:287430/150‑1 (MQ=255)
                               aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                                       >  2:561744/1‑150 (MQ=255)
                               aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                                       <  1:1154131/150‑1 (MQ=255)
                               aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                                       <  1:805275/150‑1 (MQ=255)
                               aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                                       <  2:200136/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)
                                      gAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgat                                                                                                                >  1:701610/1‑150 (MQ=255)
                                        aaGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTTGTTACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgatcg                                                                                                              <  2:1159392/150‑1 (MQ=255)
                                        aaGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgatcg                                                                                                              <  2:543871/150‑1 (MQ=255)
                                        aaGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTATTTCCTGAAAacaaataaaaaaatc                                                                                                              <  1:191950/150‑8 (MQ=255)
                                           tgtCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGt                                                                                                           <  1:602886/150‑1 (MQ=255)
                                            gtCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGta                                                                                                          <  2:618848/150‑1 (MQ=255)
                                               cTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaata                                                                                                       >  2:730636/1‑150 (MQ=255)
                                                   aCCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAAt                                                                                                   <  2:855090/150‑1 (MQ=255)
                                                   aCCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATAGCGTAATAAAAt                                                                                                   >  2:445669/1‑150 (MQ=255)
                                                          aaaGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattat                                                                                            <  2:1688123/150‑1 (MQ=255)
                                                            aGCGGCGCTAAAGGCCCGGCATCTGGTAACGGAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGTATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattatta                                                                                          <  1:164779/150‑1 (MQ=255)
                                                                  gcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtg                                                                                    <  2:457706/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:505210/1‑150 (MQ=255)
                                                                                       aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                               >  2:484998/1‑150 (MQ=255)
                                                                                       aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                               >  1:370660/1‑150 (MQ=255)
                                                                                        acgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTATTAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                              >  2:46468/1‑150 (MQ=255)
                                                                                           taacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATga                                                           <  1:739468/150‑1 (MQ=255)
                                                                                             acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                                         >  1:7173/1‑150 (MQ=255)
                                                                                             acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                                         <  2:1182842/150‑1 (MQ=255)
                                                                                             acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                                         <  2:1222176/150‑1 (MQ=255)
                                                                                                aGCCTGTAATCTCTTCTTAAAAGCACAGTATCTTAGATCCCTGGCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGTATAAATTTGGGTGCAATGAGAAtt                                                      <  2:1256257/150‑2 (MQ=255)
                                                                                                 gCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATgc                                                     >  1:1215971/1‑150 (MQ=255)
                                                                                                 gCCTGTAATCTCTGCTTAAAAGCACAGAATCAAAGATCCCTGCCATTTGGCGGGGATTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGAAATGATAATgc                                                     >  1:1336039/1‑150 (MQ=255)
                                                                                                   cTGTAATCTCTTCTTAAAAGCACAGAATCTCAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATgcgc                                                   <  1:334440/150‑1 (MQ=255)
                                                                                                    tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGATATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCa                                                  >  1:1334549/1‑150 (MQ=255)
                                                                                                    tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCa                                                  <  1:484998/150‑1 (MQ=255)
                                                                                                      tAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAAc                                                >  1:1320689/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)
                                                                                                            ctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgta                                          <  2:1313065/150‑1 (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)
                                                                                                                  aaaaGCACAGAATCTAAGATCCCTGCCATTTGGTGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGtaagtaag                                    >  2:879961/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                                 >  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)
                                                                                                                               cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt                       <  2:233618/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)
                                                                                                                                 atgATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTcc                     <  2:480885/148‑1 (MQ=255)
                                                                                                                                 aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTcc                     >  2:1341067/1‑150 (MQ=255)
                                                                                                                                 aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTcc                     >  2:746982/1‑150 (MQ=255)
                                                                                                                                 aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTcc                     <  1:547016/150‑1 (MQ=255)
                                                                                                                                  aaaTCCCCGCCAAATGGCAGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTATTAAAATCTATTTTTATTTTTGTGAAGAATAAATTTGGGTGCTATGAGAATGCGAAACTCCGTAATTAATGAGTGAATAATTTccc                    >  1:566936/3‑150 (MQ=255)
                                                                                                                                   gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccg                   >  1:1412835/1‑150 (MQ=255)
                                                                                                                                     tCCCTGCCATTTGGCGGGGATTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccgcc                 <  2:880969/150‑1 (MQ=255)
                                                                                                                                      cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcct                <  1:1612580/150‑2 (MQ=255)
                                                                                                                                      cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcct                <  1:476845/150‑2 (MQ=255)
                                                                                                                                       ccTGCCATTTGGTGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcctt               <  2:1136143/150‑3 (MQ=255)
                                                                                                                                         tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgccttac             <  1:57902/150‑5 (MQ=255)
                                                                                                                                           ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgccttactt           <  1:366379/150‑7 (MQ=255)
                                                                                                                                                 ggCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGactcttac     <  2:370660/150‑1 (MQ=255)
                                                                                                                                                   cGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGactcttactc   >  1:1018747/1‑150 (MQ=255)
                                                                                                                                                    ggggATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGactcttactct  <  2:1390655/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTCTTACTCT  >  NC_000913/3720157‑3720454

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