Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,951 A→G 32.0% intergenic (+40/‑61) quiA → / → oprB‑III quinate dehydrogenase/carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,049,9510AG32.0% 87.8 / 64.1 79intergenic (+40/‑61)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base G (15/10);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-03

GACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGC  >  NC_002947/4049803‑4050098
                                                                                                                                                    |                                                                                                                                                   
gACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                                                                    <  2:215078/150‑1 (MQ=255)
   cTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCACTCCTACAAGGcc                                                                                                                                                 >  1:398456/1‑150 (MQ=255)
    tACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGccc                                                                                                                                                <  2:2148357/150‑1 (MQ=255)
      cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCt                                                                                                                                              <  2:1724063/150‑1 (MQ=255)
          tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTCCAAGCCCCCTGTAg                                                                                                                                          <  2:1203038/150‑1 (MQ=255)
          tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAg                                                                                                                                          >  2:917435/1‑150 (MQ=255)
          tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                          >  2:2085833/1‑150 (MQ=255)
          tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                          <  1:1712644/150‑1 (MQ=255)
            gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGGCTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGa                                                                                                                                        >  1:528575/1‑150 (MQ=255)
            gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTACGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                        <  1:347577/150‑1 (MQ=39)
              cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                      >  2:1488340/1‑150 (MQ=255)
               gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                      >  1:960321/1‑149 (MQ=255)
               gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                      <  2:960321/149‑1 (MQ=255)
               gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGCGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCg                                                                                                                                     >  2:607714/1‑150 (MQ=255)
                aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGGCTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                    <  2:528575/150‑1 (MQ=255)
                aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                    <  1:1495871/150‑1 (MQ=255)
                aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                    <  1:96849/150‑1 (MQ=255)
                 aGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg                                                                                                                                   <  1:2085833/150‑1 (MQ=255)
                  gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGc                                                                                                                                  <  1:414601/150‑1 (MQ=255)
                    ccGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCtt                                                                                                                                >  2:673185/1‑150 (MQ=255)
                          aGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                    >  1:726014/1‑140 (MQ=255)
                          aGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                    <  2:726014/140‑1 (MQ=255)
                                aCGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGaagca                                                                                                                    <  1:917435/150‑4 (MQ=37)
                                aCGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                    <  1:1488340/150‑1 (MQ=255)
                                  gtggtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGccc                                                                                                                  <  2:398456/150‑1 (MQ=255)
                                    ggtgATCACTTCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACCAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCt                                                                                                                >  2:515842/1‑150 (MQ=255)
                                     gtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg                                                                                                                                   <  1:1040050/130‑1 (MQ=255)
                                     gtgATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg                                                                                                                                   >  2:1040050/1‑130 (MQ=255)
                                               ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACTCGGCACCTGTAGGAGCGGGCTCGCCCGCGAATGCCCCTCAAAGACAAcc                                                                                                     <  1:538844/150‑1 (MQ=255)
                                                  gcggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCga                                                                                                  <  2:284073/150‑1 (MQ=255)
                                                      cgcccgccAGTCGACTGACCGCGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGAtt                                                                                              <  1:607714/150‑1 (MQ=255)
                                                       gcccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTc                                                                                             <  2:1522784/150‑1 (MQ=255)
                                                         ccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTg                                                                                           <  2:770902/150‑1 (MQ=255)
                                                             cAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAAt                                                                                       <  1:1461334/150‑1 (MQ=255)
                                                                tCGACTGTCCGTGGCGACTATGTGATTTCTTACGCTCTGCCATAAAACCGCGTCGGCTGCTTCGCGGGCAAGCCCGCCCCTAAAACGCCCCTGTAGAAGCGGGCTTGACCGCGAATGCCCCCCAACGACAACCCGAGATACTGTACTGcc                                                                                    >  2:215474/1‑150 (MQ=255)
                                                                    cTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                    <  2:790180/114‑1 (MQ=255)
                                                                    cTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                    >  1:790180/1‑114 (MQ=255)
                                                                     tGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGc                                                                                                                                  <  2:1053051/99‑1 (MQ=25)
                                                                     tGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGc                                                                                                                                  >  1:1053051/1‑99 (MQ=25)
                                                                     tGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCg                                                                               >  1:1409396/1‑150 (MQ=255)
                                                                        ccGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                          <  1:1686963/88‑1 (MQ=255)
                                                                        ccGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                          >  2:1686963/1‑88 (MQ=255)
                                                                           tGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTc                                                                         <  1:673185/150‑1 (MQ=255)
                                                                                   aTGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                          <  1:152921/77‑1 (MQ=255)
                                                                                   aTGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                          >  2:152921/1‑77 (MQ=255)
                                                                                       gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACaa                                                                                                       >  1:257108/1‑108 (MQ=255)
                                                                                       gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACaa                                                                                                       <  2:257108/108‑1 (MQ=255)
                                                                                       gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCa                                                             >  1:617751/1‑150 (MQ=255)
                                                                                        aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCAc                                                            >  1:2094502/1‑150 (MQ=255)
                                                                                        aTTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCAc                                                            >  1:463822/1‑150 (MQ=255)
                                                                                         tttCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCAcc                                                           <  2:463822/150‑1 (MQ=255)
                                                                                              tACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgct                                                      >  1:446478/1‑150 (MQ=255)
                                                                                                 gCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgctgct                                                   >  1:568577/1‑150 (MQ=255)
                                                                                                      gCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAg                                                                                                                                          <  2:923447/58‑1 (MQ=17)
                                                                                                      gCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAg                                                                                                                                          >  1:923447/1‑58 (MQ=17)
                                                                                                       ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccct                                             <  2:617751/150‑1 (MQ=255)
                                                                                                        cGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctg                                            <  2:568577/150‑1 (MQ=255)
                                                                                                          tAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctggc                                          >  1:1531732/1‑150 (MQ=255)
                                                                                                                         tGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCt                           >  2:726444/1‑150 (MQ=255)
                                                                                                                                  ggCAAGCCCGCTCCTACCAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                  <  1:515842/150‑1 (MQ=255)
                                                                                                                                        cccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct            >  2:363952/1‑150 (MQ=255)
                                                                                                                                        cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct            >  2:2104287/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc           >  1:1624576/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACACTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtcag           >  1:2006856/1‑148 (MQ=38)
                                                                                                                                          cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGc                              <  1:2008803/130‑1 (MQ=255)
                                                                                                                                          cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGc                              >  2:2008803/1‑130 (MQ=255)
                                                                                                                                           gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt         >  1:178128/1‑150 (MQ=255)
                                                                                                                                           gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCACCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt         >  1:278448/1‑150 (MQ=255)
                                                                                                                                           gCTCCTACAGGGGACTTGTAGGAGCGGGCTGGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAACTCACACCACCGCCACCCTGCTGCGGGCCCTGGCCAGCACCACTCCCCTTGCCAGCGAACATGTCAt         >  2:315986/1‑150 (MQ=255)
                                                                                                                                           gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGCCGATCTCAt         >  1:719753/1‑150 (MQ=255)
                                                                                                                                           gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt         >  1:418390/1‑150 (MQ=255)
                                                                                                                                           gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGcagcagg                                <  1:752547/127‑6 (MQ=255)
                                                                                                                                           gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGcagcagg                                >  2:752547/1‑122 (MQ=255)
                                                                                                                                             tCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                  >  1:2292681/1‑139 (MQ=255)
                                                                                                                                             tCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                  <  2:2292681/139‑1 (MQ=255)
                                                                                                                                             tCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa       >  2:2287637/1‑150 (MQ=255)
                                                                                                                                               ttACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCTTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAcc     <  1:992203/149‑1 (MQ=255)
                                                                                                                                                 aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGAAGATCTGATGAgcca   >  2:992203/1‑149 (MQ=255)
                                                                                                                                                  cAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgc  <  2:446478/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACCTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGC  >  NC_002947/4049803‑4050098

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