Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,720,299 T→A 5.4% 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 5.4% 320.7 / 7.2 92intergenic (+38/+149)cspA/hokARNA chaperone and antiterminator, cold‑inducible/toxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (35/51);  new base A (2/3);  total (38/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

TCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTC  >  NC_000913/3720151‑3720447
                                                                                                                                                    |                                                                                                                                                    
tctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAt                                                                                                                                                      >  2:113377/1‑149 (MQ=255)
tctGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAt                                                                                                                                                      >  4:167031/1‑149 (MQ=255)
  tGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCAttt                                                                                                                                                    <  2:84210/149‑1 (MQ=255)
    cTATCCAGAACGATGGTTACAAATCTATGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTgg                                                                                                                                                  >  3:519083/1‑149 (MQ=255)
     tATCCAGAACGATGGTTACAAATCTCTGCACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGc                                                                                                                                                 >  4:921703/1‑149 (MQ=255)
      aTCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCg                                                                                                                                                <  1:909184/149‑1 (MQ=255)
       tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCCTCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAATCATAGAATCTAAGATCCCTGCCATTTGGCgg                                                                                                                                               <  4:566884/149‑1 (MQ=255)
       tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgg                                                                                                                                               <  3:710035/149‑1 (MQ=255)
       tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgg                                                                                                                                               <  2:902373/149‑1 (MQ=255)
       tCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCgg                                                                                                                                               <  1:472414/149‑1 (MQ=255)
        ccAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCggg                                                                                                                                              <  2:485238/149‑1 (MQ=255)
          aGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGa                                                                                                                                            <  2:188565/149‑1 (MQ=255)
            aaCGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAtc                                                                                                                                          <  2:783604/149‑2 (MQ=255)
                 tGGTTACAAATCTATGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGAttttttt                                                                                                                                     <  4:519083/149‑1 (MQ=255)
                  ggTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTa                                                                                                                                    <  3:884855/149‑1 (MQ=255)
                  ggTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTa                                                                                                                                    >  4:500538/1‑149 (MQ=255)
                     tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                                                                                                                                 <  2:616848/149‑1 (MQ=255)
                     tACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTAttt                                                                                                                                 >  4:922262/1‑149 (MQ=255)
                      aCAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTg                                                                                                                                >  3:711491/1‑149 (MQ=255)
                        aaaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                                                                                                                              >  3:644274/1‑149 (MQ=255)
                        aaaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGtt                                                                                                                              >  1:738909/1‑149 (MQ=255)
                         aaTCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGttt                                                                                                                             <  4:746389/149‑1 (MQ=255)
                           tctcTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTc                                                                                                                           <  4:184321/149‑1 (MQ=255)
                            ctctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCa                                                                                                                          <  4:559286/149‑1 (MQ=255)
                              ctGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAgg                                                                                                                        <  3:96224/149‑1 (MQ=255)
                               tGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGa                                                                                                                       >  1:936531/1‑149 (MQ=255)
                                  aCGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaat                                                                                                                    <  2:247944/149‑1 (MQ=255)
                                    gAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataa                                                                                                                  <  2:421825/149‑1 (MQ=255)
                                     aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaa                                                                                                                 <  3:745319/149‑1 (MQ=255)
                                     aaGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGaaataaa                                                                                                                 <  4:305694/149‑1 (MQ=255)
                                              aaGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAatcgatc                                                                                                        >  1:146598/1‑149 (MQ=255)
                                                     cTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaat                                                                                                 >  3:128383/1‑149 (MQ=255)
                                                      ttCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGtaata                                                                                                >  3:131241/1‑149 (MQ=255)
                                                        cACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATaaa                                                                                              <  4:644274/149‑1 (MQ=255)
                                                         aCCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATaaaa                                                                                             <  3:249015/149‑1 (MQ=255)
                                                          ccATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAAt                                                                                            >  3:89181/1‑149 (MQ=255)
                                                           cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATc                                                                                           <  1:633212/149‑1 (MQ=255)
                                                           cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                           <  4:238828/101‑1 (MQ=255)
                                                           cATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCCGCCAAATGGCAGGGAt                                                                                                                                           >  3:238828/1‑101 (MQ=255)
                                                            aTCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATca                                                                                          >  1:827122/1‑148 (MQ=255)
                                                                        gcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAttttt                                                                              <  4:172672/149‑1 (MQ=255)
                                                                        gcTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAttttt                                                                              <  2:758062/149‑1 (MQ=255)
                                                                          tAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                            >  1:398026/1‑149 (MQ=255)
                                                                            aaGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTga                                                                          >  1:13144/1‑149 (MQ=255)
                                                                            aaGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTga                                                                          <  1:409060/149‑1 (MQ=255)
                                                                             aGGCCCGGCAGTTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTgaa                                                                         <  4:193801/149‑1 (MQ=255)
                                                                             aGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTgaa                                                                         <  4:810241/149‑1 (MQ=255)
                                                                              ggCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCta                                                                                         >  3:87377/1‑132 (MQ=255)
                                                                              ggCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCta                                                                                         <  4:87377/132‑1 (MQ=255)
                                                                                         tGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAAt                                                                     >  4:302599/1‑141 (MQ=255)
                                                                                         tGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAAt                                                                     <  3:302599/141‑1 (MQ=255)
                                                                                         tGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTgg                                                             <  1:113377/149‑1 (MQ=255)
                                                                                         tGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTgg                                                             >  3:657023/1‑149 (MQ=255)
                                                                                         tGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTgg                                                             <  3:582554/149‑1 (MQ=255)
                                                                                           gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAttttt                                                                              <  1:214050/130‑1 (MQ=255)
                                                                                           gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAttttt                                                                              >  2:214050/1‑130 (MQ=255)
                                                                                           gtaacgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGt                                                           >  1:659431/1‑149 (MQ=255)
                                                                                               cgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                       <  2:13144/149‑1 (MQ=255)
                                                                                               cgtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCaa                                                       <  1:400381/149‑1 (MQ=255)
                                                                                                gtaacCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAAt                                                      <  1:935120/149‑1 (MQ=255)
                                                                                                      aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAAt                                                <  3:167031/149‑1 (MQ=255)
                                                                                                      aGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAAt                                                <  2:797930/149‑1 (MQ=255)
                                                                                                          tGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATgcgc                                            >  2:486709/1‑149 (MQ=255)
                                                                                                                ctctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGcc                                      <  2:578393/149‑1 (MQ=255)
                                                                                                                ctctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGcc                                      <  1:749453/149‑1 (MQ=255)
                                                                                                                 tctGCTTAAAAGCACTGACTCTAAGATCCTTGCCATTTGGATGGTATTTGTTTATTTGTTTTCAGTAAATAAATAATCGATCGCGTAATAATATCTATTATTATTTTTGTGAAGATTAAATTTGGGTGCAATGAGAATGCGGAACGCCg                                     >  2:448613/1‑149 (MQ=255)
                                                                                                                 tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCg                                     >  3:766405/1‑149 (MQ=255)
                                                                                                                 tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCg                                     >  1:703733/1‑149 (MQ=255)
                                                                                                                 tctGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCg                                     >  3:201233/1‑149 (MQ=255)
                                                                                                                   tGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAttttt                                                                              >  1:6404/1‑106 (MQ=255)
                                                                                                                   tGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTAttttt                                                                              <  2:6404/106‑1 (MQ=255)
                                                                                                                    gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaa                                  <  4:865844/149‑1 (MQ=255)
                                                                                                                    gCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaa                                  <  3:500538/149‑1 (MQ=255)
                                                                                                                      ttAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCgtaagt                                >  3:649608/1‑149 (MQ=255)
                                                                                                                            gCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCg                          <  3:894972/149‑1 (MQ=255)
                                                                                                                            gCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCg                          <  2:659431/149‑1 (MQ=255)
                                                                                                                              acaGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCggg                        <  3:922262/149‑1 (MQ=255)
                                                                                                                               caGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGa                       >  4:545310/1‑149 (MQ=255)
                                                                                                                                aGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaa                      <  4:681942/149‑1 (MQ=255)
                                                                                                                                  aaTCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaata                    <  4:916567/149‑1 (MQ=255)
                                                                                                                                  aaTCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaata                    >  2:155812/1‑149 (MQ=255)
                                                                                                                                   aTCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaataa                   <  1:659594/149‑1 (MQ=255)
                                                                                                                                     cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                            <  4:699412/90‑1 (MQ=255)
                                                                                                                                     cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTtgt                                                                            >  3:699412/1‑90 (MQ=255)
                                                                                                                                     cTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAAtt                 <  2:936531/149‑1 (MQ=255)
                                                                                                                                       aaGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTc               >  2:30121/1‑149 (MQ=255)
                                                                                                                                         aaTCCCCGCCAAATGGCAGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTccc             >  4:746781/2‑149 (MQ=255)
                                                                                                                                          aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAgg                            >  2:785358/1‑133 (MQ=255)
                                                                                                                                          aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAgg                            <  1:785358/133‑1 (MQ=255)
                                                                                                                                          aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaataat                  <  1:35657/143‑1 (MQ=255)
                                                                                                                                          aTCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGaataat                  >  2:35657/1‑143 (MQ=255)
                                                                                                                                           tCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCccgc           >  1:240794/1‑149 (MQ=255)
                                                                                                                                              cTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGa        >  3:562178/1‑149 (MQ=255)
                                                                                                                                                 ccATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGa     <  1:57587/149‑1 (MQ=255)
                                                                                                                                                 ccACTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGaa                                                 >  2:525656/1‑105 (MQ=255)
                                                                                                                                                    tttGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGactc  <  1:155812/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCTGCTATCCAGAACGATGGTTACAAATCTCTGGACGAAGGTCAGAAAGTGTCCTTCACCATCGAAAGCGGCGCTAAAGGCCCGGCAGCTGGTAACGTAACCAGCCTGTAATCTCTGCTTAAAAGCACAGAATCTAAGATCCCTGCCATTTGGCGGGGATTTTTTTATTTGTTTTCAGGAAATAAATAATCGATCGCGTAATAAAATCTATTATTATTTTTGTGAAGAATAAATTTGGGTGCAATGAGAATGCGCAACGCCGTAAGTAAGGCGGGAATAATTTCCCGCCGAAGACTC  >  NC_000913/3720151‑3720447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: