Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,476,585 G→A 6.0% intergenic (‑145/+22) 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,5850GA 6.0% 377.1 / 7.0 115intergenic (‑145/+22)yheO/fkpAputative PAS domain‑containing DNA‑binding transcriptional regulator/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base G (56/52);  new base A (2/5);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGG  >  NC_000913/3476440‑3476720
                                                                                                                                                 |                                                                                                                                       
tGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCa                                                                                                                                      >  4:66044/1‑149 (MQ=255)
 gAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAg                                                                                                                                     >  2:168580/1‑149 (MQ=255)
   aaaaaaaCACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAggc                                                                                                                                   >  2:323474/1‑149 (MQ=255)
   aaaaaaaCACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATACGGGAAACTACATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTCATAAAAATACACCCGCCTGATCAggc                                                                                                                                   >  3:200010/1‑149 (MQ=255)
      aaaaCACTCCTTTTCAGGAGCCTGTCGTTAACTTTTGAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAggcggc                                                                                                                                <  3:701786/149‑1 (MQ=255)
           acTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGGCGGCGGttt                                                                                                                           <  1:68767/149‑2 (MQ=255)
           acTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAAGTATATCCGGCCTGTAATAAACAAACCGCCGCCTGACCAGGCGGCGGttt                                                                                                                           <  2:481452/149‑2 (MQ=255)
            cTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCt                                                                                                                          <  3:212467/149‑1 (MQ=255)
             tCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCtt                                                                                                                         >  1:189726/1‑149 (MQ=255)
             tCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGGCGGCGGTtttt                                                                                                                         <  3:486448/149‑1 (MQ=255)
              ccTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTa                                                                                                                        >  2:797301/1‑149 (MQ=255)
               cTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCCAGCGCTTTAATGATGCGAGTAACTAAATCAATACAGGGGAGGATCCGCTTTCCAGCACTAACTATATCCAGCCTGTAATAAAAAAACCGCCGCCAGGTCAGCCTGCCGATCttca                                                                                                                       >  3:838475/1‑147 (MQ=255)
               cTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTaa                                                                                                                       >  1:266443/1‑149 (MQ=255)
               cTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTaa                                                                                                                       <  1:674984/149‑1 (MQ=255)
               cTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGGCGGCGGTtttttt                                                                                                                       <  1:236581/149‑3 (MQ=255)
               cTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGACGTTCCGCTTTCCAGCACTAATAATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGGGGCGGTTCTTaa                                                                                                                       >  2:1034975/1‑149 (MQ=255)
                ttttCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAAt                                                                                                                      >  3:604/1‑149 (MQ=255)
                ttttCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAAt                                                                                                                      >  4:680143/1‑149 (MQ=255)
                ttttCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAAt                                                                                                                      >  3:860605/1‑149 (MQ=255)
                  ttCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGc                                                                                                                    <  2:1000494/149‑1 (MQ=255)
                      ggAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAt                                                                                                                >  1:32784/1‑149 (MQ=255)
                           cTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAtttttt                                                                                                           >  2:708672/1‑149 (MQ=255)
                           cTGTCGTTAACTTTTCAGGGCAGGCTCATCAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAtttttt                                                                                                           <  4:634837/149‑1 (MQ=255)
                           cTGTCGTTAACTTTTCAGGGCAGGCTAATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAGACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAtttttt                                                                                                           <  2:887388/149‑1 (MQ=255)
                            tGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa                                                                                                          <  3:473287/149‑1 (MQ=255)
                               cGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCa                                                                                                       <  4:3769/149‑1 (MQ=255)
                               cGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCa                                                                                                       <  2:83315/149‑1 (MQ=255)
                                 ttAACTTTTCAGGGCAGGCTCATTAATGATGCGGGGAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGAGGTTCTTAATGCTTATTTTTTAGCAGa                                                                                                     >  4:93137/1‑149 (MQ=255)
                                   aaCTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAAt                                                                                                   >  3:362193/1‑149 (MQ=255)
                                   aaCTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAAt                                                                                                   <  4:618197/149‑1 (MQ=255)
                                   aaCTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAAt                                                                                                   <  4:359064/149‑1 (MQ=255)
                                      ttttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGc                                                                                                        <  2:45677/141‑1 (MQ=255)
                                      ttttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTg                                                                                                >  4:565776/1‑149 (MQ=255)
                                        ttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCg                                                                                              >  4:386275/1‑149 (MQ=255)
                                            gggCAGGCTCATTAATGATGCGGGTAACTGAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCtt                                                                                          >  3:994612/1‑149 (MQ=255)
                                            gggCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa                                                                                                          >  1:900418/1‑133 (MQ=255)
                                            gggCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa                                                                                                          <  2:900418/133‑1 (MQ=255)
                                              gCAGGCTCATTAATGATGCGGGTAACTAAATTATTACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTc                                                                                        <  3:866529/149‑1 (MQ=255)
                                                  gCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCAt                                                                                    <  4:604/149‑1 (MQ=255)
                                                      aTTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCCGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGc                                                                                >  1:342517/1‑149 (MQ=255)
                                                      aTTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGc                                                                                >  1:708628/1‑149 (MQ=255)
                                                      aTTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGc                                                                                >  1:126399/1‑149 (MQ=255)
                                                          atgatgCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTcc                                                                            >  1:124272/1‑149 (MQ=255)
                                                            gatgCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCgg                                                                          >  4:535506/1‑149 (MQ=255)
                                                             atgCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGc                                                                         >  1:160756/1‑149 (MQ=255)
                                                              tgCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCt                                                                        <  4:226831/149‑1 (MQ=255)
                                                               gCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCtt                                                                       <  2:189726/149‑1 (MQ=255)
                                                                 gggTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCttta                                                                     <  2:32784/149‑2 (MQ=255)
                                                                    tAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCAt                                                                  >  2:631111/1‑149 (MQ=255)
                                                                        tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTaa                                                                                                                       >  3:951133/1‑92 (MQ=255)
                                                                        tAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCAGGCCTGTAATAAAAAAACCGACGCCTGGTCAGGCGGCGGTCCTTAATGCTTCTTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGc                                                              >  2:991501/1‑149 (MQ=255)
                                                                          aaTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTc                                                                                                                           <  1:1005261/86‑1 (MQ=255)
                                                                          aaTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTc                                                                                                                           >  2:1005261/1‑86 (MQ=255)
                                                                               aTACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGgc                                                       >  2:617344/1‑149 (MQ=255)
                                                                                   aGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTg                                                   <  1:797301/149‑1 (MQ=255)
                                                                                       gAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCtt                                                                       <  4:1025970/125‑1 (MQ=255)
                                                                                       gAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCtt                                                                       >  3:1025970/1‑125 (MQ=255)
                                                                                        aGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTc                                              <  3:535506/149‑1 (MQ=255)
                                                                                        aGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTc                                              <  4:136933/149‑1 (MQ=255)
                                                                                          gTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTcac                                            <  1:617344/149‑1 (MQ=255)
                                                                                           ttCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTcaca                                           >  3:99168/1‑149 (MQ=255)
                                                                                                cTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcag                                      >  1:991162/1‑149 (MQ=255)
                                                                                                 tttCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcagc                                     >  2:559713/1‑149 (MQ=255)
                                                                                                     cAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGct                                 >  3:497999/1‑149 (MQ=255)
                                                                                                     cAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGct                                 <  2:160756/149‑1 (MQ=255)
                                                                                                       gCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGctct                               <  2:226533/149‑1 (MQ=255)
                                                                                                         aCTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTAc                             >  3:507706/1‑149 (MQ=255)
                                                                                                         aCTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTAc                             >  4:655288/1‑149 (MQ=255)
                                                                                                           tAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGt                           <  3:838961/149‑1 (MQ=255)
                                                                                                             aTTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCa                         >  1:785803/1‑149 (MQ=255)
                                                                                                               tataTCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCTGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTGCGCATCATCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaaa                       <  3:1015914/149‑1 (MQ=255)
                                                                                                               tataTCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaaa                       <  1:168580/149‑1 (MQ=255)
                                                                                                                atatCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAac                      <  4:838475/149‑1 (MQ=255)
                                                                                                                atatCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAac                      <  1:708672/149‑1 (MQ=255)
                                                                                                                  atCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCt                                                                        >  2:460581/1‑97 (MQ=255)
                                                                                                                  atCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCt                                                                        <  1:460581/97‑1 (MQ=255)
                                                                                                                    ccGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCa                  <  4:497999/149‑1 (MQ=255)
                                                                                                                     cGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAg                 <  2:825323/149‑1 (MQ=255)
                                                                                                                      ggCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAgg                >  3:157840/1‑149 (MQ=255)
                                                                                                                         cTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTa             >  4:527108/1‑149 (MQ=255)
                                                                                                                            taataaAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGaa          >  1:1039716/1‑149 (MQ=255)
                                                                                                                               taaAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTc       >  4:735098/1‑149 (MQ=255)
                                                                                                                                 aaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGt                                                 >  3:816159/1‑105 (MQ=255)
                                                                                                                                 aaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGt                                                 <  4:816159/105‑1 (MQ=255)
                                                                                                                                 aaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCgg     >  1:293326/1‑149 (MQ=255)
                                                                                                                                 aaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCgg     <  1:309908/149‑1 (MQ=255)
                                                                                                                                 aaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCgg     >  4:568679/1‑149 (MQ=255)
                                                                                                                                  aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGt    >  4:428004/1‑149 (MQ=255)
                                                                                                                                  aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGt    <  2:311779/149‑1 (MQ=255)
                                                                                                                                   aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg   >  4:216484/1‑149 (MQ=255)
                                                                                                                                   aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg   >  2:924662/1‑149 (MQ=255)
                                                                                                                                   aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg   <  3:565776/149‑1 (MQ=255)
                                                                                                                                   aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg   <  3:655288/149‑1 (MQ=255)
                                                                                                                                   aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg   <  3:66044/149‑1 (MQ=255)
                                                                                                                                    aaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTgg  <  4:941247/149‑1 (MQ=255)
                                                                                                                                     aaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCAt                                                                  >  1:441328/1‑84 (MQ=255)
                                                                                                                                        cgccgcCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTa             <  4:77039/134‑1 (MQ=255)
                                                                                                                                        cgccgcCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTa             >  3:77039/1‑134 (MQ=255)
                                                                                                                                         gccgccTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg   >  3:744340/1‑143 (MQ=255)
                                                                                                                                          ccgccTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcag                                      >  4:83806/1‑107 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
TGAAAAAAAACACTCCTTTTCAGGAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGG  >  NC_000913/3476440‑3476720

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