Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,355 G→C 10.8% intergenic (+219/‑597) tyrS → / → PP_16SC tyrosine‑‑tRNA ligase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947524,3550GC10.8% 212.5 / 14.6 72intergenic (+219/‑597)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base C (5/3);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACA  >  NC_002947/524209‑524503
                                                                                                                                                  |                                                                                                                                                    
ggtTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCACCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCgcgggcg                                                                                                                                                   >  1:985234/1‑150 (MQ=255)
      tgttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGCCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                             <  1:100296/150‑1 (MQ=255)
       gttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                              >  1:1430522/1‑148 (MQ=255)
       gttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                              <  2:1430522/148‑1 (MQ=255)
          gtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtc                                                                                                                                         >  1:1674692/1‑150 (MQ=255)
           tAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtct                                                                                                                                        >  1:1273592/1‑150 (MQ=255)
           tAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGAGCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtct                                                                                                                                        >  1:1273585/1‑150 (MQ=255)
            agggggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtctc                                                                                                                                       <  2:1621445/150‑1 (MQ=255)
            agggggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtctc                                                                                                                                       <  2:1219058/150‑1 (MQ=255)
             gggggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCa                                                                                                                                      <  2:1674692/150‑1 (MQ=255)
              ggggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtctc                                                                                                                                       <  1:1144594/148‑1 (MQ=255)
              ggggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtctc                                                                                                                                       >  2:1144594/1‑148 (MQ=255)
                    ggggCTGTTTTGAGGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCTAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCt                                                                                                                              >  2:1636185/1‑150 (MQ=255)
                          gTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAAc                                                                                                                         >  2:1517230/1‑150 (MQ=255)
                                  agGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTgg                                                                                                                 <  2:1715961/150‑1 (MQ=255)
                                           gCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                                        <  2:1770966/150‑1 (MQ=255)
                                                    tgtgTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGAGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTc                                                                                               <  2:1255982/150‑1 (MQ=255)
                                                     gtgtCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                             >  2:726999/1‑103 (MQ=25)
                                                     gtgtCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                             <  1:726999/103‑1 (MQ=25)
                                                      tgtCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGa                                                                                             >  2:299313/1‑150 (MQ=255)
                                                      tgtCCGCCACAACCCGTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTGGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGa                                                                                             <  1:108322/150‑1 (MQ=255)
                                                        tCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCTGGGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATa                                                                                           >  1:313025/1‑150 (MQ=255)
                                                           gCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACgg                                                                                        >  2:419729/1‑150 (MQ=255)
                                                                 aCCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCtttt                                                                                  >  2:304995/1‑150 (MQ=255)
                                                                  cccTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTc                                                                                 >  1:152895/1‑150 (MQ=255)
                                                                  cccTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTc                                                                                 <  1:1193710/150‑1 (MQ=255)
                                                                    cTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTccc                                                                               >  1:525107/1‑150 (MQ=255)
                                                                    cTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTccc                                                                               >  1:628534/1‑150 (MQ=255)
                                                                    cTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTccc                                                                               <  1:287142/150‑1 (MQ=255)
                                                                      tGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCt                                                                             <  1:419729/150‑1 (MQ=255)
                                                                                gtgAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCAc                                                                   <  1:1517230/150‑1 (MQ=255)
                                                                                  gagaTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACtt                                                                 >  1:930513/1‑150 (MQ=255)
                                                                                    gaTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                                               >  2:1720650/1‑150 (MQ=255)
                                                                                    gaTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                                               >  2:1587962/1‑150 (MQ=255)
                                                                                     aTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                              <  2:481889/150‑1 (MQ=255)
                                                                                        gAGCGCCGCCCGCGTCCGGCTTGAGATTTGAGCCGCCTATGAAACCGAGCGCCGCGCGGGAGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGa                                                           >  1:49907/1‑150 (MQ=18)
                                                                                        gAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGa                                                           >  2:1395070/1‑150 (MQ=255)
                                                                                        gAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGa                                                           >  2:70338/1‑150 (MQ=255)
                                                                                          gcgcCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                                         >  1:1445952/1‑150 (MQ=255)
                                                                                                   gcgTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                                <  1:455334/150‑1 (MQ=255)
                                                                                                     gTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGa                                              >  1:545537/1‑150 (MQ=255)
                                                                                                      tCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAg                                                   >  1:1507023/1‑144 (MQ=255)
                                                                                                      tCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAg                                                   <  2:1507023/144‑1 (MQ=255)
                                                                                                      tCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTg                                               >  1:1751671/1‑148 (MQ=255)
                                                                                                      tCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTg                                               <  2:1751671/148‑1 (MQ=255)
                                                                                                         ggCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGcc                                          >  2:1670221/1‑150 (MQ=255)
                                                                                                             tGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCctct                                      <  2:1273592/150‑1 (MQ=255)
                                                                                                             tGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCctct                                      <  2:1273585/150‑1 (MQ=255)
                                                                                                             tGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGGCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCctct                                      <  1:1276703/150‑1 (MQ=255)
                                                                                                              gagaCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCtctc                                     <  2:152895/150‑1 (MQ=255)
                                                                                                                gaCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGa                                                                                             <  2:1570758/92‑1 (MQ=255)
                                                                                                                gaCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGGGGCGAACACTTGGCAGCCTCGa                                                                                             >  1:1570758/1‑92 (MQ=255)
                                                                                                                 aCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGaa                                  >  1:1752947/1‑150 (MQ=255)
                                                                                                                   aTTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATc                                >  2:208994/1‑150 (MQ=255)
                                                                                                                     tCAGCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATccc                              >  1:367109/1‑150 (MQ=255)
                                                                                                                     tCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATccc                              <  2:985234/150‑1 (MQ=255)
                                                                                                                       aCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                            >  1:1460046/1‑150 (MQ=255)
                                                                                                                       aCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                            >  1:716194/1‑150 (MQ=255)
                                                                                                                        ccGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCtt                           <  2:628534/150‑1 (MQ=255)
                                                                                                                         cGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTatat                                                                        >  2:772208/1‑104 (MQ=255)
                                                                                                                         cGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTatat                                                                        <  1:772208/104‑1 (MQ=255)
                                                                                                                         cGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTta                          >  1:75367/1‑150 (MQ=255)
                                                                                                                         cGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTta                          >  2:52485/1‑150 (MQ=255)
                                                                                                                           ccTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata                        >  2:1318761/1‑150 (MQ=255)
                                                                                                                           ccTATGAAACCGAGCGCCGCGAGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata                        <  1:1895975/150‑1 (MQ=255)
                                                                                                                            cTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATaa                       >  1:412015/1‑150 (MQ=255)
                                                                                                                               tGAAACCCAGCAGCGCTCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAGCATGGGGCAGCCTCGATGCGCTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgc                    <  1:1076796/150‑1 (MQ=14)
                                                                                                                                  gaCCGAGCGCCGCCCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcc                 >  1:1227337/2‑150 (MQ=255)
                                                                                                                                   aCCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGccc                >  1:113122/1‑150 (MQ=11)
                                                                                                                                         cgccgcGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTc          >  1:754385/1‑150 (MQ=255)
                                                                                                                                         cgccgcGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTc          <  1:1587962/150‑1 (MQ=255)
                                                                                                                                          gccgcGCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTcc         >  2:292475/1‑150 (MQ=255)
                                                                                                                                             gcgcgGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc      <  1:208994/150‑1 (MQ=255)
                                                                                                                                                 gggcggCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACAGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACa  >  2:836690/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACA  >  NC_002947/524209‑524503

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: