Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 467,104 C→T 8.4% intergenic (‑34/+19) 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,1040CT 8.4% 228.1 / 9.4 83intergenic (‑34/+19)PP_5445/PP_t71hypothetical protein/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base C (34/42);  new base T (5/2);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGAAA  >  NC_002947/466957‑467250
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                                                <  2:434316/139‑1 (MQ=255)
             ggCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATcgcgc                                                                                                                                                >  1:434316/1‑139 (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                                                                                                                 <  2:2436934/150‑1 (MQ=255)
                                 ttCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGa                                                                                                                 <  1:1008039/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:1312983/1‑150 (MQ=255)
                                              cGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGa                                                                                                    <  1:5731/150‑1 (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                                                                                          >  2:1725537/1‑150 (MQ=255)
                                                        cTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGAtt                                                                                          <  1:1725537/150‑1 (MQ=255)
                                                          aCTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                        <  2:1831057/150‑1 (MQ=255)
                                                          aCTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCCGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTAt                                                                                        <  2:1235570/150‑1 (MQ=255)
                                                           cTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                    >  2:1175926/1‑121 (MQ=255)
                                                           cTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCacac                                                                                                                    <  1:1175926/121‑1 (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                                                                                   >  1:1460846/1‑150 (MQ=255)
                                                               cAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTc                                                                                   <  2:5930/150‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
GCGCTGGGTGACCGGCTTCTTGGTGGGGCTGATTTCGCTGAGATACCGTTTTAGCGCTACTTCCAAGGTCATGCGCTCGGCGCCGCTGCGCTGGATGTACGCGCCCCTGACCATTTATTTCTTGCGGGCATGAAAAAGCCGCGCAATCGCGCGGCTTTGTCGTATCTGGTGGGCCCACACGGACTCGAACCGTGGACCAAAGGATTATGAGTCCTCTGCTCTAACCAACTGAGCTATAGGCCCTCAGTAAGTGGGCGGATTATAACGAGGGTTTCGCACCTGTGCCATCTGAAA  >  NC_002947/466957‑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: