Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,720,299 T→A 7.2% intergenic (+38/+149) 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,2990TA 7.2% 418.5 / 17.1 153intergenic (+38/+149)cspA/hokARNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (63/78);  new base A (4/7);  total (68/85)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 157 total aligned reads displayed.

CTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACT  >  NC_000913/3720150‑3720446
                                                                                                                                                     |                                                                                                                                                   
ctctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAt                                                                                                                                                     >  2:1154131/1‑150 (MQ=255)
  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)
         ccAGAACGATGGGTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCGTCACCAGCGACAGCGGCGCTAAAGGCCCGGCAGCGGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgggg                                                                                                                                            <  1:336901/150‑1 (MQ=255)
          cAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGGAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                           >  2:925916/1‑150 (MQ=255)
          cAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                           >  1:506147/1‑150 (MQ=255)
          cAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGa                                                                                                                                           >  2:1612580/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                                                                                                                                 <  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)
                             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)
                                 ggACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTAAGGaaa                                                                                                                    >  1:812826/1‑150 (MQ=255)
                                  gACGAAGGTCAGAAAGAGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATATGGCGGGGATTTTTTTATTTGTTTTCAGGaaat                                                                                                                   <  1:287430/150‑1 (MQ=255)
                                      aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaat                                                                                                               >  2:561744/1‑150 (MQ=255)
                                         gTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAaataat                                                                                                            >  2:57902/1‑150 (MQ=255)
                                           cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGAACCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcg                                                                                                          <  1:246123/150‑1 (MQ=255)
                                           cAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCAGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcg                                                                                                          >  2:1529057/1‑150 (MQ=255)
                                               aaGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTTGTTACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgatcg                                                                                                      <  2:1159392/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)
                                                            cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCt                                                                                         >  2:1723060/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)
                                                                   aGCGGCGCTAAAGGCCCGGCATCTGGTAACGGAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGTATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCtattatta                                                                                  <  1:164779/150‑1 (MQ=255)
                                                                     cggcgCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGTGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTattattatt                                                                                <  2:525917/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)
                                                                                 cccGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAAt                                                                    <  2:29629/150‑1 (MQ=255)
                                                                                 cccGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTAAAAAGCACAGAATCTAACATCCCTGCCAGTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGCTCGCGTAATAAAATCTATTATTATTTTTGTgaagatt                                                                    >  2:870023/1‑148 (MQ=255)
                                                                                          ttgtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTggg                                                           <  1:925916/148‑1 (MQ=255)
                                                                                          tGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTggg                                                           <  1:45734/150‑1 (MQ=255)
                                                                                           ggTAACGTAACCAGCCTGTAATCACTGCTTAAATGCACAGAATCTAAGATCCCAGCCATTTGGCGGGGATTATATAATTTGTTATCAGGAAATAAATAATCGATCGCGTAATAAAATATATTATTATTTTTGTGAAGAATAAAATTggct                                                          >  1:116447/1‑148 (MQ=255)
                                                                                              aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                       >  1:370660/1‑150 (MQ=255)
                                                                                              aacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCa                                                       >  2:484998/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                                                 >  1:7173/1‑150 (MQ=255)
                                                                                                    acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                                 <  2:1222176/150‑1 (MQ=255)
                                                                                                    acCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATgaga                                                 <  2:1182842/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)
                                                                                                             tAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAAc                                        >  1:1320689/1‑150 (MQ=255)
                                                                                                                  tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGTATTTTTTTATTTTTTTTCAGGTAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgt                                   <  2:291659/150‑1 (MQ=255)
                                                                                                                  tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgt                                   <  2:1054636/150‑1 (MQ=255)
                                                                                                                   ctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgta                                  <  2:1313065/150‑1 (MQ=255)
                                                                                                                   ctGCTTAAAAGCACAGAATCTAAGATCCCAGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATACATTTGGGTGCAATGAGAATGCGCAACGCCgta                                  <  1:70411/150‑1 (MQ=255)
                                                                                                                     gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGAAATGAGAATGCGCAACGCCgtaag                                >  2:1176230/1‑150 (MQ=255)
                                                                                                                         aaaaGCACAGAATCTAAGATCCCTGCCATTTGGTGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGtaagtaag                            >  2:879961/1‑150 (MQ=255)
                                                                                                                         aaaaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGATTAAATTTGGGTGCAATGAGAATGCGCAACGCCGtaagtaag                            >  2:274409/1‑150 (MQ=255)
                                                                                                                          aaaGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTTTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaagtaatg                           >  2:823917/1‑148 (MQ=255)
                                                                                                                            aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCg                         <  2:450433/150‑1 (MQ=255)
                                                                                                                            aGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCg                         >  1:491013/1‑150 (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)
                                                                                                                                      cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt               <  2:233618/150‑1 (MQ=255)
                                                                                                                                      cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAttt               <  2:445060/150‑1 (MQ=255)
                                                                                                                                      aaaaaaTCCCAGCCAAAAGGCAGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATACGCAACGCCGTAAGTAAGGCGGGAATAAttt               >  2:132123/6‑150 (MQ=255)
                                                                                                                                       aaaaaTCCCCGCCAAATGGCAGGGATTTTTTTATGTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATATTTTc              >  1:1048514/5‑150 (MQ=255)
                                                                                                                                        aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTcc             >  2:1341067/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)
                                                                                                                                           aTCCCTGCCATTTGGCGTGGATTTTTTTAATTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAAATTGGGTGCAATGAGAATGCGCAACCACGTAAGAAATGCGTGAATAATTTCacgc          >  2:480778/1‑150 (MQ=255)
                                                                                                                                           aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccgc          <  1:405281/150‑1 (MQ=255)
                                                                                                                                            tCCCTGCCATTTGGCGGGGATTTTTTTATATGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccgcc         <  2:880969/150‑1 (MQ=255)
                                                                                                                                             cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcct        <  1:476845/150‑2 (MQ=255)
                                                                                                                                             cccTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgcct        <  1:1612580/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)
                                                                                                                                                  ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAAATTATTCccgccttactt   <  1:360063/150‑7 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACT  >  NC_000913/3720150‑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: