Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,101 A→G 7.0% 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,1010AG 7.0% 248.1 / 6.5 86intergenic (‑31/+22)PP_5445/PP_t71hypothetical protein/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (36/44);  new base G (4/2);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

CGGCGCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGAAA  >  NC_002947/466953‑467250
                                                                                                                                                    |                                                                                                                                                     
cGGCGCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCaa                                                                                                                                                      >  1:2436934/1‑150 (MQ=255)
    gcgcTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgc                                                                                                                                                  <  1:618437/150‑1 (MQ=255)
        tGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCgg                                                                                                                                              >  1:388206/1‑150 (MQ=255)
               ccGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGCCCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                       <  1:1093998/150‑1 (MQ=255)
               ccGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                       <  2:1387222/150‑1 (MQ=255)
               ccGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                       <  1:376125/150‑1 (MQ=255)
                cGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCGCGGCGCCGCTGCGCTGGCTGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCg                                                                                                                                      >  2:155129/1‑150 (MQ=255)
                 ggCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcgc                                                                                                                                                >  1:434316/1‑139 (MQ=255)
                 ggCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcgc                                                                                                                                                <  2:434316/139‑1 (MQ=255)
                    ttcttGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATc                                                                                                                                  >  2:1417554/1‑150 (MQ=255)
                        tggtggGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCtggt                                                                                                                              <  2:1945459/150‑1 (MQ=255)
                           gggggCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTggg                                                                                                                           <  1:16338/149‑1 (MQ=255)
                               gCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCca                                                                                                                       >  1:1381924/1‑150 (MQ=255)
                                    tttCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACgg                                                                                                                  >  1:1371733/1‑150 (MQ=255)
                                     ttCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGa                                                                                                                 <  1:1008039/150‑1 (MQ=255)
                                     ttCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGa                                                                                                                 <  2:2436934/150‑1 (MQ=255)
                                          tGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGa                                                                                                            >  2:2116584/1‑150 (MQ=255)
                                            agaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATc                                                                                                                                                    <  2:1710115/108‑1 (MQ=255)
                                            agaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATc                                                                                                                                                    >  1:1710115/1‑108 (MQ=255)
                                            agaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAAc                                                                                                          >  1:1831057/1‑150 (MQ=255)
                                            agaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAAc                                                                                                          <  2:388206/150‑1 (MQ=255)
                                             gaTACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAAcc                                                                                                         <  2:1381924/150‑1 (MQ=255)
                                                aCCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                      >  2:5731/1‑150 (MQ=255)
                                                aCCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                      >  1:301101/1‑150 (MQ=255)
                                                 ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                      >  1:657090/1‑149 (MQ=255)
                                                 ccGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTg                                                                                                      <  2:657090/149‑1 (MQ=255)
                                                  cGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGa                                                                                                    <  1:5731/150‑1 (MQ=255)
                                                  cGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGa                                                                                                    >  1:1312983/1‑150 (MQ=255)
                                                  cGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGa                                                                                                    <  1:1439508/150‑1 (MQ=255)
                                                       tAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                               <  1:1192828/150‑1 (MQ=255)
                                                         gcgcTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                               <  2:40397/148‑1 (MQ=255)
                                                         gcgcTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCaaa                                                                                               >  1:40397/1‑148 (MQ=255)
                                                           gataattCCAGGGGCATGCGATCGGCGCTGCTGCGCTAGATGTGCGGGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGAt                                                                                           <  2:1721270/145‑1 (MQ=255)
                                                            cTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGAtt                                                                                          <  1:1725537/150‑1 (MQ=255)
                                                            cTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGAtt                                                                                          >  2:1725537/1‑150 (MQ=255)
                                                              aCTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                        <  2:1831057/150‑1 (MQ=255)
                                                              aCTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCCGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                        <  2:1235570/150‑1 (MQ=255)
                                                               cTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                    <  1:1175926/121‑1 (MQ=255)
                                                               cTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                    >  2:1175926/1‑121 (MQ=255)
                                                                ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGTCCCACACGGACTCGAACCGTGGACCAAAGGATTATGa                                                                                      <  2:1117749/150‑1 (MQ=255)
                                                                ttCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGa                                                                                      <  2:1117686/150‑1 (MQ=255)
                                                                  ccAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGt                                                                                    >  1:932839/1‑150 (MQ=255)
                                                                   cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTc                                                                                   <  2:4477/150‑1 (MQ=255)
                                                                   cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTc                                                                                   <  2:5930/150‑1 (MQ=255)
                                                                   cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTc                                                                                   >  1:1460846/1‑150 (MQ=255)
                                                                               cTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAAc                                                                       <  2:699297/150‑1 (MQ=255)
                                                                                tCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAAcc                                                                      >  2:86267/1‑150 (MQ=255)
                                                                                     gccgcTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTg                                                                 >  2:1016358/1‑150 (MQ=255)
                                                                                      ccgcTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGa                                                                <  1:900291/150‑1 (MQ=255)
                                                                                       cgctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAAACGTGGCTCACAGGAAGAGGAGGCCTCGGCGCTAACCAACTGAg                                                               >  2:2092047/1‑150 (MQ=255)
                                                                                        gctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAAc                                                                       <  2:2135197/141‑1 (MQ=255)
                                                                                        gctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAAc                                                                       >  1:2135197/1‑141 (MQ=255)
                                                                                        gctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGc                                                              >  2:2115422/1‑150 (MQ=255)
                                                                                        gctgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGc                                                              >  2:2115578/1‑150 (MQ=255)
                                                                                          tgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCta                                                            >  1:338109/1‑150 (MQ=255)
                                                                                          tgcgctgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCta                                                            >  2:2238908/1‑150 (MQ=255)
                                                                                               tgGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGc                                                       <  1:2238908/150‑1 (MQ=255)
                                                                                                 gATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGccc                                                     >  1:946186/1‑150 (MQ=255)
                                                                                                    gTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCa                                                  <  1:1016358/150‑1 (MQ=255)
                                                                                                            cccTGACCATTTATTTCGGGCGGGCATGACGAAGCCGCTCAACCGCGCGGCTTTGTCGTATCTAGTGGACAGACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTgg                                          <  1:1615232/150‑1 (MQ=255)
                                                                                                              cTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCttt                                                                                                                                          <  1:453207/52‑1 (MQ=25)
                                                                                                              cTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCGATTGCGCGGCttt                                                                                                                                          >  2:453207/1‑52 (MQ=25)
                                                                                                              cTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                       <  2:1854605/55‑1 (MQ=255)
                                                                                                              cTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTc                                                                                                                                       >  1:1854605/1‑55 (MQ=255)
                                                                                                              cTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGc                                        <  2:1460846/150‑1 (MQ=255)
                                                                                                                 aCCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTgg                                          >  2:186581/1‑145 (MQ=255)
                                                                                                                 aCCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTgg                                          <  1:186581/145‑1 (MQ=255)
                                                                                                                     tttatttCTTGCGGCCATGAAAAAGCCGCGAAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAACGGATTATGAGTCCCCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATtat                                 <  1:422974/150‑1 (MQ=255)
                                                                                                                        atttCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAAc                              >  1:1102977/1‑150 (MQ=255)
                                                                                                                             ttGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                        <  2:1911577/87‑1 (MQ=255)
                                                                                                                             ttGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                        >  1:1911577/1‑87 (MQ=255)
                                                                                                                               gCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTCTCGGCTCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGtt                       <  1:155129/150‑1 (MQ=255)
                                                                                                                                   gCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGc                   >  2:131836/1‑150 (MQ=255)
                                                                                                                                    cATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCa                  >  2:1420430/1‑150 (MQ=255)
                                                                                                                                       gacaaaGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCt               >  2:868207/4‑150 (MQ=255)
                                                                                                                                       gAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCt               <  1:2115422/150‑1 (MQ=255)
                                                                                                                                       gAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCt               <  1:2115578/150‑1 (MQ=255)
                                                                                                                                         aaaaGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGTCCCGACACGTGGACTAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCtgt             <  1:2386220/150‑1 (MQ=255)
                                                                                                                                          aaaGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAAGATTATAAGTACTCTTCTCTAACCAACGTGGCTATAGGCCCTCAGTAAGTGGGCGGAAGATAAAGAAGGTTTCGCACCtgtg            >  2:323583/1‑150 (MQ=255)
                                                                                                                                          aaaGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCtgtg            >  2:1699285/1‑150 (MQ=255)
                                                                                                                                            aGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTggg                                         <  2:509034/119‑1 (MQ=255)
                                                                                                                                            aGCCGCGCGATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTggg                                         >  1:509034/1‑119 (MQ=255)
                                                                                                                                              ccGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCAt        <  1:2428872/150‑1 (MQ=255)
                                                                                                                                               cgcgcgATTGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATc       >  1:2069062/1‑150 (MQ=255)
                                                                                                                                               cgcgCAACCGCGCGGCTTTGTCGTATCTGGTGGGTCCTCACGGCCTCGAACCGTGGACCAAAGGATTATGAGTCCTGTGCGCTAACCAACTTAGCTATGGGCTCTCCGTATGTGCGCGCATTATAGCGAGGGTTTCTCACCTGTGCCATc       <  2:398766/150‑1 (MQ=255)
                                                                                                                                                  gcAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAg                           >  1:2211784/1‑127 (MQ=255)
                                                                                                                                                  gcAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAg                           <  2:2211784/127‑1 (MQ=255)
                                                                                                                                                  gcAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCtag    >  1:398766/1‑148 (MQ=255)
                                                                                                                                                    aaTCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGaaa  <  1:1420430/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGCGCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGAAA  >  NC_002947/466953‑467250

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