Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,150,143 C→A 6.3% intergenic (+98/‑299) PP_3654 → / → PP_3655 leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,150,1430CA 6.3% 191.8 / 2.8 63intergenic (+98/‑299)PP_3654/PP_3655leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base A (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

CTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCTTTT  >  NC_002947/4149994‑4150288
                                                                                                                                                     |                                                                                                                                                 
cTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAggccggcc                                                                                                                                                   >  2:940530/1‑150 (MQ=255)
 tCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCc                                                                                                                                                  >  2:345807‑M2/1‑149 (MQ=255)
  cGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCGAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCca                                                                                                                                                 >  2:529051‑M2/1‑148 (MQ=255)
   gACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAg                                                                                                                                                <  1:244755/150‑1 (MQ=255)
   gACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAg                                                                                                                                                <  1:1176513/150‑1 (MQ=255)
     ccTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCTGGTTGACCCGCGAAGAGGCCGGCCCAGGc                                                                                                                                              >  2:298112/1‑150 (MQ=255)
             atgaAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTc                                                                                                                                      >  2:1130045/1‑150 (MQ=255)
                   gaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTc                                                                                                                                      <  1:1034425/144‑1 (MQ=255)
                   gaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTc                                                                                                                                      >  2:1034425/1‑144 (MQ=255)
                     agTGATCCAGGATCGGTTGCTCCCCCTGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCAGGAAAACCCGCCACTGCAAAAGTCCCGTGGTCGGTCGCGCGAGCGGGGTGACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCtcctc                                                                                                                              >  2:1041653/1‑146 (MQ=255)
                     agTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGccc                                                                                                                              <  1:940530/150‑1 (MQ=255)
                      gTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCt                                                                                                                             >  1:702046/1‑150 (MQ=255)
                                aTCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACgcgc                                                                                                                   >  2:945008/1‑150 (MQ=255)
                                    gTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACAt                                                                                                               >  1:694347/1‑150 (MQ=255)
                                     ttGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATc                                                                                                              <  1:1130045/150‑1 (MQ=255)
                                       gCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGc                                                                                                            <  1:945008/150‑1 (MQ=255)
                                        cTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCa                                                                                                           >  1:575723/1‑150 (MQ=255)
                                         tCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCAc                                                                                                          >  2:966821/1‑150 (MQ=255)
                                           cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGc                                                                                                                                >  1:229542/1‑126 (MQ=255)
                                           cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGc                                                                                                                                <  2:229542/126‑1 (MQ=255)
                                                            tctTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCcaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaagg                                                                                       <  2:491476‑M2/150‑61 (MQ=255)
                                                            tctTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGAt                                                                                       >  1:1181113/1‑150 (MQ=255)
                                                                 tgCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGt                                                                                  >  1:402132/1‑150 (MQ=255)
                                                                  gCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGtt                                                                                 <  1:966821/150‑1 (MQ=255)
                                                                    cGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGc                                                                               <  2:702046/150‑1 (MQ=255)
                                                                      gCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCa                                                                             >  2:1129857/1‑150 (MQ=255)
                                                                         ccTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCAccc                                                                          >  1:955236/1‑150 (MQ=255)
                                                                                 tgTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCgccgcc                                                                  >  1:351409/2‑150 (MQ=255)
                                                                                  gTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc                                                                 >  1:462381/1‑150 (MQ=255)
                                                                                    aaaCCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCcaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacattgg                                                               >  1:491476‑M2/1‑66 (MQ=255)
                                                                                    aaaCCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCct                                                               <  2:462381/150‑1 (MQ=255)
                                                                                       cccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCcaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagc                                                                                 <  2:224836‑M2/129‑67 (MQ=255)
                                                                                       cccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCcaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagc                                                                                 >  1:224836‑M2/1‑63 (MQ=255)
                                                                                        ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCg                                                                                                             >  2:286235/1‑100 (MQ=255)
                                                                                        ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCg                                                                                                             <  1:286235/100‑1 (MQ=255)
                                                                                        ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAc                                                           <  2:1181113/150‑1 (MQ=255)
                                                                                        ccGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAc                                                           >  2:1017293/1‑150 (MQ=255)
                                                                                        ccGCGCCTGCAAAGGTGACGTAGCCCGTCGAAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAc                                                           >  2:850433/1‑150 (MQ=255)
                                                                                         cgcgCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGGGGCCGGCAcaagaggaaacagctacatggggagcaaccgatcctggatcacttc                                                                                                     <  2:983948‑M2/107‑47 (MQ=255)
                                                                                         cgcgCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGGGGCCGGCAcaagaggaaacagctacatggggagcaaccgatcctggatcacttc                                                                                                     >  1:983948‑M2/1‑61 (MQ=255)
                                                                                         cgcgCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCACTCGTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAcc                                                          <  1:980573/150‑1 (MQ=255)
                                                                                         cgcgCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacattgggaata                                                          <  1:1126973‑M2/150‑90 (MQ=255)
                                                                                           cgcCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGACAAGCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTAcccc                                                        >  1:917955/1‑150 (MQ=255)
                                                                                                  aaaGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCt                                                 >  1:1176103/1‑150 (MQ=255)
                                                                                                     ggTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTca                                              <  2:694347/150‑1 (MQ=255)
                                                                                                           gTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAcaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacattgggaatactggtcaggtgctcgtca                                        >  2:1126973‑M2/1‑43 (MQ=255)
                                                                                                            tAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGGTGACTCCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCCGCt                                       >  1:1157071/1‑150 (MQ=255)
                                                                                                             aGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTg                                      <  1:1041653/150‑1 (MQ=255)
                                                                                                             aGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTg                                      <  2:402132/150‑1 (MQ=255)
                                                                                                              gCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGGCTCGTCCCTCTCCCCTCTACCCCAAGGCCTTCAGCCAGCTgg                                     >  2:886507/1‑150 (MQ=39)
                                                                                                                   tCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCActc                                <  2:575723/150‑1 (MQ=255)
                                                                                                                       aGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCCGCTGGCACTCTCCt                            >  2:4820/1‑150 (MQ=255)
                                                                                                                        ggAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCtt                           <  1:1129857/150‑1 (MQ=255)
                                                                                                                           gCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc                                                                 <  1:960218/109‑1 (MQ=255)
                                                                                                                           gCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGccc                                                                 >  2:960218/1‑109 (MQ=255)
                                                                                                                            cGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTa                       <  2:955236/150‑1 (MQ=255)
                                                                                                                              ggTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTAcc                     >  2:472233/1‑150 (MQ=255)
                                                                                                                               gTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGcc                                                  <  1:709207/120‑1 (MQ=255)
                                                                                                                               gTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGcc                                                  >  2:709207/1‑120 (MQ=255)
                                                                                                                                   aCCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCAc                <  1:1017293/150‑1 (MQ=255)
                                                                                                                                     ccGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCAccc              >  1:226097/1‑150 (MQ=255)
                                                                                                                                     ccGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCAccc              >  2:721823/1‑150 (MQ=255)
                                                                                                                                         gAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCtt          <  2:226097/150‑1 (MQ=255)
                                                                                                                                            gagGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAACCGTCACTCTCCTTGCTACCGTCACCCGCTTAAg       <  1:886507/150‑1 (MQ=255)
                                                                                                                                                 cggccCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCtttt  <  1:298112/150‑1 (MQ=255)
                                                                                                                                                   gccCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCtt    <  1:978258/146‑1 (MQ=255)
                                                                                                                                                   gccCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCtt    >  2:978258/1‑146 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCTTTT  >  NC_002947/4149994‑4150288

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: