Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,705 A→C 8.0% intergenic (‑85/+20) 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,7050AC 8.0% 217.1 / 9.1 85intergenic (‑85/+20)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base A (27/51);  new base C (3/4);  total (30/55)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

AAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACC  >  NC_000913/3491559‑3491850
                                                                                                                                                   |                                                                                                                                                 
aaaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCGCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACTCCGGGAGAGCAGAAGCGCGGG‑CAAGATTg                                                                                                                                                <  1:1417576/150‑1 (MQ=255)
aaaTTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTg                                                                                                                                                <  1:167583/150‑1 (MQ=255)
    tacgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCcg                                                                                                                                            <  1:81115/150‑1 (MQ=255)
      cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGTGGG‑CAAGATTGCCcgcg                                                                                                                                          >  1:819775/1‑150 (MQ=255)
      cgattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACATGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAATCTTGCCcgcg                                                                                                                                          <  1:429000/150‑1 (MQ=255)
       gattacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAAAATTGccgcgc                                                                                                                                          >  2:1095198/1‑146 (MQ=255)
        attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAATCTTGCCcgcgct                                                                                                                                        <  2:1061757/150‑2 (MQ=255)
        attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAATCTTGCCcgcgct                                                                                                                                        <  1:1270234/150‑2 (MQ=255)
        attacgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCcgcgcg                                                                                                                                        >  1:1196510/1‑150 (MQ=255)
          tacgCTGAAGCTCCCAGAGTCGGGACTTTTGCTTATCTGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCgag                                                                                                                                      <  2:1684458/150‑1 (MQ=255)
           acgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATGGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCgaga                                                                                                                                     <  2:389692/150‑1 (MQ=255)
            cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGaa                                                                                                                                    <  1:1564429/150‑1 (MQ=255)
            cgCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGaa                                                                                                                                    >  2:1161824/1‑150 (MQ=255)
            cgCTGAAGCGCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAATCTTGCCcgcgcttctg                                                                                                                                    <  2:1201650/150‑6 (MQ=255)
               tGAAGCTCCCAGAGACGGTACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGCACGAGTCCATTTGCCACAAGTAGAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGCGCGGG‑GCAGATTGCCAGCGCGAGtaaag                                                                                                                                 >  1:256915/1‑148 (MQ=255)
                gAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCa                                                                                                                                >  2:1139679/1‑150 (MQ=255)
                 aaGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCAt                                                                                                                               <  1:903059/150‑1 (MQ=255)
                    cTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGG‑GCAGCTTTGCCGCGGGCGAAATCATTaa                                                                                                                            >  2:1595399/1‑149 (MQ=255)
                    cTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTAc                                                                                                                            >  2:1020443/1‑150 (MQ=255)
                       ccAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGc                                                                                                                         <  2:900815/150‑1 (MQ=255)
                        cAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTACACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCAGGGTTACGCCGGTAGAGCAGAAGCGCTGG‑CAATATTGCCCGCCCCACAAATAAATACAGCa                                                                                                                        >  1:1109257/1‑150 (MQ=255)
                         agagACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAg                                                                                                                       <  2:7740/150‑1 (MQ=255)
                            gaCGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGAc                                                                                                                    <  2:41533/150‑1 (MQ=255)
                              cGGGTCTTTTGCTTTTCGGTGAGTTTCTTCACTTGATGCCTCGCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACtt                                                                                                                  <  2:1492760/150‑1 (MQ=255)
                                  aCTTTTGCTTATCGGTGAGTTTCTTCCCTTGATGCCTCCCTGAACGTGTCCCTTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGCGCGGG‑CAAGATAGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGc                                                                                                              <  2:97394/150‑1 (MQ=255)
                                   cTTTTGCTTATCGGTGCGCTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTCCCACAAGAATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                             <  2:1082891/150‑1 (MQ=255)
                                    ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCACTTCCGGGTTACGCCGGGAGAGCAGAAGCGCCGG‑CAAGATTGCCCGCGCGAGAAATGAATACCGCAGGACAGTAGCgg                                                                                                            >  1:119260/1‑150 (MQ=255)
                                     tttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGa                                                                                                           <  1:1305107/150‑1 (MQ=255)
                                       tGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGTGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGaaa                                                                                                         <  2:819775/150‑1 (MQ=255)
                                         cTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTTGCCGCGCGAGAAATCAATACGGCACGCCGTTTGCGGGAAgg                                                                                                       >  2:1614505/1‑150 (MQ=255)
                                           tATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTGTAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGAt                                                                                                     <  2:545468/150‑1 (MQ=255)
                                              cGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAAc                                                                                                  <  2:537627/150‑1 (MQ=255)
                                                gTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTa                                                                                                <  1:452669/150‑1 (MQ=255)
                                                  gAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTAcc                                                                                              <  1:1098120/150‑1 (MQ=255)
                                                         ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGGCCGCGCGCGCAATCATTACGGCAGGCCATTAGCGGAAATGATTACTACCGGgtttt                                                                                       >  1:907602/1‑149 (MQ=255)
                                                         ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTg                                                                                       <  2:1673349/150‑1 (MQ=255)
                                                           cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAc                                                                                     <  2:1056737/150‑1 (MQ=255)
                                                                  tGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAACAATGGCCCCGCGAGAAATCAATAACGCAGGAATTTAGCGGAAAGGATAACAAAAGGGTTTTACGGcacaa                                                                              >  2:431388/1‑149 (MQ=255)
                                                                     cTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCt                                                                           <  2:1139746/150‑1 (MQ=255)
                                                                       cccTGAACGTGTCCAATTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTgg                                                                         <  1:862210/150‑1 (MQ=255)
                                                                        ccTGAACGTGTCCATTTGCCACAAGTATTAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCCGGAGAAATCATTAAGGCAGGACGTTAGCGGTAAGGATAACTACCGGTTTTGAAGGCACCTTCTGGt                                                                        >  1:104393/1‑150 (MQ=255)
                                                                         cTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATAATTACGGCAGGACTTTAGCGGAAAGGATAACAACCGGTTTTGACGGCCCATTCTGGtt                                                                       >  2:1173483/1‑149 (MQ=255)
                                                                            aaCGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACgg                                                                    >  1:1017640/1‑150 (MQ=255)
                                                                             tcttgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGATAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGa                                                                   <  1:1695724/147‑1 (MQ=255)
                                                                              cGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAc                                                                  <  1:689917/150‑1 (MQ=255)
                                                                               gtgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTTGGGGAAAGGATAACTACCGGTTTTGGCGGCACCTTCTGGGACGGAcc                                                                 >  1:898552/1‑150 (MQ=255)
                                                                               gtgtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAcc                                                                 <  2:405165/150‑1 (MQ=255)
                                                                                         gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGa                                                       <  1:1391350/150‑1 (MQ=255)
                                                                                         gCCACAAGTATAAGCGGCAAATCCGGGTCACGTCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCAAATTCTGGTACGGACCAACGTCATGa                                                       <  1:1190281/150‑1 (MQ=255)
                                                                                          ccACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCACGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTTCGGACCAACGTCATGAg                                                      >  2:1560367/1‑150 (MQ=255)
                                                                                            ccaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTc                                                    <  1:1095198/149‑1 (MQ=255)
                                                                                                 tataAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCAc                                               <  2:1416170/150‑1 (MQ=255)
                                                                                                     aGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTgg                                           <  1:1020443/150‑1 (MQ=255)
                                                                                                      gCGGCAAATCCGGGTTACGCCGGGAGAGCAGACGGGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTCCCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTggg                                          <  1:832036/150‑1 (MQ=255)
                                                                                                      gCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTggg                                          <  1:611511/150‑1 (MQ=255)
                                                                                                          cAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGAATGCCCGCGCGAGAAATCATTACGGCAGGACTTAAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAAt                                      <  1:1518005/150‑1 (MQ=255)
                                                                                                            aaTCCGGGTTACGACGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                    <  1:689495/150‑1 (MQ=255)
                                                                                                                cGGGTTACGCCGGGAGAGCAGAAGCACGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTTCGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                <  1:1560367/150‑1 (MQ=255)
                                                                                                                cGGGGTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                <  2:138084/150‑1 (MQ=255)
                                                                                                                 gtgTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCg                               <  1:1324938/148‑1 (MQ=255)
                                                                                                                 ggggTATGCCGGGAGAGCAGAAGCGCGGG‑CAAGATGGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAATGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCg                               <  2:1150256/150‑1 (MQ=255)
                                                                                                                  gggtACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACATTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                              <  1:466976/147‑1 (MQ=255)
                                                                                                                  ggTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                              <  1:1161824/150‑1 (MQ=255)
                                                                                                                  ggTTACGCCGGGAGAGCAGAAGCGCGGG‑AAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCgg                              <  2:438869/150‑1 (MQ=255)
                                                                                                                        gCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTCCGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGt                        <  2:1509804/150‑1 (MQ=255)
                                                                                                                           gggAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCa                     <  2:1283644/150‑1 (MQ=255)
                                                                                                                            ggAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                    <  2:1618888/150‑1 (MQ=255)
                                                                                                                             gagaGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTCGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATg                   <  1:1139679/150‑1 (MQ=255)
                                                                                                                              cgagCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTGTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGc                  <  2:352648/149‑1 (MQ=255)
                                                                                                                               gagCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGcc                 <  1:431388/150‑1 (MQ=255)
                                                                                                                                  cAGAAGCGCGGG‑CAAGATTGCCCGCGCGCGAAATAAATAAAGCACGACATTAGCGGAAAAGATAACTAACGGTTTTGACGGCACAATATTGTTAGGACCCACGGCAAGAGTCCGGACCAGGGGAATCTTTTAGGCAACGTCCATGCCttt              >  2:1094169/1‑150 (MQ=255)
                                                                                                                                  cAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCttt              <  2:119260/150‑1 (MQ=255)
                                                                                                                                    gAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca            <  1:773185/150‑1 (MQ=255)
                                                                                                                                    ctcgcgcGGG‑CAATCTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca            >  2:840859/4‑150 (MQ=255)
                                                                                                                                      cgcgcGGG‑CAATCGTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGCCGGCACATTCTGGTACGGCCCAACGGCAGGAGTTGGCACCAGGGAAATCGTATCGGCCCAGACCATGCCTTTcccc          >  2:736157/2‑150 (MQ=255)
                                                                                                                                        cgcgGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                        >  2:309111/1‑102 (MQ=255)
                                                                                                                                        cgcgGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATg                                                        <  1:309111/102‑1 (MQ=255)
                                                                                                                                         gcgGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACt       >  1:217670/1‑150 (MQ=255)
                                                                                                                                         gcgGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTGTCACCACt       >  1:1090731/1‑150 (MQ=255)
                                                                                                                                          cgGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGACACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCCGCAACGTCCATGCCTTTCACCACtt      >  1:1132360/1‑150 (MQ=255)
                                                                                                                                            gg‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa    >  1:197174/1‑150 (MQ=255)
                                                                                                                                            gg‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa    >  1:1019163/1‑150 (MQ=255)
                                                                                                                                             g‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc   >  2:1257487/1‑150 (MQ=255)
                                                                                                                                              gctaGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGGGTTGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc   >  2:1410000/4‑150 (MQ=255)
                                                                                                                                               cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc  >  1:1661155/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAATTACGATTACGCTGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGG‑CAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACC  >  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: ‑