Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,065 A→G 8.0% intergenic (+45/+58) 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,0650AG 8.0% 323.4 / 12.5 114intergenic (+45/+58)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (52/47);  new base G (3/6);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 120 total aligned reads displayed.

GATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAA  >  NC_002947/5541916‑5542213
                                                                                                                                                     |                                                                                                                                                    
gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                                      <  1:2819310/150‑1 (MQ=255)
gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                                      <  1:1676365/150‑1 (MQ=255)
 aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                                                     <  1:965358/150‑1 (MQ=255)
  tGCCTTCCTGCATGGGCCGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAcc                                                                                                                                                    <  1:222127/150‑1 (MQ=255)
    ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCGGCGAATACGAGGGTGAACCCa                                                                                                                                                  >  2:3037593/1‑150 (MQ=255)
    ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCGGCGAATACGAGGGTGAACCCa                                                                                                                                                  >  2:3038077/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)
      ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAcc                                                                                                                                                >  2:3300004/1‑150 (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: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)
                   aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                   <  1:605955/150‑1 (MQ=25)
                   aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                   <  1:2784615/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)
                        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)
                             ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTTCACCCGc                                                                                                                         <  2:4701742/150‑1 (MQ=255)
                              cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGGATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCt                                                                                                                        >  2:1789231/1‑150 (MQ=255)
                                          gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTc                                                                                                            >  2:4215708/1‑150 (MQ=255)
                                           aCATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCGCGCAAACGAGGGTGTACCCACCCACGCAGTCGCCGGTGAACCCCCACCGGCAGGGGTCg                                                                                                           >  2:2083588/1‑150 (MQ=255)
                                             atcatcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                                         >  2:2716692/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)
                                                                     aCCCCGCCCACGCGCGAACTGCCGTTCGCTCATTTATACTTGGCCTCCAGGCGGGGTTTCACCCGCTAATACGCGGGTGAACCCACCTACGCATGCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                 <  2:2106762/150‑1 (MQ=255)
                                                                     aCCCAGCACACGCGCTAACTGCCGTTCGCTCATTGATATGTGGCCTGCTGGAGGGGTTTCACCCGCGAATACGCGGGTTAACCCACCTACGCTTTGGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                                 <  2:2106754/150‑1 (MQ=255)
                                                                       ccAGCACCAGCGCTAAGTGCCGTTCGCTCCTTGATAGGTGGCCTGCTGGAGCGTGTTCACCCGCTAATACGTGGGTCACCCCACCGCCGCTTTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                               <  2:637311/150‑1 (MQ=255)
                                                                         aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                   >  1:4113068/1‑96 (MQ=255)
                                                                         aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                   <  2:4113068/96‑1 (MQ=255)
                                                                                  gccaagtGCGGTTCGCTCGTTGATATGTGGCCTGCTGGCGGGGGTTCCCCCGCTAATCGGTGGTTCACCCCACCGACGCTTTGGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                    <  1:1287917/144‑1 (MQ=255)
                                                                                     aaCTGCCGTTCGCTCCTTGTTATGTGGCCTGCTGGAGCGGTTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcatc                                                                 <  1:5113883/150‑1 (MQ=38)
                                                                                        tGCGGTTCGCTCCTTGATATGTGGCCTGCTGGCGCGGTTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                              <  1:4136881/150‑1 (MQ=38)
                                                                                        tGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCCCTCCAGCACGGTTCGCCCTCAGCCTTTGATTACTTGACGATGACCGTTCCGTTAatcatc                                                              >  2:3738588/1‑150 (MQ=38)
                                                                                        tGCCGTGCGCTCTTTTATACGTGGCCTGCCGGCGCGGGTCCCCCCGAGAATCCGCGGGACACCCCACCGCCGCTTTGGGGGGCTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                              <  1:3370679/150‑1 (MQ=255)
                                                                                         gCCGTTCGCTCCTTGATAAGTGGCCTGCTGGAGGGGGTTCACCCGCTAATACGTGGTTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                             <  1:4679971/150‑1 (MQ=39)
                                                                                         gCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                             <  2:3743065/150‑1 (MQ=255)
                                                                                                         tctgTGGCCGGCTGGCGGGGTTCCCCCCGCTATTCCGGGGTTCAACCCCCCGCCGCATTCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGaa                                             <  1:2083588/147‑1 (MQ=255)
                                                                                                            gTGGCCTGCTGGCGGGTTGTCACCCCCGATTACGTGGGTAAACCCCCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGt                                          <  1:1407211/150‑1 (MQ=255)
                                                                                                              ggCCTGCAGGAGCGGTTCCACCCGCTAATCCGCGGGTCAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                        <  2:1121030/150‑1 (MQ=255)
                                                                                                                    ccggcgCGGGTTCACCCGCGAATACGCGGGTGAACCCCCCGACGCATTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                                  <  2:4223301/145‑1 (MQ=255)
                                                                                                                    cAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGTCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                                  <  2:4620492/150‑1 (MQ=255)
                                                                                                                      ggtgCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                <  2:330540/147‑1 (MQ=255)
                                                                                                                       gcgCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                               <  2:1631683/148‑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)
                                                                                                                             gTTCACCCTCGTATTCGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                         >  2:1398056/1‑150 (MQ=255)
                                                                                                                             gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                         >  1:5431820/1‑150 (MQ=255)
                                                                                                                              ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                        >  1:997555/1‑150 (MQ=255)
                                                                                                                              ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                        >  1:997557/1‑150 (MQ=255)
                                                                                                                                cACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGAAGAAGCTGTAGTcac                      >  1:4760007/1‑150 (MQ=25)
                                                                                                                                 aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                     >  2:5131970/1‑150 (MQ=25)
                                                                                                                                  cccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGCTGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                    >  2:2495275/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCGGGGACGGTGCAGAAGTAGCTGTAGTcaccac                   >  2:3515020/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                   >  1:1023051/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                   >  2:495465/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                   >  1:625855/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                   >  2:1290767/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                   >  2:2173043/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                   >  1:1596616/1‑150 (MQ=255)
                                                                                                                                   ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                   >  1:2244828/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)
                                                                                                                                           aCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTcc           >  2:5493535/1‑150 (MQ=255)
                                                                                                                                             gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg         >  1:830839/1‑150 (MQ=255)
                                                                                                                                             gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg         >  2:3999457/1‑150 (MQ=255)
                                                                                                                                                 gttaaCCCACCGACCCATTCGCGGGTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg     <  1:4161076/147‑1 (MQ=255)
                                                                                                                                                 gTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg     <  2:5168799/150‑1 (MQ=255)
                                                                                                                                                 gTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg     <  1:4145889/150‑1 (MQ=255)
                                                                                                                                                    aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa  <  2:216671/150‑1 (MQ=255)
                                                                                                                                                    aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa  >  2:3159924/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAA  >  NC_002947/5541916‑5542213

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: