Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,040,850 T→A 5.4% intergenic (+19/+15) PP_5292 → / ← PP_5293 catabolite repression control protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,040,8500TA 5.4% 240.5 / 4.4 74intergenic (+19/+15)PP_5292/PP_5293catabolite repression control protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base A (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

ATCGGTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAAC  >  NC_002947/6040701‑6040997
                                                                                                                                                     |                                                                                                                                                    
aTCGGTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAt                                                                                                                                                      <  1:394619/150‑1 (MQ=255)
     ttCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCATGCCGCGTCAGACACGCTTCTCCCAGCATGCACCGCTGATCGCCAACTATGACTGGACGCTAACCCTTTAAAAGGCCATTGGGGCTTCATTGCAg                                                                                                                                                 >  1:342250/1‑150 (MQ=255)
           tACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAg                                                                                                                                                 <  1:1193000/144‑1 (MQ=255)
           tACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGATCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAg                                                                                                                                                 >  2:1193000/1‑144 (MQ=255)
           tACCAGATCCTCACCCCGGACCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCAACACAGCCACGCTTCTCCCAGCTTGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCaa                                                                                                                                           >  1:199659/1‑150 (MQ=255)
                gATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      >  1:1578888/1‑150 (MQ=255)
                   ccTCACCCGGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAtt                                                                                                                                   <  1:703752/150‑1 (MQ=255)
                    cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAttt                                                                                                                                  >  1:506584/1‑150 (MQ=255)
                      cACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGa                                                                                                                                >  1:729243/1‑150 (MQ=255)
                         cccGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCg                                                                                                                             <  2:606784/150‑1 (MQ=255)
                            gggCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGgcc                                                                                                                          <  2:1578888/150‑1 (MQ=255)
                              gCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGgccgc                                                                                                                        >  2:338640/1‑150 (MQ=255)
                                cTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAtt                                                                                                                                                     >  1:1244729/1‑119 (MQ=255)
                                cTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAtt                                                                                                                                                     <  2:1244729/119‑1 (MQ=255)
                                cTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc                                                                                                                                              >  2:226652/1‑126 (MQ=255)
                                cTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc                                                                                                                                              <  1:226652/126‑1 (MQ=255)
                                   cgccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCGATggcct                                                                                                                                     >  2:1304655/1‑129 (MQ=255)
                                   cgccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCGATggcct                                                                                                                                     <  1:1304655/132‑4 (MQ=255)
                                    gccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatc                                                                                                                  >  1:573817/1‑150 (MQ=255)
                                             tGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAg                                                                                                         >  2:1471205/1‑150 (MQ=255)
                                                gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGa                                                                                                      >  2:620271/1‑150 (MQ=255)
                                                  aaCGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATa                                                                                                    >  2:1091936/1‑150 (MQ=255)
                                                         gCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTg                                                                                             >  1:1513539/1‑150 (MQ=255)
                                                         gCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGATTAGCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTg                                                                                             <  1:15071/150‑1 (MQ=255)
                                                          ccTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGc                                                                                            <  1:338640/150‑1 (MQ=255)
                                                            tGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGccc                                                                                          <  2:493164/150‑1 (MQ=255)
                                                               cgcgTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCAAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCa                                                                                       <  2:99896/150‑1 (MQ=255)
                                                                      gCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCg                                                                                                              >  2:27825/1‑120 (MQ=255)
                                                                      gCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCg                                                                                                              <  1:27825/120‑1 (MQ=255)
                                                                      gCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGccc                                                                               <  2:573817/150‑1 (MQ=255)
                                                                         aCGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAt                                                                                                                                    <  1:1550870/95‑1 (MQ=255)
                                                                         aCGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAt                                                                                                                                    >  2:1550870/1‑95 (MQ=255)
                                                                           gCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCAcc                                                                          >  2:796248/1‑150 (MQ=255)
                                                                                  ccAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAAt                                                                                                                                          >  1:418871/1‑80 (MQ=255)
                                                                                  ccAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAAt                                                                                                                                          <  2:418871/80‑1 (MQ=255)
                                                                                         gCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCgatgatgat                                                            >  2:167989/1‑150 (MQ=255)
                                                                                          cACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAtt                                                                                                                                   <  1:1286516/79‑1 (MQ=255)
                                                                                          cACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAtt                                                                                                                                   >  2:1286516/1‑79 (MQ=255)
                                                                                             cGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCa                                                                                       <  1:339148/120‑1 (MQ=255)
                                                                                             cGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCa                                                                                       >  2:339148/1‑120 (MQ=255)
                                                                                             cGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAat                                                        >  1:342150/1‑150 (MQ=255)
                                                                                                      tcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTCGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAatat                                                      <  1:81717/143‑1 (MQ=255)
                                                                                                      tcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTCGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAatat                                                      >  2:81717/1‑143 (MQ=255)
                                                                                                               aCTGGACGCTGACAATATAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTg                                      >  2:1381199/1‑150 (MQ=255)
                                                                                                                cTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCAcc                                         <  2:461291/146‑1 (MQ=255)
                                                                                                                cTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCAcc                                         >  1:461291/1‑146 (MQ=255)
                                                                                                                cTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGa                                     <  1:796248/150‑1 (MQ=255)
                                                                                                                 tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCAt                                                                                     <  1:1462111/101‑1 (MQ=255)
                                                                                                                 tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCAt                                                                                     >  2:1462111/1‑101 (MQ=255)
                                                                                                                 tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAg                                    >  2:1281971/1‑150 (MQ=255)
                                                                                                                 tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAg                                    >  2:849738/1‑150 (MQ=255)
                                                                                                                   gACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTg                                                                                             <  1:1537379/92‑1 (MQ=255)
                                                                                                                   gACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTg                                                                                             >  2:1537379/1‑92 (MQ=255)
                                                                                                                   gACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGcc                                                                                           <  1:258073/94‑1 (MQ=255)
                                                                                                                   gACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGcc                                                                                           >  2:258073/1‑94 (MQ=255)
                                                                                                                   gACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCa                                  <  1:495780/150‑1 (MQ=255)
                                                                                                                      gCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATgg                               <  1:167989/150‑1 (MQ=255)
                                                                                                                      gCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATgg                               >  1:1483391/1‑150 (MQ=255)
                                                                                                                        tGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGGTCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCa                             >  2:1629845/1‑150 (MQ=255)
                                                                                                                           gatcttAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAAc                          <  1:184099/145‑1 (MQ=255)
                                                                                                                           ccATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAAc                          <  1:1281971/150‑1 (MQ=255)
                                                                                                                                 aaaaGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCAcc                    >  1:242509/1‑150 (MQ=255)
                                                                                                                                 aaaaGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATAGGGTAGCGATAAGGCTTACTCACATTCGCCCGCATCACCTGGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACaaaatc                    >  2:184099/1‑145 (MQ=255)
                                                                                                                                  aaaGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGAt                                                                                                                               >  1:599363/1‑43 (MQ=255)
                                                                                                                                  aaaGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGAt                                                                                                                               <  2:599363/43‑1 (MQ=255)
                                                                                                                                    aGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGc                 >  1:778449/1‑150 (MQ=255)
                                                                                                                                         aTTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCAcc                                         <  1:1364293/121‑1 (MQ=255)
                                                                                                                                         aTTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCAcc                                         >  2:1364293/1‑121 (MQ=255)
                                                                                                                                             gggCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCgatga                                                                <  1:346652/94‑1 (MQ=255)
                                                                                                                                             gggCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCgatga                                                                >  2:346652/1‑94 (MQ=255)
                                                                                                                                             gggCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAac        >  1:1357808/1‑150 (MQ=255)
                                                                                                                                              ggCTGCATTGTAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAaca       <  2:729243/150‑1 (MQ=255)
                                                                                                                                              ggCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCaaa                           <  2:1514091/130‑1 (MQ=255)
                                                                                                                                              ggCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCaaa                           >  1:1514091/1‑130 (MQ=255)
                                                                                                                                               gCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAacac      >  1:327911/1‑150 (MQ=255)
                                                                                                                                                  gCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCATCACGTCGCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGaa   <  1:1103085/150‑1 (MQ=255)
                                                                                                                                                   cATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAAc  <  2:1357808/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGGTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAAC  >  NC_002947/6040701‑6040997

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