Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,937,440:1 +C 7.2% intergenic (+81/‑28) fucI → / → fucK L‑fucose isomerase/L‑fuculokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,937,4401.C 7.2% 441.0 / 28.5 111intergenic (+81/‑28)fucI/fucKL‑fucose isomerase/L‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base . (58/45);  new base C (4/4);  total (62/49)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 116 total aligned reads displayed.

CACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGAT  >  NC_000913/2937291‑2937589
                                                                                                                                                      |                                                                                                                                                     
cagtgCATGGATATTGAAGTCCAGGAGTTCCGCGATTGCCAGAACTACTGTAGGTTGTACAAGCTTTTCTTCATTACTCCCTGATGTTTTGCCCGGTCGATCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATg                                                                                                                                                        <  1:855543/146‑1 (MQ=255)
cACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTATTACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATg                                                                                                                                                        >  2:342182/1‑150 (MQ=255)
 aCGGCATGGATATTGAAGGCCAGGAGTAACGCGCTTGCCAGAACTACGGGCCGTTGTCCAAGCGTTAATACAGTACTCCCTGATGTGCTGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGt                                                                                                                                                      <  1:305247/150‑1 (MQ=255)
  cGGCATGGATATTGCAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTTATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGtt                                                                                                                                                     <  1:1196927/150‑1 (MQ=255)
     cATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATc                                                                                                                                                  >  2:1249374/1‑150 (MQ=255)
        ggaacttGAAGGCCAGGAATACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTATTACTTTACTCCCTGATGTGATGGCCGGTCGGTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAACATTGCTTATCCgg                                                                                                                                               <  1:1596634/145‑1 (MQ=255)
         gATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAACATTGCTTATCCGGc                                                                                                                                              <  1:1131271/150‑1 (MQ=255)
         gATAATGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTACTACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAACATTGCTTATCCGGc                                                                                                                                              <  1:310263/150‑1 (MQ=255)
          atatTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCt                                                                                                                                             >  2:736510/1‑150 (MQ=255)
          atatTGAAGGCCAGGATTACAGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAACATTGCTTATCCGGCt                                                                                                                                             <  2:136251/150‑1 (MQ=255)
             ttGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCtata                                                                                                                                          >  1:568000/1‑150 (MQ=255)
              tGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTatat                                                                                                                                         >  1:750776/1‑150 (MQ=255)
                  ggCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCa                                                                                                                                     >  2:283218/1‑150 (MQ=255)
                   gCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTCTCCGGCTATATTGCAg                                                                                                                                    >  1:167987/1‑150 (MQ=255)
                    ccAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAgg                                                                                                                                   <  2:1016789/150‑1 (MQ=255)
                    ccAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAgg                                                                                                                                   <  2:1557135/150‑1 (MQ=255)
                    ccAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGAGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAgg                                                                                                                                   >  1:1531124/1‑150 (MQ=255)
                       ggATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCTGCTATATTGCAGGAGc                                                                                                                                >  2:1049807/1‑150 (MQ=255)
                       ggATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGc                                                                                                                                >  2:3130/1‑150 (MQ=255)
                        gATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTTATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCg                                                                                                                               <  1:1654469/150‑1 (MQ=255)
                        gATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCg                                                                                                                               >  2:1284992/1‑150 (MQ=255)
                           tACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGAtt                                                                                                                            <  1:647578/150‑1 (MQ=255)
                            aCCGCGCTTGCCAGAACTACGGTCCGTTCTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGATAGCCGGATAAGCAATGTTATCCAGCTATATTGCAGGAGCGATTa                                                                                                                           >  2:1313617/1‑150 (MQ=255)
                                  cTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAAGACATAACTCCCTGATGAGCTGCCCGGTCGCTCCGGCTACCGGGCCTGAAAAAGCAACAGTGGTTAGCCGGATAAGCATTGTTATCCGGCTTTATTGCAGGAGCGATGATGAAAc                                                                                                                     >  1:136188/1‑150 (MQ=255)
                                   ttGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACa                                                                                                                    <  2:1474952/150‑1 (MQ=255)
                                        aGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACaagaag                                                                                                               <  2:1391059/150‑1 (MQ=255)
                                         gAACTACGGCCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGAAGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGt                                                                                                              <  2:182557/150‑1 (MQ=255)
                                            cTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTAt                                                                                                           <  1:264620/150‑1 (MQ=255)
                                            cTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTAt                                                                                                           >  2:1682110/1‑150 (MQ=255)
                                              aCGGTCCGTTGTACAAGCGTTAATACATTACTCCCCGATGTGATGCCCGTTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCAGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATcc                                                                                                         <  1:1135565/150‑1 (MQ=255)
                                               cGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCt                                                                                                        <  2:1151091/150‑1 (MQ=255)
                                                ggTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTg                                                                                                       <  2:1256584/150‑1 (MQ=255)
                                                 gTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTTTCCTgg                                                                                                      >  2:782738/1‑150 (MQ=255)
                                                 gTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTgg                                                                                                      >  2:1199653/1‑150 (MQ=255)
                                                          aCAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACt                                                                                             <  2:862855/150‑1 (MQ=255)
                                                           cAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGTTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACtg                                                                                            >  1:1586646/1‑150 (MQ=255)
                                                             aGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACtgtg                                                                                          <  2:1405256/150‑1 (MQ=255)
                                                              gCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTgg                                                                                         >  1:729105/1‑150 (MQ=255)
                                                              gCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTgg                                                                                         <  2:364214/150‑1 (MQ=255)
                                                               cGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGTAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGgc                                                                                        >  1:222766/1‑150 (MQ=255)
                                                               cGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGgc                                                                                        >  1:126933/1‑150 (MQ=255)
                                                                ggtAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGgcg                                                                                       <  1:1475579/148‑1 (MQ=255)
                                                                 ttAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGgcgc                                                                                      <  1:144261/150‑1 (MQ=255)
                                                                 ttAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGgcgc                                                                                      >  2:228542/1‑150 (MQ=255)
                                                                     tACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGAcc                                                                                  <  2:391738/150‑1 (MQ=255)
                                                                       cATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCaa                                                                                <  1:896120/150‑1 (MQ=255)
                                                                         ttACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATg                                                                              <  2:960594/150‑1 (MQ=255)
                                                                         ttACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACCGTGGCGCGACCAATg                                                                              <  2:167987/150‑1 (MQ=255)
                                                                              cccTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGAGCGCCCAATGTCAgg                                                                         >  2:1709155/1‑150 (MQ=255)
                                                                               ccTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAggg                                                                        <  1:1733542/150‑1 (MQ=255)
                                                                               ccTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAggg                                                                        <  1:342182/150‑1 (MQ=255)
                                                                                   atgtgatgCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCAt                                                                    >  1:1678511/1‑150 (MQ=255)
                                                                                      tgatgCCCGGTCGCTCCGTCTACGGGGCGTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATcgc                                                                 <  2:436881/150‑1 (MQ=255)
                                                                                       gatgCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATcgcg                                                                >  1:221901/1‑150 (MQ=255)
                                                                                       gatgCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATcgcg                                                                <  1:1612607/150‑1 (MQ=255)
                                                                                          gCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGGCCAATGTCAGGGCCATCGCGGtt                                                             >  2:1203806/1‑150 (MQ=255)
                                                                                            ccGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTaa                                                           >  2:712210/1‑150 (MQ=255)
                                                                                               gTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCg                                                        <  1:788627/150‑1 (MQ=255)
                                                                                                    tCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAgg                                                   <  2:1638660/150‑1 (MQ=255)
                                                                                                    gccGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAgg                                                   <  2:136659/149‑1 (MQ=255)
                                                                                                    gccGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAgg                                                   <  2:1537885/149‑1 (MQ=255)
                                                                                                         cTACCGGGACTGAACAATAAATAGTGGTTAGCCGGATAAGCAATGTTATACGGCTATATTGCAGTATGGATTTTGAAAAAAGATGTTAACCTGGTAATCTACTGTGGCGCGACCAATGTCAGGCACAACGCGTTTAAACGGCAGGGCaaa                                              >  2:69224/1‑150 (MQ=255)
                                                                                                           aCCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAAt                                            >  2:447341/1‑150 (MQ=255)
                                                                                                           aCCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAAt                                            >  2:1508682/1‑150 (MQ=255)
                                                                                                             cGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAAttg                                          >  2:1465254/1‑150 (MQ=255)
                                                                                                                gCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGAAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAAttgttg                                       <  2:1250527/150‑1 (MQ=255)
                                                                                                                 ccTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGc                                      <  1:1133124/150‑1 (MQ=255)
                                                                                                                   tGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGccc                                    >  1:956042/1‑150 (MQ=255)
                                                                                                                    gAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCcg                                   >  2:1187643/1‑150 (MQ=255)
                                                                                                                    gAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCcg                                   >  2:679985/1‑150 (MQ=255)
                                                                                                                    gAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCcg                                   >  1:940550/1‑150 (MQ=255)
                                                                                                                     aaCAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCcgc                                  >  2:767786/1‑150 (MQ=255)
                                                                                                                     aaCAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCcgc                                  >  2:4479/1‑150 (MQ=255)
                                                                                                                     aaCAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCcgc                                  >  2:781783/1‑150 (MQ=255)
                                                                                                                      aCAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAACAATTTTGccctgc                                 <  2:933312/150‑4 (MQ=255)
                                                                                                                         atcaagAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCt                              <  2:1092581/148‑1 (MQ=255)
                                                                                                                         agcaagAGTGGTTAGCCGGCTAAGCAATGTTATCCGGCTATATTGCAGGGGCGATTTTGAAACAAGAAGTTATACTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCACAATTGTTGCCCGCGCCt                              >  1:213329/1‑150 (MQ=255)
                                                                                                                         agcaagAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCACCt                              >  1:646488/1‑150 (MQ=255)
                                                                                                                            aagAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCaa                           >  2:391282/1‑150 (MQ=255)
                                                                                                                                  ggTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTa                     >  2:1030057/1‑150 (MQ=255)
                                                                                                                                  ggTTAGCCGGATAAGCAATGTTATCCGGCTATATTGAAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCACCGCCTa                     >  2:1062484/1‑150 (MQ=255)
                                                                                                                                   gTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTaa                    <  2:729105/150‑1 (MQ=255)
                                                                                                                                    ttAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAAt                   <  1:283218/150‑1 (MQ=255)
                                                                                                                                     tAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGGCCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATg                  <  1:1203806/150‑1 (MQ=255)
                                                                                                                                     tAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCAGCGCCTCAACGCCTAATg                  >  1:1350637/1‑150 (MQ=255)
                                                                                                                                      aGCCGGATAACATTGCTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGc                 >  2:1676779/1‑150 (MQ=255)
                                                                                                                                       gCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGcc                <  1:117970/150‑1 (MQ=255)
                                                                                                                                       gCCGGATAACATTGCTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGcc                >  2:1263828/1‑150 (MQ=255)
                                                                                                                                         cGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAg              <  1:879997/150‑1 (MQ=255)
                                                                                                                                          ggATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGc             >  1:768793/1‑150 (MQ=255)
                                                                                                                                          ggATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGc             >  1:1273863/1‑150 (MQ=255)
                                                                                                                                          ggATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGc             <  2:890848/150‑1 (MQ=255)
                                                                                                                                          ggATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGc             >  2:1110667/1‑150 (MQ=255)
                                                                                                                                                gCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATc       >  1:371426/1‑150 (MQ=255)
                                                                                                                                                 cAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATcg      <  2:1676367/150‑1 (MQ=255)
                                                                                                                                                  aaTGTTATCCGGCTATATTGCAGGAGCGATTTTGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCACCGCCTAATGCCAGCGATATcgc     >  2:235601/1‑150 (MQ=255)
                                                                                                                                                  aaTGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCGCAACGCCTAATGCCAGCGATATcgc     <  1:1199653/150‑1 (MQ=255)
                                                                                                                                                  aaTGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCATAGTACTCGACTGTCGCGCGACCAATTTCAGGGCCATCGCGGATAATCGCCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCAAGCGATATcgc     >  2:1610288/1‑150 (MQ=255)
                                                                                                                                                    tGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGa   >  2:440417/1‑150 (MQ=255)
                                                                                                                                                     gTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGAt  <  1:1249374/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CACGGCATGGATATTGAAGGCCAGGATTACCGCGCTTGCCAGAACTACGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGAT  >  NC_000913/2937291‑2937589

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: