Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,939,331 A→G 11.6% intergenic (+21/‑37) fucU → / → fucR L‑fucose mutarotase/l‑fucose operon activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,939,3310AG11.6% 264.1 / 23.4 103intergenic (+21/‑37)fucU/fucRL‑fucose mutarotase/l‑fucose operon activator
Reads supporting (aligned to +/- strand):  ref base A (41/50);  new base G (9/3);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 6.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-02

Only 100 of 103 total aligned reads displayed.

TCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGT  >  NC_000913/2939184‑2939475
                                                                                                                                                   |                                                                                                                                                
tCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATg                                                                                                                                                <  2:605572/150‑1 (MQ=255)
 ccTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGa                                                                                                                                               >  2:1348773/1‑150 (MQ=255)
 ccTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGa                                                                                                                                               >  2:230267/1‑150 (MQ=255)
 ccTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAAAGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCTCATGa                                                                                                                                               >  1:25373/1‑150 (MQ=255)
   tGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAcg                                                                                                                                             <  2:431411/150‑1 (MQ=255)
     aCATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAcgcg                                                                                                                                           <  1:1012067/150‑1 (MQ=255)
      catcatCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgg                                                                                                                                          >  1:728622/1‑150 (MQ=255)
          atcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgt                                                                                                                                      <  1:1270100/150‑2 (MQ=255)
          atcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcgg                                                                                                                                      <  1:1598210/150‑1 (MQ=255)
           tcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGt                                                                                                                                     <  1:1718742/150‑1 (MQ=255)
           tcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATAGTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGt                                                                                                                                     >  1:696060/1‑150 (MQ=255)
             cGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtacc                                                                                                                                   <  1:211059/150‑5 (MQ=255)
             cGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtacc                                                                                                                                   <  1:1000086/150‑5 (MQ=255)
                aTCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCg                                                                                                                                <  1:838887/150‑1 (MQ=255)
                aTCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCg                                                                                                                                <  1:752655/150‑1 (MQ=255)
                   aaTCGTTTTGCGTTTTATGAACGGGCGCAAAACGCATTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                             >  2:651638/1‑150 (MQ=255)
                   aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACTGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                             <  1:1500109/150‑1 (MQ=255)
                   aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                             >  1:1419904/1‑150 (MQ=255)
                   aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTGTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTGTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                             >  1:457127/1‑150 (MQ=255)
                       gTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgatg                                                                                                                         >  2:1068426/1‑150 (MQ=255)
                        ttttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATgg                                                                                                                        <  2:45289/150‑1 (MQ=255)
                        ttttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATgg                                                                                                                        <  1:624076/150‑1 (MQ=255)
                         tttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                       <  1:683739/150‑1 (MQ=255)
                         tttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                       <  2:728622/150‑1 (MQ=255)
                          ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATgggg                                                                                                                      <  1:1613328/150‑1 (MQ=255)
                                tttATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGaaaa                                                                                                                >  2:623507/1‑150 (MQ=255)
                                  tATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAat                                                                                                              <  2:1665255/150‑1 (MQ=255)
                                  tATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAat                                                                                                              >  1:437529/1‑150 (MQ=255)
                                   ttGAACGGGCGCAAAAAGCCTTTGCGATAGTTATCACCGGGGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                             <  2:834743/149‑1 (MQ=255)
                                   aTGAACGGGCGCAAAAAGGCTTTGCGAGCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                             <  2:163626/150‑1 (MQ=255)
                                    ttaaCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAatat                                                                                                            <  1:599015/148‑1 (MQ=255)
                                    tGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTTAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAatat                                                                                                            <  1:1410862/150‑1 (MQ=255)
                                       aCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaa                                                                                                         <  2:683459/150‑1 (MQ=255)
                                         gtgcgcAAAAAGCCTTTGCGATCGTTATCACAGGCCAACGAGCTAAGTACGGGAATATTCTTTTAAAAAAATGGGTAACACAGTAATCTCATACCGGTACGCCCGCATGCCGAGGGCGGTTATCGAATGATGGGGTGAAAAATATCaaac                                                                                                       >  2:770113/3‑149 (MQ=255)
                                          gacgcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGc                                                                                                      >  1:1193024/3‑150 (MQ=255)
                                           gcgcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCg                                                                                                     <  2:812959/150‑1 (MQ=255)
                                            cgcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCTTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCgg                                                                                                    <  1:722920/150‑1 (MQ=255)
                                            cgcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGTTGGGGTGAAAAATATGAAAGCgg                                                                                                    <  1:230267/150‑1 (MQ=255)
                                            cgcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCgg                                                                                                    <  1:915665/150‑1 (MQ=255)
                                             gcAACAAGCCTTTGCGATCGTTATCACAGGCGAACGAGAGAAGTACTGGTATATTCTTTTAAAAAAAGGGGCCACCCCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGc                                                                                                   <  2:725584/150‑1 (MQ=255)
                                                    gCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAg                                                                                            <  1:1080737/150‑1 (MQ=255)
                                                       tttGCGATCGTTATCACAGTCGATCGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTTACACAGTAATATCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAATAATATGAAAGCGGCACGCCAGCaa                                                                                         >  2:708222/1‑150 (MQ=255)
                                                        ttGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAg                                                                                        <  2:1391916/150‑1 (MQ=255)
                                                        ttGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGAAACACCGTAATCTCATCCCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGGGAAAAATATGAAAGCGGCACGCCAGCAAg                                                                                        <  2:1656744/150‑1 (MQ=255)
                                                         tGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGc                                                                                       <  2:1339806/150‑1 (MQ=255)
                                                           cGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGa                                                                                     >  1:273593/1‑150 (MQ=255)
                                                           cGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGAGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGa                                                                                     >  2:115288/1‑150 (MQ=255)
                                                                gTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCg                                                                               >  2:335871/1‑150 (MQ=255)
                                                                  tATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGa                                                                              <  2:701774/150‑1 (MQ=255)
                                                                   aTCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAc                                                                             >  2:828584/1‑150 (MQ=255)
                                                                    tCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTCACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAcc                                                                            <  2:1320979/150‑1 (MQ=255)
                                                                    tCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAcc                                                                            <  1:1348773/150‑1 (MQ=255)
                                                                    tCACAGGCGAACGAGCGAAGTAAGGGAATATTCTTTTAAAAAAAGGGGTAAAACCGTAATCTAATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAcc                                                                            <  2:122111/150‑1 (MQ=255)
                                                                     cacaGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGAAACACAGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTAATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACct                                                                           >  2:679173/1‑150 (MQ=255)
                                                                      acaGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCCCGCCAGCAAGCGATAGTCGACctg                                                                          >  1:1290485/1‑150 (MQ=255)
                                                                         ggCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctg                                                                       >  1:622926/1‑150 (MQ=255)
                                                                          gCGAACGAGCGAAGGACGGGAATATTCTTTTAAAAAAAGGGGGAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgc                                                                      <  1:830128/150‑1 (MQ=255)
                                                                           cGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgct                                                                     <  1:1091333/150‑1 (MQ=255)
                                                                             aaCGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGa                                                                   >  2:1660899/1‑150 (MQ=255)
                                                                                    gAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAc                                                            >  1:115428/1‑150 (MQ=255)
                                                                                         aCGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGcc                                                       <  2:140156/150‑1 (MQ=255)
                                                                                         aCGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGcc                                                       <  2:488320/150‑1 (MQ=255)
                                                                                           gggAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTg                                                     <  1:473495/150‑1 (MQ=255)
                                                                                            ggAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGa                                                    >  1:762409/1‑150 (MQ=255)
                                                                                            ggAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGa                                                    >  2:806537/1‑150 (MQ=255)
                                                                                            ggAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGa                                                    >  1:1194924/1‑150 (MQ=255)
                                                                                                tatTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCAc                                                >  1:628829/1‑150 (MQ=255)
                                                                                                  ttCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTCCTGCTGAACCATACCAGCCTGACCACgg                                              >  2:1034729/1‑150 (MQ=255)
                                                                                                    cTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                            >  1:706194/1‑150 (MQ=255)
                                                                                                    cTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                            >  1:1654769/1‑150 (MQ=255)
                                                                                                    cTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTAAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                            >  1:1124056/1‑150 (MQ=255)
                                                                                                     ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAg                                           <  1:1050435/150‑1 (MQ=255)
                                                                                                     ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAg                                           >  1:1573418/1‑150 (MQ=255)
                                                                                                     ttttAAAAAAAGGGGTAACACCGTAATCTCAAACCGGTACGCCCGGATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCAGCTGAACCATACCAGCCTGCCCACGGAAg                                           >  2:56092/1‑150 (MQ=255)
                                                                                                          aaaaaaGGGGTAACACCGTAATCTCATACCGTTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACATGCTGCTGAACAATACCAGCCTGACCACGGAAGctctc                                      >  1:1168143/1‑150 (MQ=255)
                                                                                                               aGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGa                                 <  2:1486143/150‑1 (MQ=255)
                                                                                                                ggggTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGaa                                <  2:1478918/150‑1 (MQ=255)
                                                                                                                ggggTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGaa                                >  2:622656/1‑150 (MQ=255)
                                                                                                                   gTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGACTGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAg                             <  1:889138/150‑1 (MQ=255)
                                                                                                                       cacCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCCGCCTGACCACGGAAGCTCTCTCTGAACAGCTaa                         >  2:153409/1‑150 (MQ=255)
                                                                                                                        acCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaaa                        >  1:1169403/1‑150 (MQ=255)
                                                                                                                         ccGTAATCTCATACCTGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATAGGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACAACGGAAGCTCTCTCTGAACAGCTAAAg                       <  2:1031474/150‑1 (MQ=255)
                                                                                                                          cGTAATCTCATACCTGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAgg                      <  2:25373/150‑1 (MQ=255)
                                                                                                                            tAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTc                    <  1:1310201/150‑1 (MQ=255)
                                                                                                                              aTCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAg                  <  2:273593/150‑1 (MQ=255)
                                                                                                                                  cATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAACAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGCACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTaaa              >  1:71959/1‑150 (MQ=255)
                                                                                                                                  cATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTaaa              >  1:627326/1‑150 (MQ=255)
                                                                                                                                    tACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGa            <  1:623507/150‑1 (MQ=255)
                                                                                                                                     cgataaccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGaa           >  2:725072/8‑150 (MQ=255)
                                                                                                                                      gataaccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGCCCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGaaa          >  2:813208/7‑150 (MQ=255)
                                                                                                                                      gataaccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGaaa          >  2:539961/7‑150 (MQ=255)
                                                                                                                                      ccGGTACGCCAGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGCCCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCGCTCTCTGCACAGCTAAAGGTCAGTAAAGaaa          >  1:1593894/1‑150 (MQ=255)
                                                                                                                                        taaccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCCCGGAAGCTCTCGCTGAACAGCTAAAGGTCAGTAAAGAAAcc        >  1:1283267/5‑150 (MQ=255)
                                                                                                                                         aaccGCCCGCGTCATGCGGGCGGTTATCGAATGGTGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCATGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCa       >  2:831669/4‑150 (MQ=255)
                                                                                                                                         aaccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCa       >  1:1492357/4‑150 (MQ=255)
                                                                                                                                          accGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCCGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt      >  2:1130270/3‑150 (MQ=255)
                                                                                                                                          accGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt      >  1:971159/3‑150 (MQ=255)
                                                                                                                                          accGCCCGCGTCATGCGGGCGGTTATCGAATCATGGGGTGAAAAATATGAAAGCGGCACGCCCGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAACCCAt      >  2:243812/3‑150 (MQ=255)
                                                                                                                                             gCCCGCATGACGTGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAGGTCAGTAAAGAAACCATtcgt  >  2:789234/1‑150 (MQ=255)
                                                                                                                                             gCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcg   <  1:1245594/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGT  >  NC_000913/2939184‑2939475

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: