Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,040,850 T→A 10.7% 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,8500TA10.7% 215.9 / 16.4 75intergenic (+19/+15)PP_5292/PP_5293catabolite repression control protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (3/5);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

TTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAACA  >  NC_002947/6040706‑6040998
                                                                                                                                                |                                                                                                                                                    
ttCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAg                                                                                                                                                 <  1:34826/150‑1 (MQ=255)
 tggtgtACTAGATCCTCACCCCGCGCCTTCGCCGCTTCGTATTCAACGCTTGCTGGCCGCGTCAGCCACGCATCTCTCAACATTCATCGCTGATCTTCGACTATGGCTGGACGCTAACCCACTAGAAGGCCATTGGGGCTGCATTGCAGc                                                                                                                                                <  2:90000/145‑1 (MQ=255)
 tCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGc                                                                                                                                                >  1:105415/1‑150 (MQ=255)
  cGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGcc                                                                                                                                               >  1:267680/1‑150 (MQ=255)
   gACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGccc                                                                                                                                              <  2:291223/150‑1 (MQ=255)
    aCTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGcccc                                                                                                                                             >  1:325839/1‑150 (MQ=255)
      tACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCaa                                                                                                                                           >  1:1116193/1‑150 (MQ=255)
      tACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCaa                                                                                                                                           >  1:979915/1‑150 (MQ=255)
       aCCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAAt                                                                                                                                          <  1:1199874/150‑1 (MQ=255)
       aCCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATCCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAAt                                                                                                                                          <  1:10820/150‑1 (MQ=255)
        ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGATTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa                                                                                                                                                  <  1:647688/141‑1 (MQ=255)
        ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGATTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa                                                                                                                                                  >  2:647688/1‑141 (MQ=255)
        ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATg                                                                                                                                         <  2:150670/150‑1 (MQ=255)
        ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATg                                                                                                                                         <  2:267680/150‑1 (MQ=255)
              ccTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAtt                                                                                                                                   <  2:1242458/150‑1 (MQ=255)
              ccTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGCGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAttt                                                                                                                                  <  2:171292/150‑1 (MQ=255)
                tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTg                                                                                                                                 <  2:1055445/150‑1 (MQ=255)
                tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTg                                                                                                                                 >  1:993751/1‑150 (MQ=255)
                tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                         <  2:818847/142‑1 (MQ=255)
                tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                         >  1:818847/1‑142 (MQ=255)
                 cACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGa                                                                                                                                <  2:1203002/150‑1 (MQ=255)
                  aCCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGAt                                                                                                                               >  2:475304/1‑150 (MQ=255)
                    cccGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCg                                                                                                                             >  1:137905/1‑150 (MQ=255)
                              cgccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGat                                                                                                                   <  2:391647/150‑1 (MQ=255)
                                    ttCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCAGCCAGATCATCgg                                                                                                             >  1:730981/1‑150 (MQ=255)
                                     tCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCggg                                                                                                            <  2:1116193/150‑1 (MQ=255)
                                     tCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCggg                                                                                                            <  2:966615/150‑1 (MQ=255)
                                       gTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTa                                                                                                          >  1:595249/1‑150 (MQ=255)
                                       gTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTa                                                                                                          >  1:346329/1‑150 (MQ=255)
                                           gaaaCGCTCGCCTGCCGCGACAGCCACGTTTATCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGa                                                                                                      >  2:149864/3‑150 (MQ=255)
                                              aCGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGAACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATaa                                                                                                   <  1:7268/150‑1 (MQ=255)
                                                  tCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCt                                                                                               <  1:1241076/150‑1 (MQ=255)
                                                   cGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCtt                                                                                              <  2:993751/150‑1 (MQ=255)
                                                                   cACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGc                                                                              >  1:1137454/1‑150 (MQ=255)
                                                                    aCGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCa                                                                             >  2:352334/1‑150 (MQ=255)
                                                                      gCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCAcc                                                                           >  2:417799/1‑150 (MQ=255)
                                                                      gCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                         <  1:1017291/88‑1 (MQ=255)
                                                                      gCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                         >  2:1017291/1‑88 (MQ=255)
                                                                                 cATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCgatgat                                                                <  2:979915/150‑1 (MQ=255)
                                                                                        cGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAat                                                         >  1:558130/1‑150 (MQ=255)
                                                                                        cGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAat                                                         >  1:1200876/1‑150 (MQ=255)
                                                                                          cTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAatat                                                       >  2:375013/1‑150 (MQ=255)
                                                                                              tcgtcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATatcatc                                                   <  2:949227/150‑1 (MQ=255)
                                                                                               cgtcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATAtcatca                                                  <  2:137905/150‑1 (MQ=255)
                                                                                                 tcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATatcat                                                    <  2:830062/146‑1 (MQ=255)
                                                                                                 tcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATatcat                                                    >  1:830062/1‑146 (MQ=255)
                                                                                                   gACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCACCAATGGTTATTTGATCGGGCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCCCATTCGTCCGCACCGCCGCGATGATGATCAATATCATCACTGc                                              <  2:290410/150‑1 (MQ=255)
                                                                                                         gACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGt                                        <  1:149864/150‑1 (MQ=255)
                                                                                                          aCTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCggg                                                                                                            >  1:1174122/1‑81 (MQ=255)
                                                                                                          aCTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCggg                                                                                                            <  2:1174122/81‑1 (MQ=255)
                                                                                                           cTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGa                                      <  2:696062/150‑1 (MQ=255)
                                                                                                           cTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGa                                      <  2:1137454/150‑1 (MQ=255)
                                                                                                            tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatcatc                                                                                                               >  1:789664/1‑76 (MQ=255)
                                                                                                            tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatcatc                                                                                                               >  1:789471/1‑76 (MQ=255)
                                                                                                            tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatcatc                                                                                                               <  2:789664/76‑1 (MQ=255)
                                                                                                            tGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatcatc                                                                                                               <  2:789471/76‑1 (MQ=255)
                                                                                                             ggACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGc                                    >  1:1007705/1‑150 (MQ=255)
                                                                                                                cGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATg                                 <  2:346329/150‑1 (MQ=255)
                                                                                                                 gCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATgg                                <  2:558130/150‑1 (MQ=255)
                                                                                                                 gCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATgg                                <  2:1200876/150‑1 (MQ=255)
                                                                                                                   tGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTACGGGTGTTGATGTCCGTGAAGCGCCCGCACCGCCGCGCTACTGATCACCATCAACACCTTCACCGTGAGCATGGCa                              >  1:665558/1‑150 (MQ=255)
                                                                                                                     aCCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCACTCATTCGCACGAACCGAGCAGCTGATGTTCATAATCATGACGTCCACCGTGAGCATGGCaaa                            >  2:140803/1‑150 (MQ=255)
                                                                                                                         tCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGc                        <  1:724710/150‑1 (MQ=255)
                                                                                                                            aaaaGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCaa                                                          >  2:119723/1‑113 (MQ=255)
                                                                                                                            aaaaGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCaa                                                          <  1:119723/113‑1 (MQ=255)
                                                                                                                                    aTTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAACCAACACCATCACGAt             >  2:713860/1‑150 (MQ=255)
                                                                                                                                    aTTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGAt             >  2:722132/1‑150 (MQ=255)
                                                                                                                                          gCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCa                                                                                       <  1:945064/70‑1 (MQ=255)
                                                                                                                                          gCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCa                                                                                       >  2:945064/1‑70 (MQ=255)
                                                                                                                                          gCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCa                                   <  2:1179114/122‑1 (MQ=255)
                                                                                                                                          gCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCa                                   >  1:1179114/1‑122 (MQ=255)
                                                                                                                                            tGCATTGCAGACCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATACGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGa     <  2:105415/150‑1 (MQ=255)
                                                                                                                                              cATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAAc   <  1:352334/150‑1 (MQ=255)
                                                                                                                                               aTTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAACa  >  2:916802/1‑150 (MQ=255)
                                                                                                                                               aTTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAACa  >  2:545868/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAACA  >  NC_002947/6040706‑6040998

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