Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,098 A→G 9.3% intergenic (‑31/+22) PP_5445 ← / ← PP_t71 hypothetical protein/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947467,0980AG 9.3% 206.6 / 9.0 75intergenic (‑31/+22)PP_5445/PP_t71hypothetical protein/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base G (4/3);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GCGCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCT  >  NC_002947/466952‑467243
                                                                                                                                                  |                                                                                                                                                 
gcgcgcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGTCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAAt                                                                                                                                                 >  7:209007/1‑149 (MQ=255)
gcgcgcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGAt                                                                                                                                                 >  5:285734/1‑149 (MQ=255)
gcgcgcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAAt                                                                                                                                                 >  1:327840/1‑149 (MQ=255)
   cgcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATtgc                                                                                                                                              >  8:159872/1‑149 (MQ=255)
     cTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcgc                                                                                                                                            <  1:217685/149‑1 (MQ=255)
            aCCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCttt                                                                                                                                      <  7:159872/148‑1 (MQ=255)
            aCCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCttt                                                                                                                                      <  6:285734/148‑1 (MQ=255)
            aCCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCttt                                                                                                                                      <  6:1645/148‑1 (MQ=255)
            aCCGGCTTCTTGGTGGGGCTGATTTCGATGAGATACCGTTTTATCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCttt                                                                                                                                      <  8:275255/148‑1 (MQ=255)
             ccGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGt                                                                                                                                    <  4:90313/149‑1 (MQ=255)
                   tcttGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATc                                                                                                                              >  1:7981/1‑149 (MQ=255)
                             gCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGcc                                                                                                                     >  7:137621/1‑148 (MQ=255)
                             gCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGcc                                                                                                                     >  5:27266/1‑148 (MQ=255)
                                 aTTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                <  7:50773/149‑1 (MQ=255)
                                  tttCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACg                                                                                                               <  6:56523/149‑1 (MQ=255)
                                      gCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACt                                                                                                           >  5:206511/1‑149 (MQ=255)
                                       cTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTc                                                                                                          >  6:160947/1‑149 (MQ=255)
                                          agaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                   <  5:72153/121‑1 (MQ=255)
                                          agaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                   >  6:72153/1‑121 (MQ=255)
                                             tACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCg                                                                                                    >  8:128991/1‑149 (MQ=255)
                                               ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGTCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                  <  8:209007/149‑1 (MQ=255)
                                               ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGAc                                                                                                            >  2:161082/1‑139 (MQ=255)
                                               ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGAc                                                                                                            <  1:161082/139‑1 (MQ=255)
                                               ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                  <  2:58345/149‑1 (MQ=255)
                                               ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                  <  7:188950/149‑1 (MQ=255)
                                                cGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCTCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTgg                                                                                                 <  1:146301/149‑1 (MQ=255)
                                                 gTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGa                                                                                                <  8:137621/149‑1 (MQ=255)
                                                 gTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGa                                                                                                <  4:165034/149‑1 (MQ=255)
                                                      aGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCa                                                                                             >  7:300297/1‑147 (MQ=255)
                                                         gcTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGa                                                                                        >  5:214486/1‑149 (MQ=255)
                                                             cTTCCAAGGTCCTGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCAGGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                    <  6:245469/149‑1 (MQ=255)
                                                             cTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGAATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                    <  5:218478/149‑1 (MQ=255)
                                                              ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                    >  6:222264/1‑148 (MQ=255)
                                                              ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                    >  8:180012/1‑148 (MQ=255)
                                                                ccAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGcc                                                                                                                     <  5:156170/113‑1 (MQ=255)
                                                                ccAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGcc                                                                                                                     >  6:156170/1‑113 (MQ=255)
                                                                ccAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGa                                                                                  <  4:318370/148‑1 (MQ=255)
                                                                 cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAg                                                                                          >  7:62026/1‑139 (MQ=255)
                                                                 cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAg                                                                                          <  8:62026/139‑1 (MQ=255)
                                                                   aGGTCCTGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGGGGACCAAAGGATTATGAGTcc                                                                              <  6:328442/149‑1 (MQ=255)
                                                                    ggTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCct                                                                             <  6:206511/149‑1 (MQ=255)
                                                                        aTGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTg                                                                          <  2:28064/148‑1 (MQ=255)
                                                                           cgcTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGctc                                                                       <  2:327840/148‑1 (MQ=255)
                                                                                     cgctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGACCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGa                                                            <  1:49088/149‑1 (MQ=255)
                                                                                          cgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                <  6:109911/92‑1 (MQ=255)
                                                                                          cgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                >  5:109911/1‑92 (MQ=255)
                                                                                           gctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACAAGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGGTCTAACCAACTgctcta                                                        >  3:234866/1‑142 (MQ=255)
                                                                                            ctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTggg                                                                                                                       <  1:290037/83‑1 (MQ=255)
                                                                                            ctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTggg                                                                                                                       >  2:290037/1‑83 (MQ=255)
                                                                                             tgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAgg                                                    >  8:270537/1‑149 (MQ=255)
                                                                                             tgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAgg                                                    <  8:300297/149‑1 (MQ=255)
                                                                                                   tACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTAGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTTAGCTATATGCCCTCa                                              >  5:197984/1‑149 (MQ=255)
                                                                                                         ccccTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGt                                        <  7:128991/149‑1 (MQ=255)
                                                                                                          cccTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTg                                       >  7:226685/1‑149 (MQ=255)
                                                                                                            cTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTggg                                     <  6:27266/149‑1 (MQ=255)
                                                                                                               aCCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCg                                   <  5:160947/148‑1 (MQ=255)
                                                                                                                  atttatttCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTATGATCTAACCAACTGAGATATAGGCCCTCAGTAAGTGGGCGGAt                                >  7:135591/1‑148 (MQ=255)
                                                                                                                   tgtatttCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATta                              >  1:296207/3‑149 (MQ=255)
                                                                                                                          cTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGa                        <  7:270537/148‑1 (MQ=255)
                                                                                                                          cTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGa                        <  5:222264/148‑1 (MQ=255)
                                                                                                                           ttGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCt                                                         >  2:225863/1‑114 (MQ=255)
                                                                                                                           ttGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCt                                                         <  1:225863/114‑1 (MQ=255)
                                                                                                                           ttGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAgg                      >  3:304954/1‑149 (MQ=255)
                                                                                                                               gggCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGc                                    >  8:141845/1‑131 (MQ=255)
                                                                                                                               gggCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGc                                    <  7:141845/131‑1 (MQ=255)
                                                                                                                                ggCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTc                 <  6:214486/149‑1 (MQ=255)
                                                                                                                                 gCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCg                <  7:180012/149‑1 (MQ=255)
                                                                                                                                   aTGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCa              >  2:168619/1‑149 (MQ=255)
                                                                                                                                     gacaaaGCCGCGCGATTGCGCGGCTTTTTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTaa                                          <  5:260449/116‑1 (MQ=39)
                                                                                                                                     gacaaaGCCGCGCGATTGCGCGGCTTTTTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTaa                                          >  6:260449/4‑119 (MQ=39)
                                                                                                                                      aaaaaGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCt           >  6:285019/1‑149 (MQ=255)
                                                                                                                                        aaaGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCtgt         >  6:145939/1‑149 (MQ=255)
                                                                                                                                        aaaGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCtgt         >  1:151889/1‑149 (MQ=255)
                                                                                                                                             cgcgCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCAt    >  3:43821/1‑149 (MQ=255)
                                                                                                                                               cgcAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCt  <  2:151889/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCT  >  NC_002947/466952‑467243

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