Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,720,300 T→A 7.0% intergenic (+39/+148) 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,3000TA 7.0% 428.5 / 14.1 153intergenic (+39/+148)cspA/hokARNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (65/77);  new base A (3/8);  total (68/85)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 156 total aligned reads displayed.

CTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACT  >  NC_000913/3720152‑3720446
                                                                                                                                                    |                                                                                                                                                  
ctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAttt                                                                                                                                                   <  1:561168/150‑1 (MQ=255)
ctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAttt                                                                                                                                                   >  2:687463/1‑150 (MQ=255)
ctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAttt                                                                                                                                                   >  1:855090/1‑150 (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)
        cAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGGAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                           >  2:925916/1‑150 (MQ=255)
         aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAATGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAt                                                                                                                                          <  2:1353573/150‑1 (MQ=255)
           aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCGCCAAATGGCAGGGAtctt                                                                                                                                       <  2:898778/150‑1 (MQ=255)
            aCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtctt                                                                                                                                       <  2:1382442/150‑1 (MQ=255)
              gATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtcttag                                                                                                                                     <  1:11378/150‑3 (MQ=255)
               aTGGTTACAAATCTCTGGACGAAGGGCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTGAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttttt                                                                                                                                    <  1:1670257/150‑1 (MQ=255)
                tGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAAATCCCCGCCAAATGGCAGGGAtcttagat                                                                                                                                   <  2:1137524/150‑5 (MQ=255)
                  gTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAtt                                                                                                                                 <  2:1304410/150‑1 (MQ=255)
                  gTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAtt                                                                                                                                 <  1:516639/150‑1 (MQ=255)
                    tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTg                                                                                                                               <  1:687463/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)
                          tctcTGGAAGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                         <  2:1347928/150‑1 (MQ=255)
                           ctctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAg                                                                                                                        <  2:1429060/150‑1 (MQ=255)
                           ctctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAg                                                                                                                        >  2:1058681/1‑150 (MQ=255)
                            tctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGATAGCCGCGCTAAAGGCCCGGCTGCTGGTATCGTAACCAGCCTGTAATCTGTGCTTAAAAGCACATAATCTAATATCCCTGGCATTTGTCGGGGATTTTTTTATTTGTTTTCAgg                                                                                                                       <  2:498010/150‑1 (MQ=255)
                             ctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTGTTTATTTGTTTTCAGGa                                                                                                                      >  2:112163/1‑150 (MQ=255)
                               ggACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                    >  2:958353/1‑150 (MQ=255)
                               ggACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGATCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAATAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaa                                                                                                                    >  2:1134301/1‑150 (MQ=255)
                                gACGAAGGTCAGAAAGAGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATATGGCGGGGATTTTTTTATTTGTTTTCAGGaaat                                                                                                                   <  1:287430/150‑1 (MQ=255)
                                    aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                               >  2:561744/1‑150 (MQ=255)
                                    aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                               <  2:200136/150‑1 (MQ=255)
                                    aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                               <  1:805275/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)
                                             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:733027/1‑150 (MQ=255)
                                                        aCCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAAt                                                                                           <  2:855090/150‑1 (MQ=255)
                                                        aCCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATAGCGTAATAAAAt                                                                                           >  2:445669/1‑150 (MQ=255)
                                                             cGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAATATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtatt                                                                                      <  2:68566/150‑1 (MQ=255)
                                                              gAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCCGCCTGTAATCTCTGCTTAAAAGCACAGAAGCTAAGATCGCTGCCATTTGGCGGGGATTGTTTTATTTGTTTTCAGGATATAAATAATCGCTCGCGTAAGAAAATCtatta                                                                                     >  1:118506/1‑150 (MQ=255)
                                                               aaaGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattat                                                                                    <  2:1688123/150‑1 (MQ=255)
                                                               aaaGCGGCGCAAAAGGGCCGTCAGCTGGTAACGTAACCATCCTGTAATCTCTGCTTATAATCAAAGAATATAAGATCCCTGCCATTTGGCGGGGTTTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattat                                                                                    <  2:1314823/150‑1 (MQ=255)
                                                                      cgcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAttttt                                                                             >  1:954124/1‑150 (MQ=255)
                                                                       gcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtg                                                                            <  1:1289635/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)
                                                                               cccGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTAAAAAGCACAGAATCTAACATCCCTGCCAGTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGCTCGCGTAATAAAATCTATTATTATTTTTGTgaagatt                                                                    >  2:870023/1‑148 (MQ=255)
                                                                                      gCTGGTAACGTAACCAGCTTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTg                                                             <  2:1140379/150‑1 (MQ=255)
                                                                                        tGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTggg                                                           <  1:45734/150‑1 (MQ=255)
                                                                                         ggTAACGTAACCAGCCTGTAATCACTGCTTAAATGCACAGAATCTAAGATCCCAGCCATTTGGCGGGGATTATATAATTTGTTATCAGGAAATAAATAATCGATCGCGTAATAAAATATATTATTATTTTTGTGAAGAATAAAATTggct                                                          >  1:116447/1‑148 (MQ=255)
                                                                                            aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                       >  2:484998/1‑150 (MQ=255)
                                                                                            aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                       >  1:505210/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)
                                                                                                 aacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgag                                                  <  2:434390/150‑1 (MQ=255)
                                                                                                  acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                                 <  2:1222176/150‑1 (MQ=255)
                                                                                                  acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                                 >  1:7173/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)
                                                                                                           tAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAAc                                        >  1:1320689/1‑150 (MQ=255)
                                                                                                             aTCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGc                                      >  1:1621577/1‑150 (MQ=255)
                                                                                                                tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGTATTTTTTTATTTTTTTTCAGGTAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgt                                   <  2:291659/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)
                                                                                                                       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)
                                                                                                                          aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCg                         >  1:491013/1‑150 (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)
                                                                                                                                  aTCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaataat                 <  1:730636/150‑1 (MQ=255)
                                                                                                                                    cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt               <  2:233618/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)
                                                                                                                                       aaaTCCCCGCCAAATGGCAGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTATTAAAATCTATTTTTATTTTTGTGAAGAATAAATTTGGGTGCTATGAGAATGCGAAACTCCGTAATTAATGAGTGAATAATTTccc            >  1:566936/3‑150 (MQ=255)
                                                                                                                                        gATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccg           >  1:1412835/1‑150 (MQ=255)
                                                                                                                                         aTCCCTGCCATTTGGCGTGGATTTTTTTAATTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAAATTGGGTGCAATGAGAATGCGCAACCACGTAAGAAATGCGTGAATAATTTCacgc          >  2:480778/1‑150 (MQ=255)
                                                                                                                                         aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccgc          <  1:405281/150‑1 (MQ=255)
                                                                                                                                          tCCCTGCCATTTGGCGGGGATTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccgcc         <  2:880969/150‑1 (MQ=255)
                                                                                                                                           cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcct        <  1:1612580/150‑2 (MQ=255)
                                                                                                                                           cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcct        <  1:476845/150‑2 (MQ=255)
                                                                                                                                              tGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgccttac     <  1:57902/150‑5 (MQ=255)
                                                                                                                                               gCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGa    >  2:1160325/1‑150 (MQ=255)
                                                                                                                                                ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgccttactt   <  1:360063/150‑7 (MQ=255)
                                                                                                                                                 cATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGact  >  2:355275/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACT  >  NC_000913/3720152‑3720446

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: