Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,704 G→T 8.0% intergenic (‑84/+21) yhfG ← / ← ppiA putative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,491,7040GT 8.0% 223.9 / 8.2 85intergenic (‑84/+21)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base G (25/51);  new base T (3/4);  total (30/55)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

AAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACC  >  NC_000913/3491559‑3491850
                                                                                                                                                  |                                                                                                                                                  
aaaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCGCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACTCCGGGAGAGCAGAAGCGCGGGCAAGATTg                                                                                                                                                <  1:1417576/150‑1 (MQ=255)
aaaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTg                                                                                                                                                <  1:167583/150‑1 (MQ=255)
    tacgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCcg                                                                                                                                            <  1:81115/150‑1 (MQ=255)
      cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGTGGGCAAGATTGCCcgcg                                                                                                                                          >  1:819775/1‑150 (MQ=255)
      cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACATGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcg                                                                                                                                          <  1:429000/150‑1 (MQ=255)
       gattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAAAATTGccgcgc                                                                                                                                          >  2:1095198/1‑146 (MQ=255)
        attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgct                                                                                                                                        <  2:1061757/150‑2 (MQ=255)
        attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgct                                                                                                                                        <  1:1270234/150‑2 (MQ=255)
        attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCcgcgcg                                                                                                                                        >  1:1196510/1‑150 (MQ=255)
          tacgCTGAAGCTCCCAGAGTCGGGACTTTTGCTTATCTGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCgag                                                                                                                                      <  2:1684458/150‑1 (MQ=255)
           acgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATGGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCgaga                                                                                                                                     <  2:389692/150‑1 (MQ=255)
            cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGaa                                                                                                                                    <  1:1564429/150‑1 (MQ=255)
            cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGaa                                                                                                                                    >  2:1161824/1‑150 (MQ=255)
            cgCTGAAGCGCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAATCTTGCCcgcgcttctg                                                                                                                                    <  2:1201650/150‑6 (MQ=255)
               tGAAGCTCCCAGAGACGGTACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGCACGAGTCCATTTGCCACAAGTAGAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGCGCGGGGCAGATTGCCAGCGCGAGtaaag                                                                                                                                 >  1:256915/1‑148 (MQ=255)
                gAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCa                                                                                                                                >  2:1139679/1‑150 (MQ=255)
                 aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                               <  1:903059/150‑1 (MQ=255)
                    cTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGGCAGCTTTGCCGCGGGCGAAATCATTaa                                                                                                                            >  2:1595399/1‑149 (MQ=255)
                    cTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTAc                                                                                                                            >  2:1020443/1‑150 (MQ=255)
                       ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGc                                                                                                                         <  2:900815/150‑1 (MQ=255)
                        cAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTACACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCAGGGTTACGCCGGTAGAGCAGAAGCGCTGGCAATATTGCCCGCCCCACAAATAAATACAGCa                                                                                                                        >  1:1109257/1‑150 (MQ=255)
                         agagACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAg                                                                                                                       <  2:7740/150‑1 (MQ=255)
                            gaCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGAc                                                                                                                    <  2:41533/150‑1 (MQ=255)
                              cGGGTCTTTTGCTTTTCGGTGAGTTTCTTCACTTGATGCCTCGCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACtt                                                                                                                  <  2:1492760/150‑1 (MQ=255)
                                  aCTTTTGCTTATCGGTGAGTTTCTTCCCTTGATGCCTCCCTGAACGTGTCCCTTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGCGCGGGCAAGATAGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGc                                                                                                              <  2:97394/150‑1 (MQ=255)
                                   cTTTTGCTTATCGGTGCGCTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTCCCACAAGAATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                             <  2:1082891/150‑1 (MQ=255)
                                    ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCACTTCCGGGTTACGCCGGGAGAGCAGAAGCGCCGGCAAGATTGCCCGCGCGAGAAATGAATACCGCAGGACAGTAGCgg                                                                                                            >  1:119260/1‑150 (MQ=255)
                                     tttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGa                                                                                                           <  1:1305107/150‑1 (MQ=255)
                                       tGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGTGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGaaa                                                                                                         <  2:819775/150‑1 (MQ=255)
                                         cTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTTGCCGCGCGAGAAATCAATACGGCACGCCGTTTGCGGGAAgg                                                                                                       >  2:1614505/1‑150 (MQ=255)
                                           tATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTGTAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGAt                                                                                                     <  2:545468/150‑1 (MQ=255)
                                              cGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAc                                                                                                  <  2:537627/150‑1 (MQ=255)
                                                gTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTa                                                                                                <  1:452669/150‑1 (MQ=255)
                                                  gAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTAcc                                                                                              <  1:1098120/150‑1 (MQ=255)
                                                         ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGGCCGCGCGCGCAATCATTACGGCAGGCCATTAGCGGAAATGATTACTACCGGgtttt                                                                                       >  1:907602/1‑149 (MQ=255)
                                                         ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTg                                                                                       <  2:1673349/150‑1 (MQ=255)
                                                           cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAc                                                                                     <  2:1056737/150‑1 (MQ=255)
                                                                  tGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACAATGGCCCCGCGAGAAATCAATAACGCAGGAATTTAGCGGAAAGGATAACAAAAGGGTTTTACGGcacaa                                                                              >  2:431388/1‑149 (MQ=255)
                                                                     cTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCt                                                                           <  2:1139746/150‑1 (MQ=255)
                                                                       cccTGAACGTGTCCAATTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                         <  1:862210/150‑1 (MQ=255)
                                                                        ccTGAACGTGTCCATTTGCCACAAGTATTAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCCGGAGAAATCATTAAGGCAGGACGTTAGCGGTAAGGATAACTACCGGTTTTGAAGGCACCTTCTGGt                                                                        >  1:104393/1‑150 (MQ=255)
                                                                         cTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATAATTACGGCAGGACTTTAGCGGAAAGGATAACAACCGGTTTTGACGGCCCATTCTGGtt                                                                       >  2:1173483/1‑149 (MQ=255)
                                                                            aaCGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACgg                                                                    >  1:1017640/1‑150 (MQ=255)
                                                                             tcttgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGATAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGa                                                                   <  1:1695724/147‑1 (MQ=255)
                                                                              cGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAc                                                                  <  1:689917/150‑1 (MQ=255)
                                                                               gtgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTTGGGGAAAGGATAACTACCGGTTTTGGCGGCACCTTCTGGGACGGAcc                                                                 >  1:898552/1‑150 (MQ=255)
                                                                               gtgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAcc                                                                 <  2:405165/150‑1 (MQ=255)
                                                                                         gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGa                                                       <  1:1391350/150‑1 (MQ=255)
                                                                                         gCCACAAGTATAAGCGGCAAATCCGGGTCACGTCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCAAATTCTGGTACGGACCAACGTCATGa                                                       <  1:1190281/150‑1 (MQ=255)
                                                                                          ccACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCACGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTTCGGACCAACGTCATGAg                                                      >  2:1560367/1‑150 (MQ=255)
                                                                                            ccaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTc                                                    <  1:1095198/149‑1 (MQ=255)
                                                                                                 tataAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCAc                                               <  2:1416170/150‑1 (MQ=255)
                                                                                                     aGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTgg                                           <  1:1020443/150‑1 (MQ=255)
                                                                                                      gCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGGGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTCCCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTggg                                          <  1:832036/150‑1 (MQ=255)
                                                                                                      gCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTggg                                          <  1:611511/150‑1 (MQ=255)
                                                                                                          cAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGAATGCCCGCGCGAGAAATCATTACGGCAGGACTTAAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAAt                                      <  1:1518005/150‑1 (MQ=255)
                                                                                                            aaTCCGGGTTACGACGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                    <  1:689495/150‑1 (MQ=255)
                                                                                                                cGGGTTACGCCGGGAGAGCAGAAGCACGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTTCGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                <  1:1560367/150‑1 (MQ=255)
                                                                                                                cGGGGTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                <  2:138084/150‑1 (MQ=255)
                                                                                                                 gtgTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCg                               <  1:1324938/148‑1 (MQ=255)
                                                                                                                 ggggTATGCCGGGAGAGCAGAAGCGCGGGCAAGATGGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAATGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCg                               <  2:1150256/150‑1 (MQ=255)
                                                                                                                  gggtACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACATTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                              <  1:466976/147‑1 (MQ=255)
                                                                                                                  ggTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                              <  1:1161824/150‑1 (MQ=255)
                                                                                                                  ggTTACGCCGGGAGAGCAGAAGCGCGGGAAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                              <  2:438869/150‑1 (MQ=255)
                                                                                                                        gCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTCCGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGt                        <  2:1509804/150‑1 (MQ=255)
                                                                                                                           gggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCa                     <  2:1283644/150‑1 (MQ=255)
                                                                                                                            ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                    <  2:1618888/150‑1 (MQ=255)
                                                                                                                             gagaGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTCGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATg                   <  1:1139679/150‑1 (MQ=255)
                                                                                                                              cgagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTGTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                  <  2:352648/149‑1 (MQ=255)
                                                                                                                               gagCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGcc                 <  1:431388/150‑1 (MQ=255)
                                                                                                                                  cAGAAGCGCGGGCAAGATTGCCCGCGCGCGAAATAAATAAAGCACGACATTAGCGGAAAAGATAACTAACGGTTTTGACGGCACAATATTGTTAGGACCCACGGCAAGAGTCCGGACCAGGGGAATCTTTTAGGCAACGTCCATGCCttt              >  2:1094169/1‑150 (MQ=255)
                                                                                                                                  cAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCttt              <  2:119260/150‑1 (MQ=255)
                                                                                                                                    gAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca            <  1:773185/150‑1 (MQ=255)
                                                                                                                                    ctcgcgcGGGCAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca            >  2:840859/4‑150 (MQ=255)
                                                                                                                                      cgcgcGGGCAATCGTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGCCGGCACATTCTGGTACGGCCCAACGGCAGGAGTTGGCACCAGGGAAATCGTATCGGCCCAGACCATGCCTTTcccc          >  2:736157/2‑150 (MQ=255)
                                                                                                                                        cgcgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                        >  2:309111/1‑102 (MQ=255)
                                                                                                                                        cgcgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                        <  1:309111/102‑1 (MQ=255)
                                                                                                                                         gcgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACt       >  1:217670/1‑150 (MQ=255)
                                                                                                                                         gcgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTGTCACCACt       >  1:1090731/1‑150 (MQ=255)
                                                                                                                                          cgGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGACACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCCGCAACGTCCATGCCTTTCACCACtt      >  1:1132360/1‑150 (MQ=255)
                                                                                                                                            ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa    >  1:197174/1‑150 (MQ=255)
                                                                                                                                            ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa    >  1:1019163/1‑150 (MQ=255)
                                                                                                                                             gCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc   >  2:1257487/1‑150 (MQ=255)
                                                                                                                                              gctaGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGGGTTGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc   >  2:1410000/4‑150 (MQ=255)
                                                                                                                                               cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc  >  1:1661155/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACC  >  NC_000913/3491559‑3491850

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: