Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,067 C→T 6.9% intergenic (+47/+56) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,542,0670CT 6.9% 353.3 / 11.8 116intergenic (+47/+56)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base C (53/54);  new base T (2/6);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 124 total aligned reads displayed.

TGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTG  >  NC_002947/5541918‑5542216
                                                                                                                                                     |                                                                                                                                                     
tGCCTTCCTGCATGGGCCGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAcc                                                                                                                                                       <  1:222127/150‑1 (MQ=255)
  ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCGGCGAATACGAGGGTGAACCCa                                                                                                                                                     >  2:3037593/1‑150 (MQ=255)
  ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCa                                                                                                                                                     <  1:2077067/150‑1 (MQ=255)
  ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCa                                                                                                                                                     <  1:2076984/150‑1 (MQ=255)
   cTTCCTGCATGGGCAGCCGGTTGTCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAc                                                                                                                                                    <  2:637945/150‑1 (MQ=255)
    ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGAACCGCGAATACGATGGTGAACCCagc                                                                                                                                                   >  2:4161076/1‑148 (MQ=255)
    ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGAGCCCCGCAAAAGAGGGGGTACACAcc                                                                                                                                                   >  1:1116608/1‑150 (MQ=255)
    ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACACAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAcc                                                                                                                                                   <  2:3285406/150‑1 (MQ=255)
       cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                                <  1:3300004/150‑1 (MQ=255)
       cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                                >  2:5487073/1‑150 (MQ=255)
       cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                                >  2:5487005/1‑150 (MQ=255)
        tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCACGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACg                                                                                                                                               <  1:4935399/150‑1 (MQ=255)
        tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACg                                                                                                                                               <  1:3073708/150‑1 (MQ=255)
        tGCATGGGCAGCCGGTGGCCCAGGCAGCGTTCGACACCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACg                                                                                                                                               <  1:3073445/150‑1 (MQ=255)
         gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCCCGAATACGAGGGTGTACCCACCCACGc                                                                                                                                              >  2:2116987/1‑150 (MQ=255)
          cATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCa                                                                                                                                             <  2:4679773/150‑1 (MQ=255)
            tGGGCAGCCGGTGGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAtt                                                                                                                                           <  1:193373/150‑1 (MQ=255)
               gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGAACCGCGAAAAAGAGGGTGTAACCGCCGACGCAttagg                                                                                                                                        >  2:1407211/1‑147 (MQ=255)
               gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                        <  1:3766627/150‑1 (MQ=255)
                cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgcg                                                                                                                                       <  1:1542484/150‑1 (MQ=255)
                cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAACCACCGACGCATTcgcg                                                                                                                                       >  2:4476792/1‑150 (MQ=255)
                 aGCCGGTTGACCTGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACATGCGCGAACTGCCGTTCGCTCCTTGATCAGTGGCCTGCAGGAGCTGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  2:3474629/150‑1 (MQ=21)
                 aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                      <  1:1039449/150‑1 (MQ=25)
                 aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACTAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  2:2783882/150‑1 (MQ=255)
                 aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  2:3345262/150‑1 (MQ=255)
                 aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  1:4611341/150‑1 (MQ=255)
                 aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCACCACAAGCGCGAACTGCCGTTCGCTCCTTGATAACTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  2:3345407/150‑1 (MQ=255)
                   ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCCCGAATAACAGGGTGAACCCACCGACCCAATAGCGGGt                                                                                                                                    >  2:4390191/1‑150 (MQ=255)
                      gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGaa                                                                                                                                 >  2:4145889/1‑150 (MQ=255)
                       ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGAACCCCCAAAAACAGGGGGAAACCACCCACCCAATAGCCGGTGAAc                                                                                                                                >  1:4931295/1‑150 (MQ=255)
                       ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                                <  1:1568255/150‑1 (MQ=25)
                       ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                                <  1:1460192/150‑1 (MQ=25)
                       ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAc                                                                                                                                <  2:1928662/150‑1 (MQ=255)
                         gACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAccc                                                                                                                              >  2:418649/1‑150 (MQ=255)
                            cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGGATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCt                                                                                                                           >  2:1789231/1‑150 (MQ=255)
                                     ttCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCCggg                                                                                                                  >  2:5044502/1‑150 (MQ=255)
                                        gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTc                                                                                                               >  2:4215708/1‑150 (MQ=255)
                                         aCATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCGCGCAAACGAGGGTGTACCCACCCACGCAGTCGCCGGTGAACCCCCACCGGCAGGGGTCg                                                                                                              >  2:2083588/1‑150 (MQ=255)
                                           aacatcGTCGCCACCTCCCGCACCACCCACCACACGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCTGGAGCGGGTTCACCCGCGAATACGTGGGTTAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                                            <  2:2686957/148‑1 (MQ=37)
                                                     ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGAACCCCCGATAACAGGGTGTACCCCCCGCCGCCATCGCCGGTGTAACCCCGCCTGCAGGGTTGGCCGGCCGCCgtt                                                                                                  >  2:1477022/1‑150 (MQ=255)
                                                       accgcccGCACGCCCCAGCCCACGCGCTACCGGCGGTTCGCTCCGTGATAAGGGGCCTGCAGGAGGGGGTTCACCGGCGAATACGAGGGTGACCCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGa                                                                                                <  2:3018209/145‑1 (MQ=25)
                                                       aCCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCGCGAATACGAGGGTGAACCCACCGACCCAATCGCGGGTGAACCCGCTCCTGCAGGGTTAGCCGTCAGCCTTTGa                                                                                                >  2:5055149/1‑150 (MQ=255)
                                                           cccGAAAGCCCCAGCACACGCGCTAACTGCCGTTCGCCCCTTGAAAAGTGGCCTGCAGGCGCGGGTTAACCCGCGATTACGAGGGTGAACCAACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAc                                                                                            <  2:570887/149‑1 (MQ=255)
                                                             cgcaccaCCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGTTATGTGGCCTGCAGGCGGGGGTTCACCCGCTATTCCGCGGGTTAACCCACCGACGCATTGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                          <  2:1052564/144‑1 (MQ=255)
                                                               cAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCCCGCAAAACCGGGTGTGCCCACCGACCCCGGCGCCGGGGTACCCCCGCCGGCCGGGTTCGCCCGGAGCCGTTGTGTACTTgt                                                                                        >  2:1433423/1‑149 (MQ=255)
                                                                aacaCCCCGAACACGCGCGAACTGCCGTTCGCTCCTTCATACGGGGCCTGCAGGCGCGGGTACACCCGCTAATACGCGGGTGACCCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                                       <  2:2835113/147‑1 (MQ=37)
                                                                   ccccAGCACACGCGCGAACTGCGGTTCGCTCCTTTAAACGTGGCCTGCAGGTGCGGTTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                    <  2:2106707/149‑1 (MQ=255)
                                                                   aCCCCGCCCACGCGCGAACTGCCGTTCGCTCATTTATACTTGGCCTCCAGGCGGGGTTTCACCCGCTAATACGCGGGTGAACCCACCTACGCATGCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                    <  2:2106762/150‑1 (MQ=255)
                                                                     ccAGCACCAGCGCTAAGTGCCGTTCGCTCCTTGATAGGTGGCCTGCTGGAGCGTGTTCACCCGCTAATACGTGGGTCACCCCACCGCCGCTTTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                                  <  2:637311/150‑1 (MQ=255)
                                                                       aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      <  2:4113068/96‑1 (MQ=255)
                                                                       aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                      >  1:4113068/1‑96 (MQ=255)
                                                                                gccaagtGCGGTTCGCTCGTTGATATGTGGCCTGCTGGCGGGGGTTCCCCCGCTAATCGGTGGTTCACCCCACCGACGCTTTGGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                       <  1:1287917/144‑1 (MQ=255)
                                                                                      tGCGGTTCGCTCCTTGATATGTGGCCTGCTGGCGCGGTTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                 <  1:4136881/150‑1 (MQ=38)
                                                                                      tGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCCCTCCAGCACGGTTCGCCCTCAGCCTTTGATTACTTGACGATGACCGTTCCGTTAatcatc                                                                 >  2:3738588/1‑150 (MQ=38)
                                                                                       gCCGTTCGCTCCTTGATAAGTGGCCTGCTGGAGGGGGTTCACCCGCTAATACGTGGTTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                                <  1:4679971/150‑1 (MQ=39)
                                                                                       gCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                                <  2:3743065/150‑1 (MQ=255)
                                                                                              gccccTTGATACGGGGCCTGCAGGCGGGGGTCCACCGGCGAATACGCGGGTGAACCCACCGACGCATTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGAtgtg                                                         <  2:1116608/147‑1 (MQ=255)
                                                                                                          gTGGCCTGCTGGCGGGTTGTCACCCCCGATTACGTGGGTAAACCCCCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGt                                             <  1:1407211/150‑1 (MQ=255)
                                                                                                            ggCCTGCAGGAGCGGTTCCACCCGCTAATCCGCGGGTCAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                           <  2:1121030/150‑1 (MQ=255)
                                                                                                                  cAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGTCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                                     <  2:4620492/150‑1 (MQ=255)
                                                                                                                    ggtgCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                   <  2:330540/147‑1 (MQ=255)
                                                                                                                    ggAGCGGGTTCACCCGCGAATAGGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                   <  2:329873/150‑1 (MQ=255)
                                                                                                                       gCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                >  1:4306570/1‑150 (MQ=255)
                                                                                                                       gCGGGTCCACCCGCGAATACGCGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                <  1:5125120/150‑1 (MQ=255)
                                                                                                                         ggtttCCCCGGCTAATCCGTGGGTTAACCCACCGCCGCTTGCGGGGTTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                              <  1:2396575/147‑1 (MQ=255)
                                                                                                                         gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                              >  1:2044872/1‑150 (MQ=255)
                                                                                                                           gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                            >  1:5431820/1‑150 (MQ=255)
                                                                                                                            ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                           >  1:997555/1‑150 (MQ=255)
                                                                                                                                cccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGCTGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                       >  2:2495275/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCGGGGACGGTGCAGAAGTAGCTGTAGTcaccac                      >  2:3515020/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:1023051/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:76643/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:625855/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:1290767/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:1596616/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:1467427/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:495465/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:5073629/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:2173043/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:2244828/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAGTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  1:2829943/1‑150 (MQ=255)
                                                                                                                                 ccGCGAATACGAGGGTGAACCCACCGACGCAGTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                      >  2:2406062/1‑150 (MQ=255)
                                                                                                                                  cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                     >  2:5325079/1‑150 (MQ=255)
                                                                                                                                  cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc                     >  2:581068/1‑150 (MQ=255)
                                                                                                                                         aCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTcc              >  1:4399320/1‑150 (MQ=255)
                                                                                                                                           gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg            >  1:830839/1‑150 (MQ=255)
                                                                                                                                           gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg            >  2:3999457/1‑150 (MQ=255)
                                                                                                                                               gttaaCCCACCGACTAATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg        <  1:2116987/147‑1 (MQ=255)
                                                                                                                                               gttaaCCCACCGACCCATTCGCGGGTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg        <  1:4161076/147‑1 (MQ=255)
                                                                                                                                               gTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg        <  2:5168799/150‑1 (MQ=255)
                                                                                                                                               gTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg        <  1:3999457/150‑1 (MQ=255)
                                                                                                                                                  aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa     <  2:216671/150‑1 (MQ=255)
                                                                                                                                                  aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa     >  2:3159924/1‑150 (MQ=255)
                                                                                                                                                    cccACCGACGCATTCGCGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt   <  1:5487073/150‑1 (MQ=255)
                                                                                                                                                    cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt   >  2:4707999/1‑150 (MQ=255)
                                                                                                                                                    cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt   <  1:5487005/150‑1 (MQ=255)
                                                                                                                                                    cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAAGATGACCTTGCCTTTCATCATAGAGATGTGGCCCGGTAAGGTGAATAAGAAGCTGTAGTCAACACAGGATTCCAGCTTGGAAgt   >  2:1768150/1‑150 (MQ=255)
                                                                                                                                                     ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg  >  2:2305812/1‑150 (MQ=255)
                                                                                                                                                     ccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg  >  2:3218979/1‑150 (MQ=255)
                                                                                                                                                     ccACCGACGCATTCGCGGGTGAACCCGCTCCTACAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtg  >  1:4870263/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTG  >  NC_002947/5541918‑5542216

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