Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,939,333 G→C 5.2% 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,3330GC 5.2% 548.9 / 17.1 155intergenic (+23/‑35)fucU/fucRL‑fucose mutarotase/l‑fucose operon activator
Reads supporting (aligned to +/- strand):  ref base G (68/78);  new base C (5/3);  total (74/81)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

Only 100 of 156 total aligned reads displayed.

CTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAT  >  NC_000913/2939186‑2939481
                                                                                                                                                   |                                                                                                                                                    
cTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGa                                                                                                                                                     >  3:248996/1‑149 (MQ=255)
 tGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                                                                    >  1:1175416/1‑149 (MQ=255)
 tGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                                                                    >  2:1344229/1‑149 (MQ=255)
 tGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                                                                    <  3:563214/149‑1 (MQ=255)
 tGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                                                                    >  2:24292/1‑149 (MQ=255)
 tGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                                                                    >  1:534226/1‑149 (MQ=255)
        atcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcg                                                                                                                                             <  3:1059550/149‑1 (MQ=255)
          ccGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgta                                                                                                                                           <  1:890258/149‑3 (MQ=255)
              aTCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATc                                                                                                                                       >  2:165231/1‑149 (MQ=255)
              aTCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATc                                                                                                                                       <  1:1344229/149‑1 (MQ=255)
               tCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCg                                                                                                                                      >  1:1456894/1‑149 (MQ=255)
                cAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGa                                                                                                                                     <  4:1282244/149‑1 (MQ=255)
                 aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGaa                                                                                                                                    >  4:54400/1‑149 (MQ=255)
                  aTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                                   >  3:384037/1‑149 (MQ=255)
                  aTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                                   >  3:320978/1‑149 (MQ=255)
                      ttttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgatg                                                                                                                               <  1:408206/149‑1 (MQ=255)
                      ttttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgatg                                                                                                                               <  2:384173/149‑1 (MQ=255)
                        ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGCAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                             <  1:467683/149‑1 (MQ=255)
                        ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                             <  2:593473/149‑1 (MQ=255)
                        ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                             <  2:794486/149‑1 (MQ=255)
                        ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                             <  2:1159381/149‑1 (MQ=255)
                          gCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGt                                                                                                                           >  3:115454/1‑149 (MQ=255)
                            gTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGa                                                                                                                         <  3:967030/149‑1 (MQ=255)
                            gTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGa                                                                                                                         >  4:1323595/1‑149 (MQ=255)
                                 aTGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAat                                                                                                                    <  3:95437/149‑1 (MQ=255)
                                  tGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                   <  1:1106387/149‑1 (MQ=255)
                                  tGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                   <  4:854495/149‑1 (MQ=255)
                                   gAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAGAatat                                                                                                                  >  2:616032/1‑149 (MQ=255)
                                   gAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAatat                                                                                                                  <  1:617458/149‑1 (MQ=255)
                                    aaCGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATg                                                                                                                 >  4:1202850/1‑149 (MQ=255)
                                       gggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                                                                              <  4:248996/149‑1 (MQ=255)
                                       gggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                                                                              >  2:240321/1‑149 (MQ=255)
                                         gcgcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGc                                                                                                            >  4:463779/1‑149 (MQ=255)
                                               aaaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACg                                                                                                      >  1:1026189/1‑149 (MQ=255)
                                                   ccTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAg                                                                                                  >  1:956312/1‑149 (MQ=255)
                                                          gATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGa                                                                                           <  2:956312/149‑1 (MQ=255)
                                                            tCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                         <  1:1376722/149‑1 (MQ=255)
                                                            tCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                         <  4:11864/149‑1 (MQ=255)
                                                            tCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                         <  1:465419/149‑1 (MQ=255)
                                                            tCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                         >  2:204543/1‑149 (MQ=255)
                                                            tCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                         <  4:1244230/149‑1 (MQ=255)
                                                            tCGGTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACACGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGAGGCACGCCAGCAAGCGATa                                                                                         >  2:1234994/1‑149 (MQ=255)
                                                                tATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCg                                                                                     <  4:814217/149‑1 (MQ=255)
                                                                  tCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAc                                                                                   <  2:883489/149‑1 (MQ=255)
                                                                   cacaGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAcc                                                                                  >  3:766438/1‑149 (MQ=255)
                                                                    acaGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACct                                                                                 >  2:310645/1‑149 (MQ=255)
                                                                    acaGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACct                                                                                 >  2:1169375/1‑149 (MQ=255)
                                                                    acaGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACct                                                                                 >  2:506988/1‑149 (MQ=255)
                                                                       ggCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAc                                                                                   <  3:933186/144‑1 (MQ=255)
                                                                       ggCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgct                                                                              <  4:678485/149‑1 (MQ=255)
                                                                             cGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGaa                                                                        <  1:1166553/149‑1 (MQ=255)
                                                                              gAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAAc                                                                       <  3:463779/149‑1 (MQ=255)
                                                                                   aaGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAc                                                                  <  3:1232586/149‑1 (MQ=255)
                                                                                    aGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCg                                                                                                           >  4:793520/1‑107 (MQ=255)
                                                                                    aGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCg                                                                                                           <  4:1348521/107‑1 (MQ=255)
                                                                                        cGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGcc                                                             >  3:820719/1‑149 (MQ=255)
                                                                                                ttCTTTTAAAAAAGGGGTGAACACCTTAATCTCATCCCGGTACGCCCCCATTACGCGGGCGGTTATCGAAATTTGGGTCGAAAAATATGAAATCGGCACGCCAGCAAGCGATAGTCGACCTGCTGTTCAACCATACAAGCCTGCCCACg                                                     >  2:740065/1‑149 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                  >  3:665273/1‑149 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                  >  4:1284299/1‑149 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                  <  1:1234994/149‑1 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                  >  4:374535/1‑149 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                  <  3:187324/149‑1 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCAACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                  >  4:950788/1‑149 (MQ=255)
                                                                                                    tttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAg                                                                                                             >  3:614437/1‑89 (MQ=255)
                                                                                                    tttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAg                                                                                                             <  4:614437/89‑1 (MQ=255)
                                                                                                       aaaaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctc                                              >  3:559574/1‑149 (MQ=255)
                                                                                                         aaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGc                                                                                                            <  2:399949/85‑1 (MQ=255)
                                                                                                         aaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGc                                                                                                            >  1:399949/1‑85 (MQ=255)
                                                                                                          aaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctct                                           <  1:1456451/149‑1 (MQ=255)
                                                                                                          aaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctct                                           <  4:766438/149‑1 (MQ=255)
                                                                                                          aaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctct                                           <  3:776951/149‑1 (MQ=255)
                                                                                                           aaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCa                                                                     <  1:859949/122‑1 (MQ=255)
                                                                                                           aaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCa                                                                     >  2:859949/1‑122 (MQ=255)
                                                                                                           aaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATa                                                                   >  1:1311115/1‑124 (MQ=255)
                                                                                                           aaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctctc                                          >  2:743857/1‑149 (MQ=255)
                                                                                                                  tAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAg                                   <  2:1092438/149‑1 (MQ=255)
                                                                                                                   aaCACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGc                                  >  3:401063/1‑149 (MQ=255)
                                                                                                                   aaCACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGc                                  <  2:1059085/149‑1 (MQ=255)
                                                                                                                      acCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATCCCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaa                               <  1:851054/149‑1 (MQ=255)
                                                                                                                      acCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaa                               >  2:314163/1‑149 (MQ=255)
                                                                                                                       ccGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaaa                              <  2:572922/149‑1 (MQ=255)
                                                                                                                        cGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctctc                                          >  2:648884/1‑136 (MQ=255)
                                                                                                                        cGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctctc                                          <  1:648884/136‑1 (MQ=255)
                                                                                                                          tAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGt                           <  1:165231/149‑1 (MQ=255)
                                                                                                                          tAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGt                           >  2:1017537/1‑149 (MQ=255)
                                                                                                                            aTCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCa                         <  1:1293449/149‑1 (MQ=255)
                                                                                                                               tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGACGct                                               <  3:1135858/124‑1 (MQ=255)
                                                                                                                               tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGct                                               >  4:1135858/1‑124 (MQ=255)
                                                                                                                               tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTa                      >  1:786794/1‑149 (MQ=255)
                                                                                                                                cATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTaa                     >  2:956439/1‑149 (MQ=255)
                                                                                                                                     cGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGaaa                <  3:931681/149‑1 (MQ=255)
                                                                                                                                      taaccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAAc               >  1:184144/5‑149 (MQ=255)
                                                                                                                                      ggTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAAc               >  3:1207091/1‑149 (MQ=255)
                                                                                                                                         ccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt            >  4:1086931/2‑149 (MQ=255)
                                                                                                                                         ccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt            >  2:964237/2‑149 (MQ=255)
                                                                                                                                         aCGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt            >  1:774143/1‑149 (MQ=255)
                                                                                                                                             ccGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCATGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcgt        >  4:892553/1‑149 (MQ=255)
                                                                                                                                              cGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCCGAACAGCTAAAGGTCAGTAAAGAAACCATtcgtc       <  3:192694/149‑1 (MQ=255)
                                                                                                                                                cATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTcgc     <  3:1284299/149‑1 (MQ=255)
                                                                                                                                                   gACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAt  >  2:1277595/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAT  >  NC_000913/2939186‑2939481

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