Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,476,586 T→C 6.5% intergenic (‑146/+21) yheO ← / ← fkpA putative PAS domain‑containing DNA‑binding transcriptional regulator/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,476,5860TC 6.5% 698.5 / 25.9 215intergenic (‑146/+21)yheO/fkpAputative PAS domain‑containing DNA‑binding transcriptional regulator/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (105/94);  new base C (12/2);  total (119/96)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 217 total aligned reads displayed.

CATGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACAC  >  NC_000913/3476438‑3476734
                                                                                                                                                    |                                                                                                                                                    
cATGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTc                                                                                                                                                     <  2:1529551/150‑1 (MQ=255)
  tGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCTGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAg                                                                                                                                                   <  2:2247548/150‑1 (MQ=255)
    aaaaaaaaCACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGAGGAGGTTCCGCTTTACAGCACTAATTATATCCGGCCTGTAATAAAAAATGCGCCGCCTGGTGAggc                                                                                                                                                 >  2:1133821/1‑150 (MQ=255)
        aaaaCACTCCTTTTCAGGAGCCTGTCGTGAACTTTTCAGGGCAGGCGCCTTAATTATGCGGGGAACTAAAATAATAAAGCGGCGGTTCCGCGTTACAGCACTAATTATATACGGCCGGTAATAAAAAAACCGCCGCCGGGTCAggcgggg                                                                                                                                             >  1:1445284/1‑148 (MQ=255)
            cacTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTc                                                                                                                                         <  2:1031680/150‑1 (MQ=255)
             acTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCt                                                                                                                                        <  2:3711061/150‑1 (MQ=255)
               tCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGACTGTAATAAAAAAACCGCCGCCGGGCCAGGCGGCGGTTCTTa                                                                                                                                      >  2:2335409/1‑150 (MQ=255)
                           gCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAACTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAttttt                                                                                                                          <  2:2200661/150‑1 (MQ=255)
                                   ttAACTTTACAGGGCAGGATAATTAATGATGCGGGTAAATAAATTAATACAGCGTAGGTTCCGCTTTCCAGAACTAATTATAACCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGaa                                                                                                                  <  1:1325314/150‑1 (MQ=255)
                                    tAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTCATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAAt                                                                                                                 <  2:2105668/150‑1 (MQ=255)
                                     aaCTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATc                                                                                                                >  2:1261842/1‑150 (MQ=255)
                                     aaCTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATc                                                                                                                >  1:3752431/1‑150 (MQ=255)
                                     aaCTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCAGGCCTGTAATAAAAAAACCGCAGCCTGGTCAGGCGGAGGTAATTAATGCTTATTTTTTAGCAGAATc                                                                                                                >  2:2824143/1‑150 (MQ=255)
                                      aCTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCt                                                                                                               >  1:891340/1‑150 (MQ=255)
                                        ttttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGc                                                                                                             >  2:2350640/1‑150 (MQ=255)
                                         tttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACCGCGGAGGTTCCGCTTTCCAGCACTAAGTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCg                                                                                                            <  2:705328/150‑1 (MQ=255)
                                         tttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCg                                                                                                            <  1:1649001/150‑1 (MQ=255)
                                         tttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCg                                                                                                            <  2:326414/150‑1 (MQ=255)
                                          ttCTGGGCTGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGTTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCgg                                                                                                           <  1:2411323/150‑1 (MQ=255)
                                          ttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCgg                                                                                                           >  2:1253661/1‑150 (MQ=255)
                                           tCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGc                                                                                                          >  2:1327810/1‑150 (MQ=255)
                                            cAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCt                                                                                                         >  1:1475028/1‑150 (MQ=255)
                                              gggCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCttt                                                                                                       <  2:177678/150‑1 (MQ=255)
                                              gggCAGGCTCATTAATGATGCGGGAAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCttt                                                                                                       <  2:1915262/150‑1 (MQ=255)
                                               ggCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTc                                                                                                      <  1:611864/150‑1 (MQ=255)
                                                 cAGGCTCATTAATGATGCGGGTAACTAAATTAATACCGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTTTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGc                                                                                                    <  1:285134/150‑1 (MQ=255)
                                                  aGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCa                                                                                                   >  2:367326/1‑150 (MQ=255)
                                                     cTCATTAATGATGCGGGTAACTAAATTAATACCGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCa                                                                                                <  1:2950675/150‑1 (MQ=255)
                                                     cTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCa                                                                                                <  2:2906236/150‑1 (MQ=255)
                                                     cTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCa                                                                                                <  2:1430721/150‑1 (MQ=255)
                                                     cTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCa                                                                                                <  2:3736898/150‑1 (MQ=255)
                                                     cTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCa                                                                                                <  2:395251/150‑1 (MQ=255)
                                                     cTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCa                                                                                                <  2:1987000/150‑1 (MQ=255)
                                                         ttAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCtt                                                                                            >  1:216591/1‑150 (MQ=255)
                                                           aaTGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTcc                                                                                          >  1:500960/1‑150 (MQ=255)
                                                             tgatgCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCgg                                                                                        <  2:3615001/150‑1 (MQ=255)
                                                                       aaCTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTATTTATATCCGGCATGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCa                                                                              <  1:3638904/150‑1 (MQ=255)
                                                                       aaCTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTTTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCa                                                                              <  2:1144723/150‑1 (MQ=255)
                                                                       aaCTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCa                                                                              <  2:3760593/150‑1 (MQ=255)
                                                                       aaCTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCa                                                                              <  2:1198955/150‑1 (MQ=255)
                                                                       aaCTAAATTAATACAGCGGAGGTACCGCTTTCCACCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCa                                                                              <  2:1026540/150‑1 (MQ=255)
                                                                        aCTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAg                                                                             <  2:833026/150‑1 (MQ=255)
                                                                          tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGcc                                                                           >  1:2690704/1‑150 (MQ=255)
                                                                          tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGcc                                                                           <  1:1658275/150‑1 (MQ=255)
                                                                          tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGcc                                                                           >  1:2093993/1‑150 (MQ=255)
                                                                          tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGcc                                                                           >  2:1145051/1‑150 (MQ=255)
                                                                           aaaTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCt                                                                          >  1:3214832/1‑150 (MQ=255)
                                                                             aTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTc                                                                        >  1:3025421/1‑150 (MQ=255)
                                                                                  tACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGgcgc                                                                   >  2:3714444/1‑150 (MQ=255)
                                                                                    cAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTg                                                                 >  1:59518/1‑150 (MQ=255)
                                                                                      gCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGt                                                               >  2:3000490/1‑150 (MQ=255)
                                                                                      gCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGt                                                               >  2:696415/1‑150 (MQ=255)
                                                                                        ggAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGt                                                                                                                                                      >  1:1519269/1‑61 (MQ=255)
                                                                                        ggAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGttt                                                             <  1:2316389/150‑1 (MQ=255)
                                                                                        ggAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGttt                                                             <  1:2273331/150‑1 (MQ=255)
                                                                                            gTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTcaca                                                         <  1:3925709/150‑1 (MQ=255)
                                                                                             ttCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACAt                                                        >  1:2040381/1‑150 (MQ=255)
                                                                                              tCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATc                                                       >  2:260121/1‑150 (MQ=255)
                                                                                               ccGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATcc                                                      <  1:3345669/150‑1 (MQ=255)
                                                                                               ccGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATcc                                                      <  1:1378700/150‑1 (MQ=255)
                                                                                                 gCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcag                                                    >  1:601324/1‑150 (MQ=255)
                                                                                                 gCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcag                                                    >  2:1586870/1‑150 (MQ=255)
                                                                                                 gCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcag                                                    >  1:3928700/1‑150 (MQ=255)
                                                                                                 gCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcag                                                    >  1:1200862/1‑150 (MQ=255)
                                                                                                     tCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCagcagc                                                <  2:968837/150‑1 (MQ=255)
                                                                                                        aGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGctct                                             >  2:2553113/1‑150 (MQ=255)
                                                                                                          cACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTAc                                           <  1:2567459/150‑1 (MQ=255)
                                                                                                           aCTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACg                                          <  2:667221/150‑1 (MQ=255)
                                                                                                           aCTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACg                                          <  2:19516/150‑1 (MQ=255)
                                                                                                            caaaTTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGt                                         <  2:2404432/148‑1 (MQ=255)
                                                                                                            cTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGt                                         >  1:1030283/1‑150 (MQ=255)
                                                                                                              atttATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTTCGTCa                                       <  1:693189/148‑1 (MQ=255)
                                                                                                              aaTTATAGCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCa                                       <  1:693112/150‑1 (MQ=255)
                                                                                                               aTTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaa                                      >  2:3437887/1‑150 (MQ=255)
                                                                                                               aTTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaa                                      >  2:535265/1‑150 (MQ=255)
                                                                                                               aTTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaa                                      >  2:2528322/1‑150 (MQ=255)
                                                                                                                ttATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaaa                                     <  1:1261931/150‑1 (MQ=255)
                                                                                                                   tatCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAacac                                  >  1:1144778/1‑150 (MQ=255)
                                                                                                                     tCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCa                                >  1:2409086/1‑150 (MQ=255)
                                                                                                                       cGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAgg                              >  2:752341/1‑150 (MQ=255)
                                                                                                                       cGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAgg                              >  1:2570544/1‑150 (MQ=255)
                                                                                                                        ggCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAggg                             >  2:1528974/1‑150 (MQ=255)
                                                                                                                            tGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGa                         <  1:2869843/150‑1 (MQ=255)
                                                                                                                                 taaAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                    >  2:2041106/1‑150 (MQ=255)
                                                                                                                                  ttaagaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCgg                   >  2:3318909/6‑150 (MQ=255)
                                                                                                                                   aaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGt                  <  2:2478550/150‑1 (MQ=255)
                                                                                                                                    aagaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCCGCCTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTgg                >  2:1700923/4‑150 (MQ=255)
                                                                                                                                    aaaatCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg                 >  2:79085/1‑150 (MQ=255)
                                                                                                                                    aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg                 <  1:3874159/150‑1 (MQ=255)
                                                                                                                                    aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATAAGCCATCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg                 >  1:2500858/1‑150 (MQ=255)
                                                                                                                                     aaaaCCGCCGCCTGGTCAGGCGGCGGTTATTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTgg                <  2:2128482/150‑1 (MQ=255)
                                                                                                                                       aaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGa              >  2:2696457/1‑150 (MQ=255)
                                                                                                                                         ccgccgccgGGCCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATc            <  1:1749263/150‑1 (MQ=255)
                                                                                                                                         ccgccgCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATc            >  2:2304563/1‑150 (MQ=255)
                                                                                                                                           gccgccTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATccc          >  2:2982207/1‑150 (MQ=255)
                                                                                                                                            ccgccTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATcccc         <  1:2090145/150‑1 (MQ=255)
                                                                                                                                            ccgccTGGCCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATcccc         >  2:3890976/1‑150 (MQ=255)
                                                                                                                                                cTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGaa     >  1:3908250/1‑150 (MQ=255)
                                                                                                                                                   gTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAacac  <  2:181001/150‑1 (MQ=255)
                                                                                                                                                   gTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAacac  >  2:883147/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACAC  >  NC_000913/3476438‑3476734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: