Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,939,335 C→T 6.3% intergenic (+25/‑33) 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,3350CT 6.3% 252.2 / 6.0 79intergenic (+25/‑33)fucU/fucRL‑fucose mutarotase/l‑fucose operon activator
Reads supporting (aligned to +/- strand):  ref base C (34/40);  new base T (2/3);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

ACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATC  >  NC_000913/2939189‑2939482
                                                                                                                                                   |                                                                                                                                                   
aCATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAcgc                                                                                                                                                   >  1:223853/1‑149 (MQ=255)
  atcatcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgg                                                                                                                                                 >  2:19754/1‑149 (MQ=255)
      tcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcgg                                                                                                                                             >  3:223829/1‑149 (MQ=255)
        cGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtac                                                                                                                                           <  3:343113/149‑4 (MQ=255)
        cGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgta                                                                                                                                            <  3:686337/148‑3 (MQ=255)
        cGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgta                                                                                                                                            >  4:686337/1‑146 (MQ=255)
           aTCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATc                                                                                                                                        >  1:629865/1‑149 (MQ=255)
           aTCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATc                                                                                                                                        <  1:548876/149‑1 (MQ=255)
            tCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCg                                                                                                                                       >  3:312567/1‑149 (MQ=255)
              aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGaa                                                                                                                                     >  4:408751/1‑149 (MQ=255)
              aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGaa                                                                                                                                     <  4:468196/149‑1 (MQ=255)
              aaTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCAAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGaa                                                                                                                                     >  3:121977/1‑149 (MQ=255)
               aTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                                                                    <  3:541071/149‑1 (MQ=255)
                     ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                              <  3:595978/149‑1 (MQ=255)
                         gTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGa                                                                                                                          <  1:662050/149‑1 (MQ=255)
                          ttttATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGaa                                                                                                                         >  2:143094/1‑149 (MQ=255)
                            ttATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGaaaa                                                                                                                       >  4:366037/1‑149 (MQ=255)
                               tGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTGTCGAATGATGGGGTGAAAAata                                                                                                                    <  4:20024/149‑1 (MQ=255)
                               tGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    <  3:370654/149‑1 (MQ=255)
                               tGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    <  2:778171/149‑1 (MQ=255)
                                   cGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaa                                                                                                                <  3:114488/149‑1 (MQ=255)
                                    gggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaag                                                                                                               >  2:90200/1‑148 (MQ=255)
                                    gggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                                                                               <  4:104328/149‑1 (MQ=255)
                                           aaaaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCAc                                                                                                        <  4:710843/149‑1 (MQ=255)
                                            aaaGCCTTTGCGATCGTTAGCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACg                                                                                                       >  3:517159/1‑149 (MQ=255)
                                             aaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGc                                                                                                      <  1:388808/149‑1 (MQ=255)
                                               gCCTTTGCGATCGTTATCACAGTCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTATCACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCa                                                                                                    <  2:538623/149‑1 (MQ=255)
                                               gCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCa                                                                                                    <  4:144892/149‑1 (MQ=255)
                                               gCCTTTGCGATCGTTAGCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCa                                                                                                    <  1:658122/149‑1 (MQ=255)
                                                  tttGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCCTTTAAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGc                                                                                                  >  4:320551/1‑149 (MQ=255)
                                                         tCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGaa                                                                                                                                     >  2:336802/1‑106 (MQ=255)
                                                         tCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGaa                                                                                                                                     <  1:336802/106‑1 (MQ=255)
                                                           gTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGt                                                                                        >  2:561819/1‑149 (MQ=255)
                                                            ttATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTc                                                                                       <  1:877461/149‑1 (MQ=255)
                                                            ttATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTc                                                                                       <  1:798658/149‑1 (MQ=255)
                                                                cacaGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAcc                                                                                   <  1:143094/149‑1 (MQ=255)
                                                                   aGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgc                                                                                <  1:19754/149‑1 (MQ=255)
                                                                     gCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctg                                                                              >  3:147515/1‑149 (MQ=255)
                                                                          cGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGaa                                                                         <  4:312567/149‑1 (MQ=255)
                                                                           gAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAAc                                                                        <  2:629865/149‑1 (MQ=255)
                                                                           gAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAAc                                                                        <  2:602315/149‑1 (MQ=255)
                                                                             gCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    >  4:261772/1‑103 (MQ=255)
                                                                             gCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    <  3:261772/103‑1 (MQ=255)
                                                                             gCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    >  4:325790/1‑103 (MQ=255)
                                                                             gCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    <  3:325790/103‑1 (MQ=255)
                                                                                   tACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAg                                                                >  2:710590/1‑149 (MQ=255)
                                                                                      gggAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCt                                                             <  1:378187/149‑1 (MQ=255)
                                                                                       ggAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTg                                                            >  2:541514/1‑149 (MQ=255)
                                                                                        gAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                          >  4:842553/1‑118 (MQ=255)
                                                                                        gAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                          <  3:842553/118‑1 (MQ=255)
                                                                                        gAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGa                                                           >  2:633981/1‑149 (MQ=255)
                                                                                        gAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGa                                                           >  3:634638/1‑149 (MQ=255)
                                                                                                 ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                   >  2:583943/1‑149 (MQ=255)
                                                                                                 ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                   <  3:675510/149‑1 (MQ=255)
                                                                                                 ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                   >  3:387113/1‑149 (MQ=255)
                                                                                                 ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                   <  3:709747/149‑1 (MQ=255)
                                                                                                 ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGa                                                    >  4:528353/1‑148 (MQ=255)
                                                                                                 ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGa                                                    <  3:528353/148‑1 (MQ=255)
                                                                                                 ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCACGCATGACGCGGGCGGTTATCGAAGGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGaa                                                   >  4:430382/1‑149 (MQ=255)
                                                                                                   ttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAg                                                                                                   >  2:448501/1‑99 (MQ=255)
                                                                                                   ttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAg                                                                                                   <  1:448501/99‑1 (MQ=255)
                                                                                                        aaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctct                                            >  1:347554/1‑149 (MQ=255)
                                                                                                         aaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctctc                                           <  3:408751/149‑1 (MQ=255)
                                                                                                          aaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctctct                                          <  2:223853/149‑1 (MQ=255)
                                                                                                               gTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACa                                     <  2:497485/149‑1 (MQ=255)
                                                                                                                tAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAg                                    <  3:843920/149‑1 (MQ=255)
                                                                                                                   cacCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTa                                 >  1:62385/1‑149 (MQ=255)
                                                                                                                    acCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaa                                >  1:869709/1‑149 (MQ=255)
                                                                                                                      cGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAg                              >  3:833632/1‑149 (MQ=255)
                                                                                                                      cGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAg                              <  3:877193/149‑1 (MQ=255)
                                                                                                                       gTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGCTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAgg                             >  3:77943/1‑149 (MQ=255)
                                                                                                                             tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTa                       >  2:440876/1‑149 (MQ=255)
                                                                                                                             tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTa                       >  3:283312/1‑149 (MQ=255)
                                                                                                                                   cGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGaaa                 <  4:107046/149‑1 (MQ=255)
                                                                                                                                    ggTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAAc                <  4:55202/149‑1 (MQ=255)
                                                                                                                                       aCGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCAAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt             <  4:19818/149‑1 (MQ=255)
                                                                                                                                        cGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    <  1:565349/45‑1 (MQ=25)
                                                                                                                                        cGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAata                                                                                                                    >  2:565349/1‑45 (MQ=25)
                                                                                                                                            cGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcgtc        <  1:203648/149‑1 (MQ=255)
                                                                                                                                                  aCGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGAtc  >  2:691001/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATT‑CTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATC  >  NC_000913/2939189‑2939482

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