Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,634,468 A→G 7.4% intergenic (+105/+54) PP_4096 → / ← PP_t59 hypothetical protein/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,634,4680AG 7.4% 440.5 / 15.8 150intergenic (+105/+54)PP_4096/PP_t59hypothetical protein/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base A (67/71);  new base G (6/5);  total (74/76)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

Only 100 of 157 total aligned reads displayed.

TTTGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTC  >  NC_002947/4634322‑4634616
                                                                                                                                                  |                                                                                                                                                    
tttGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAggg                                                                                                                                                   >  1:4934541/1‑150 (MQ=255)
tttGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAggg                                                                                                                                                   >  1:1479034/1‑150 (MQ=255)
tttGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAggg                                                                                                                                                   >  1:1512576/1‑150 (MQ=255)
tttGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAggg                                                                                                                                                   >  2:220091/1‑150 (MQ=255)
  tGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGcc                                                                                                                                                 >  2:227367/1‑150 (MQ=255)
    aaaTACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCtt                                                                                                                                               <  2:5380126/150‑1 (MQ=255)
       tACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCttttt                                                                                                                                            <  2:779909/150‑1 (MQ=255)
        aCCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCCGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCtttttt                                                                                                                                           >  1:87041/1‑150 (MQ=255)
         cccGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCTTCGGGGCCTTTTTTc                                                                                                                                          <  1:2180604/150‑1 (MQ=255)
         cccGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCTTCGGGGCCTTTTTTc                                                                                                                                          <  2:4676931/150‑1 (MQ=255)
          ccGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCa                                                                                                                                         >  1:1971808/1‑150 (MQ=255)
           cGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCaa                                                                                                                                        <  2:2392864/150‑1 (MQ=255)
           cGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCaa                                                                                                                                        >  1:2554421/1‑150 (MQ=255)
           cGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCaa                                                                                                                                        <  1:1747155/150‑1 (MQ=255)
           cGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCaa                                                                                                                                        <  1:1746961/150‑1 (MQ=255)
              tGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGa                                                                                                                                     <  2:2066249/150‑1 (MQ=255)
                aCGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACt                                                                                                                                   >  1:4868208/1‑150 (MQ=255)
                 cGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACtt                                                                                                                                  >  1:3812583/1‑150 (MQ=255)
                 cGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACtt                                                                                                                                  >  1:1935519/1‑150 (MQ=255)
                  gTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTa                                                                                                                                 >  1:4299407/1‑150 (MQ=255)
                      tCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGa                                                                                                                             >  1:4384593/1‑150 (MQ=255)
                        aGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACt                                                                                                                           >  2:5151098/1‑150 (MQ=255)
                            atatGTCCAGGTAGCCCACGGCGTCACGGCCTGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTTGAGGGCGGGGTTTTGAGGCAGAATGCATTCGAACTAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGa                                                                                                                       <  1:3783077/150‑1 (MQ=255)
                            atatGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGa                                                                                                                       <  2:2554421/150‑1 (MQ=255)
                            atatGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGa                                                                                                                       >  1:744916/1‑150 (MQ=255)
                            atatGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGa                                                                                                                       <  2:4384593/150‑1 (MQ=255)
                             tatGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGAc                                                                                                                      <  2:2517955/150‑1 (MQ=255)
                                 tCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGAc                                                                                                                      <  1:3145756/146‑1 (MQ=255)
                                 tCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGAc                                                                                                                      >  2:3145756/1‑146 (MQ=255)
                                      gTAGCCCACGGCGTCACGGCCTGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTg                                                                                                             <  2:87041/150‑1 (MQ=255)
                                            cACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATc                                                                                                       <  2:1024985/150‑1 (MQ=255)
                                                      cGGCCTGGGCGCCGCGTTGTGTGGTTTGTACTTGTCTGTCGTAGAGGGCGGGGTATTGTGGCAGAATGCAGTCGAACGAAAAAAGGCCCCGAAGGGCCCTTTTTCTACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGa                                                                                             <  2:181033/150‑1 (MQ=255)
                                                         ccAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCg                                                                                          <  1:2578816/150‑1 (MQ=255)
                                                         ccAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCg                                                                                          >  2:1987230/1‑150 (MQ=255)
                                                         ccAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCg                                                                                          <  2:5411874/150‑1 (MQ=255)
                                                           aGGGCGCCGCGTTGTTGTGTTTGTACTAGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCGTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGAGCAGTATTTGGAGgggg                                                                                        >  2:3819641/1‑150 (MQ=255)
                                                             ggCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCTAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCg                                                                                          <  2:2973225/146‑1 (MQ=255)
                                                              gcgcCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACCGAAACCTGTAGATCTTTATTTGGAGCGGGaaa                                                                                     >  1:2011888/1‑150 (MQ=255)
                                                               cgccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGacac                                                                                    >  1:3656237/1‑150 (MQ=255)
                                                               cgccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAAc                                                                                    >  1:3144467/1‑150 (MQ=255)
                                                                 ccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGTGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACga                                                                                  <  2:3572660/150‑1 (MQ=255)
                                                                 ccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACga                                                                                  <  1:5151098/150‑1 (MQ=255)
                                                                 ccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACga                                                                                  <  1:4511790/150‑1 (MQ=255)
                                                                 ccgcGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACga                                                                                  <  2:2692477/150‑1 (MQ=255)
                                                                    cgTTGTGGTTTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGAc                                                                               <  2:1560501/150‑1 (MQ=255)
                                                                      ttttggtgTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTc                                                                             <  2:693664/147‑1 (MQ=255)
                                                                            tgtTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTc                                                                       <  2:60236/150‑1 (MQ=255)
                                                                              tttGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCtttt                                                                                                                                             >  1:1222185/1‑78 (MQ=255)
                                                                              tttGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCtttt                                                                                                                                             <  2:1222185/78‑1 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCTGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCg                                                              >  2:202804/1‑150 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCt                                                                                                      <  1:3356447/110‑1 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTg                                                                                               <  2:2735710/117‑1 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTg                                                                                               >  1:2735710/1‑117 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCg                                                              >  2:2345210/1‑150 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCg                                                              >  2:2344893/1‑150 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCg                                                              <  2:1348348/150‑1 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCg                                                              >  2:1111386/1‑150 (MQ=255)
                                                                                     ttGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCg                                                              >  2:1111254/1‑150 (MQ=255)
                                                                                      tGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGa                                                             >  1:713103/1‑150 (MQ=255)
                                                                                       gAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTCGATCTTTATGTGGAGCGGGAAACGAGACTCGCACTCGCGCCCCCgcc                                                            >  1:1185391/1‑148 (MQ=255)
                                                                                          tttCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCtt                                                         >  1:2566509/1‑150 (MQ=255)
                                                                                           ttCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTg                                                        <  2:4832110/150‑1 (MQ=255)
                                                                                            tCGTAGAGGGCGGTGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                       <  1:3956601/150‑1 (MQ=255)
                                                                                            tCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                       <  1:491733/150‑1 (MQ=255)
                                                                                                   gggCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                                <  2:3878452/150‑1 (MQ=255)
                                                                                                     gCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGt                                              <  1:3948113/150‑1 (MQ=255)
                                                                                                       ggggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGc                                            <  2:4794749/150‑1 (MQ=255)
                                                                                                         ggTATTGAGGCATAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctc                                          <  1:2495347/150‑1 (MQ=255)
                                                                                                         ggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctc                                          <  1:2494687/150‑1 (MQ=255)
                                                                                                            aTTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAc                                       <  1:1111386/150‑1 (MQ=255)
                                                                                                            aTTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTAc                                       <  1:1111254/150‑1 (MQ=255)
                                                                                                               gAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGTAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTTTTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCaa                                    <  2:4413080/150‑1 (MQ=255)
                                                                                                                aGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAAc                                   >  2:668037/1‑150 (MQ=255)
                                                                                                                aGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAAc                                   >  2:667972/1‑150 (MQ=255)
                                                                                                                aGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAAc                                   >  1:2402003/1‑150 (MQ=255)
                                                                                                                   cAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                <  1:220091/150‑1 (MQ=255)
                                                                                                                         gCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTAtt                          <  2:1512576/150‑1 (MQ=255)
                                                                                                                           aTTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTcc                        >  1:4160310/1‑150 (MQ=255)
                                                                                                                           aTTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTcc                        >  1:4160104/1‑150 (MQ=255)
                                                                                                                            ttCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTccc                       >  1:849249/1‑150 (MQ=255)
                                                                                                                             tCGAACGAAAAAAGGCCCCTAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCaa                                                    <  2:1990175/120‑1 (MQ=255)
                                                                                                                             tCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCaa                                                    >  1:1990175/1‑120 (MQ=255)
                                                                                                                             tCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcg                      >  2:4017038/1‑150 (MQ=255)
                                                                                                                              cgaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcgc                     <  2:5240192/150‑1 (MQ=255)
                                                                                                                               gaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTGTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcgcg                    >  2:5449305/1‑150 (MQ=255)
                                                                                                                                gttgaaAAAAGGCCCTTCGGGGCCTTTTTTCGACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGcgt                   >  2:3596247/4‑150 (MQ=21)
                                                                                                                                 ttgaaAAAAGGCCCTTCGGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGcgtc                  >  1:3468915/3‑150 (MQ=255)
                                                                                                                                   gaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGcgtcgt                >  1:5311417/1‑150 (MQ=255)
                                                                                                                                      aaaaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGAt             <  2:5362295/150‑1 (MQ=255)
                                                                                                                                       aaaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATg            <  2:2638106/150‑1 (MQ=255)
                                                                                                                                        aaGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATgg           <  2:3655185/150‑1 (MQ=255)
                                                                                                                                          ggCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGt         <  1:1872389/150‑1 (MQ=255)
                                                                                                                                             cccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGcc      >  1:2659796/1‑150 (MQ=255)
                                                                                                                                              ccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCa     >  1:4814648/1‑150 (MQ=255)
                                                                                                                                               cGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCAt    >  2:4186949/1‑150 (MQ=255)
                                                                                                                                                 aaGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGa                                                                   <  2:2213978/85‑1 (MQ=255)
                                                                                                                                                 aaGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGa                                                                   >  1:2214006/1‑85 (MQ=255)
                                                                                                                                                 aaGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTAtt                          <  2:2453142/126‑1 (MQ=255)
                                                                                                                                                 aaGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTAtt                          >  1:2453142/1‑126 (MQ=255)
                                                                                                                                                 aaGGGACTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTc  <  1:5394662/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTTGAAATACCCGTTGACGTGATCAGCGATATGTCCAGGTAGCCCACGGCGTCACGGCCAGGGCGCCGCGTTGTGGTGTTTGTACTTGAATTTCGTAGAGGGCGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTC  >  NC_002947/4634322‑4634616

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