Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,939,333 G→C 11.4% intergenic (+23/‑35) 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,3330GC11.4% 299.8 / 25.0 104intergenic (+23/‑35)fucU/fucRL‑fucose mutarotase/l‑fucose operon activator
Reads supporting (aligned to +/- strand):  ref base G (42/50);  new base C (9/3);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 6.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

Only 100 of 105 total aligned reads displayed.

TCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAT  >  NC_000913/2939184‑2939481
                                                                                                                                                     |                                                                                                                                                    
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)
                   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)
                                                  aaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGcc                                                                                                    <  2:1545855/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)
                                                           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:488320/150‑1 (MQ=255)
                                                                                         aCGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGcc                                                             <  2:140156/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)
                                                                                                  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)
                                                                                                      tttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGc                                                <  1:47912/150‑1 (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)
                                                                                                                                      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)
                                                                                                                                                    tGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAt  >  1:563115/1‑150 (MQ=255)
                                                                                                                                                    tGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAt  <  2:437529/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAT  >  NC_000913/2939184‑2939481

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: