Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 964,148 T→C 6.0% intergenic (+36/‑172) ihfB → / → ycaI integration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913964,1480TC 6.0% 354.9 / 9.3 117intergenic (+36/‑172)ihfB/ycaIintegration host factor (IHF), DNA‑binding protein, beta subunit/ComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (50/60);  new base C (4/3);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

Only 100 of 117 total aligned reads displayed.

TACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACA  >  NC_000913/964004‑964295
                                                                                                                                                |                                                                                                                                                   
tACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCgg                                                                                                                                                >  2:839793/1‑150 (MQ=255)
 aCCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTTAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTAAAACTTGAACGAGAGAAAAGCACCTGTCggg                                                                                                                                               >  2:1184548/1‑150 (MQ=255)
   cGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTg                                                                                                                                             <  1:1798652/150‑1 (MQ=255)
     gACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCt                                                                                                                                           >  1:1530258/1‑150 (MQ=255)
         tAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCttttt                                                                                                                                       >  2:1633705/1‑150 (MQ=255)
          aaTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCtttttt                                                                                                                                      <  2:495452/150‑1 (MQ=255)
            tCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGAACCTGTCGGTTGCTTTTTTCa                                                                                                                                    >  1:742614/1‑150 (MQ=255)
             ccGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAt                                                                                                                                   >  1:1820568/1‑150 (MQ=255)
              cGAAGACTGGCGATAAAGTAGAACTGGAATGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAtt                                                                                                                                  <  1:210481/150‑1 (MQ=255)
              cGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAtt                                                                                                                                  >  1:127643/1‑150 (MQ=255)
              cGAAGACTGGCGATAAAGTAGAACTGGAACGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAtt                                                                                                                                  <  1:681652/150‑1 (MQ=255)
               gAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAttt                                                                                                                                 <  1:1218376/150‑1 (MQ=255)
                aaGACTGGCGATATAGTATAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAATAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGCGAAAAGCACCCGACAGGTGCTTTtctctcgtt                                                                                                                                <  1:1332831/150‑6 (MQ=255)
                aaGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtc                                                                                                                                <  2:700392/150‑1 (MQ=255)
                 aGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCCGACAGGTGCTTTtctctcgttc                                                                                                                               <  1:959085/150‑7 (MQ=255)
                  gACTGGCGATAACGTACAACTGGAAGTAAAATAAGTTCATCACATTAAACCTGGTAAATAACTTCGAGATCTCTCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGTGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                                                                              <  2:227260/150‑1 (MQ=255)
                  gACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTtctc                                                                                                                              <  2:1473915/150‑1 (MQ=255)
                      ggCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAAt                                                                                                                          >  2:509636/1‑150 (MQ=255)
                         gATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctg                                                                                                                       >  2:1701882/1‑150 (MQ=255)
                             aaGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaac                                                                                                                   <  2:550662/150‑1 (MQ=255)
                               gTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactg                                                                                                                 <  2:850736/150‑1 (MQ=255)
                                  gAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATctggaactggaa                                                                                                              >  2:434714/1‑150 (MQ=255)
                                     cTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCt                                                                                                           >  1:333031/1‑150 (MQ=255)
                                      tGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTg                                                                                                          <  2:88725/150‑1 (MQ=255)
                                         aaggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCt                                                                                                       >  1:737416/1‑150 (MQ=255)
                                          aggaaAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTc                                                                                                      >  2:1164942/1‑150 (MQ=255)
                                             aaaaTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCa                                                                                                   <  2:474371/150‑1 (MQ=255)
                                              aaaTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCag                                                                                                  >  1:1803169/1‑150 (MQ=255)
                                               aaTACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCaga                                                                                                 >  1:607454/1‑150 (MQ=255)
                                                   cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                                                                             <  1:1044716/150‑1 (MQ=255)
                                                   cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                                                                             <  2:228495/150‑1 (MQ=255)
                                                   cGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCAAATATTTACGGTTATGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGttt                                                                                             <  2:787543/150‑1 (MQ=255)
                                                     ttCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGACTCGAAGAGttttt                                                                                           >  2:842529/1‑149 (MQ=255)
                                                      tCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGa                                                                                          <  2:1104877/150‑1 (MQ=255)
                                                        cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                        >  1:806054/1‑150 (MQ=255)
                                                        cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAAc                                                                                        >  2:225663/1‑150 (MQ=255)
                                                        cTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGATTCCTCGCAGAGTATTGAAc                                                                                        >  2:1219475/1‑150 (MQ=255)
                                                            cTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTATCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGtt                                                                                    >  1:1490419/1‑150 (MQ=255)
                                                            cTTTAAACCTGGTAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGttt                                                                                   >  2:5759/1‑150 (MQ=255)
                                                             tttAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAAAGAGAGAAAAGCACCTTTCGGGTGCTTTTTTCATTTCTATAATCTTGAACTGGAAGCTGCATCGCAGAGGTTTGAACAGttt                                                                                   >  1:489495/1‑150 (MQ=255)
                                                               tAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTc                                                                                 <  1:138330/150‑1 (MQ=255)
                                                                 aaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAc                                                                               <  1:702775/150‑1 (MQ=255)
                                                                 aaCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTATCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAc                                                                               >  2:1131205/1‑150 (MQ=255)
                                                                  aCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCAcc                                                                              <  2:153703/150‑1 (MQ=255)
                                                                    cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCt                                                                            >  2:22283/1‑150 (MQ=255)
                                                                    cTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACACt                                                                            >  1:935736/1‑150 (MQ=255)
                                                                    cTGGTAAAGAACTGCGCGATCGCGCCAATATTTAAGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAAATGGAAGCTGCCACGCACAGTTTTGAACAGTTTTCACCCt                                                                            >  1:1750202/1‑150 (MQ=255)
                                                                        tAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCg                                                                        >  2:550195/1‑150 (MQ=255)
                                                                        tAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCg                                                                        <  1:181756/150‑1 (MQ=255)
                                                                         aaaGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGt                                                                       >  2:820843/1‑150 (MQ=255)
                                                                           aGACCTGCGCGATCGCGCCAATATTTTCGGTTCAGTTTTTTAGTCTAACTTGAACGAGATAAAAGCACGTGTCGGGTGCTTTTTTCATTTCTCTTATCTGTAACTTGAAGCTGCCTCGCAGCGTTTTGAACAGTTTTCACCCTTTCGTTa                                                                     <  1:947816/150‑1 (MQ=255)
                                                                             aaCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTaaa                                                                   >  2:1008125/1‑150 (MQ=255)
                                                                             aaCTGCGCGATCGCGCCAATATTTACGGGTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGGCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTaaa                                                                   <  2:1663388/150‑1 (MQ=255)
                                                                              aCTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAATCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGTACAGTTTTCACCCTTTCGTTAAAt                                                                  <  1:587538/150‑1 (MQ=255)
                                                                                 gcgcGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttct                                                               <  1:1632721/150‑1 (MQ=255)
                                                                                 gcgcGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGGCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttct                                                               <  1:977930/150‑1 (MQ=255)
                                                                                  cgcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTAGCAGAGTTTTGAACAGTTTTCACCCTTTCGTAAAAttctt                                                              >  1:1439705/1‑150 (MQ=255)
                                                                                   gcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAGCGAGAGAAAAGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc                                                             <  1:1710810/150‑1 (MQ=255)
                                                                                   gcgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAAttcttc                                                             <  1:1147498/150‑1 (MQ=255)
                                                                                    cgATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATtcttct                                                            >  1:1777411/1‑150 (MQ=255)
                                                                                     gATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTg                                                           <  2:1069306/150‑1 (MQ=255)
                                                                                         gcgcCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAATCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAata                                                       <  1:1219475/150‑1 (MQ=255)
                                                                                          cgcCAATATGTACGGGTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGTACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAatat                                                      <  1:896625/150‑1 (MQ=255)
                                                                                            ccAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGc                                                    >  1:1794050/1‑150 (MQ=255)
                                                                                            ccAATATTTACGGTTAAGTTTTATACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGc                                                    >  1:428187/1‑150 (MQ=255)
                                                                                               atatTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTc                                                 <  1:1069107/150‑1 (MQ=255)
                                                                                                 atTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCgg                                               <  1:468760/150‑1 (MQ=255)
                                                                                                       ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACg                                         <  2:607454/150‑1 (MQ=255)
                                                                                                       ggTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACg                                         <  1:37893/150‑1 (MQ=255)
                                                                                                        gTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCCCCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGc                                        >  2:10032/1‑150 (MQ=255)
                                                                                                          tAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCaa                                      >  1:1732717/1‑150 (MQ=255)
                                                                                                                ttttACTCAAACTTGAACGAGGGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTcc                                <  1:42918/150‑1 (MQ=255)
                                                                                                                 tttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTccc                               <  2:215877/150‑1 (MQ=255)
                                                                                                                  ttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCcca                              <  1:16337/150‑1 (MQ=255)
                                                                                                                  ttACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCcca                              <  2:333031/150‑1 (MQ=255)
                                                                                                                   tACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccac                             <  2:1664589/150‑1 (MQ=255)
                                                                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacc                            <  1:230747/150‑1 (MQ=255)
                                                                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacc                            <  2:802512/150‑1 (MQ=255)
                                                                                                                    aCTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacc                            <  2:127643/150‑1 (MQ=255)
                                                                                                                     cTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCccacca                           <  2:798988/150‑1 (MQ=255)
                                                                                                                       cAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCCCCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAgc                         >  1:1398186/1‑149 (MQ=255)
                                                                                                                       cAAACTTGAACGAGAGAAAAGCACCTGTCGGGAGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTAGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGa                         <  1:512979/150‑1 (MQ=255)
                                                                                                                        aaaCTTGAACGAGAGAATAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACAGTAAGCTTCCTCTCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGAc                        <  1:1703518/150‑1 (MQ=255)
                                                                                                                          aCTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCAGTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCCCCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCATAATTACCACCAGACaa                      >  2:1288738/1‑150 (MQ=255)
                                                                                                                           cTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAATGCAAAATTCCCACCAGACAAc                     >  1:527178/1‑150 (MQ=255)
                                                                                                                           cTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAAc                     >  1:1240466/1‑150 (MQ=255)
                                                                                                                                aCGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTc                <  2:1530258/150‑1 (MQ=255)
                                                                                                                                 cGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTca               <  1:1499987/150‑1 (MQ=255)
                                                                                                                                   tgaaaaaaGCACCCGACAGGTGCTTTTTTCATGTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTGTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATACCCACCAGACAACCGCTcaac             >  1:1586758/5‑150 (MQ=255)
                                                                                                                                   agagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaac             >  2:794801/1‑150 (MQ=255)
                                                                                                                                   agagaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaac             <  2:806054/150‑1 (MQ=255)
                                                                                                                                    gaaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaaca            >  1:737535/4‑150 (MQ=255)
                                                                                                                                     aaaaaaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTcaacaa           >  2:41897/3‑150 (MQ=255)
                                                                                                                                      gaAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACaaa          >  2:556085/1‑150 (MQ=255)
                                                                                                                                        aaaGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGt        >  2:546698/1‑150 (MQ=255)
                                                                                                                                         aaGCACCCGACAGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGAAAAATTCCCACCAGACAAACGCTCAACAAAGtt       >  2:1131077/1‑150 (MQ=255)
                                                                                                                                           gCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGc     <  2:1820568/150‑1 (MQ=255)
                                                                                                                                            cACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTTGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGca    <  1:1735843/150‑1 (MQ=255)
                                                                                                                                             aCCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGACTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGcac   <  1:1164942/150‑1 (MQ=255)
                                                                                                                                              ccTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGcaca  >  1:942425/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCGGGTGCTTTTTTCATTTCTCTAATCTGGAACTGGAAGCTGCCTCGCAGAGTTTTGAACAGTTTTCACCCTTTCGTTAAATTCTTCTGAATATGCCTCGGGGAACGCAAAATTCCCACCAGACAACCGCTCAACAAAGTTGCACA  >  NC_000913/964004‑964295

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: