Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,473,481:1 +G 5.6% coding (122/693 nt) nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,473,4811.G 5.6% 527.5 / 22.5 125coding (122/693 nt)nadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base . (55/63);  new base G (3/4);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 129 total aligned reads displayed.

GGGCATCGACGCTCAGGCAAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATCACCACATACTTCTCGCAGG  >  NC_002947/5473332‑5473629
                                                                                                                                                       |                                                                                                                                                     
gggCATCGACGCTCAGGCAAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACtt                                                                                                                                                        >  1:1800891/1‑150 (MQ=255)
   cagcGACTCTCAGGCCAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCa                                                                                                                                                   <  2:647525/147‑1 (MQ=255)
            tCAGGCAAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCGGGGCATTGGGCAACAGGCGTAGCTCATCCAGCACCATGAACTCGGACACTTACAGCGCGCTgc                                                                                                                                          >  1:280806/1‑150 (MQ=255)
               ggCAAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCa                                                                                                                                       >  1:1195840/1‑150 (MQ=255)
                gCAAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCaa                                                                                                                                      >  1:1967501/1‑150 (MQ=255)
                gCAAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCaa                                                                                                                                      <  2:1715934/150‑1 (MQ=255)
                   aGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTg                                                                                                                                   <  2:2807265/150‑1 (MQ=255)
                    gCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTgg                                                                                                                                  >  1:1501725/1‑150 (MQ=255)
                       aCGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTc                                                                                                                                                      >  1:1051320/1‑128 (MQ=255)
                        cggccTGCACCGCTTCACGGAACATCGCCAGGCGGTCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGa                                                                                                                              <  2:941966/147‑1 (MQ=255)
                        cGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGa                                                                                                                              <  2:2454623/150‑1 (MQ=255)
                          cccTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGata                                                                                                                            >  2:336993/1‑150 (MQ=255)
                           ccTTCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGatat                                                                                                                           >  1:1645441/1‑150 (MQ=255)
                           ccTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGatat                                                                                                                           >  1:2713471/1‑150 (MQ=255)
                           ccTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATCACCAGTGCCACTTCCAGCGCGCTGCGCAAGTGGCCGatat                                                                                                                           >  1:160843/1‑150 (MQ=255)
                                aCCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACg                                                                                                                      <  2:1703695/150‑1 (MQ=255)
                                aCCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACg                                                                                                                      >  1:2311564/1‑150 (MQ=255)
                                        aCGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGaa                                                                                                              <  2:1195840/150‑1 (MQ=255)
                                        aCGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGaa                                                                                                              <  1:2436154/150‑1 (MQ=255)
                                         cGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGt                                                                                                                  >  2:1314615/1‑145 (MQ=255)
                                         cGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGt                                                                                                                  <  1:1314615/145‑1 (MQ=255)
                                         cGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAg                                                                                                             <  1:1160783/150‑1 (MQ=255)
                                          ggACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAgg                                                                                                            <  1:2038311/150‑1 (MQ=255)
                                          ggACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAgg                                                                                                            >  2:812/1‑150 (MQ=255)
                                          ggACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAAGAAGATCGCATCCAGCCCCATGAACTACGGCAACAACAGCGCGCTGCGCAAGTGGCCGATCTGCACGGGGACGAAgg                                                                                                            >  1:2981633/1‑150 (MQ=255)
                                           gACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGt                                                                                                           <  2:1645441/150‑1 (MQ=255)
                                              cATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTGCGCAGCGCGCTGGAAGTGGCCGag                                                                                                                             <  2:295493/129‑2 (MQ=38)
                                               aTCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTacca                                                                                                       >  2:1956164/1‑150 (MQ=255)
                                                tCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTaccac                                                                                                      >  2:695541/1‑150 (MQ=255)
                                                 cGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTaccacc                                                                                                     <  2:803781/150‑1 (MQ=255)
                                                  gCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGACGGTACCACCt                                                                                                    >  1:1962187/1‑150 (MQ=255)
                                                      ggCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGaa                                                                                                <  2:1800891/150‑1 (MQ=255)
                                                       gCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTCGCCCCCCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAAt                                                                                               <  1:2493177/150‑1 (MQ=255)
                                                        cGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATg                                                                                              <  1:812/150‑1 (MQ=255)
                                                         gATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGGTATGCACGGGGTCGAAGGTACCACCTAGAATGc                                                                                             >  1:1527118/1‑150 (MQ=255)
                                                          aTCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTc                                                                                                                 >  1:823147/1‑129 (MQ=255)
                                                          aTCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGcc                                                                                            <  1:1976569/150‑1 (MQ=255)
                                                          aTCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGcc                                                                                            <  2:1670971/150‑1 (MQ=255)
                                                          aTCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGcc                                                                                            <  2:2507584/150‑1 (MQ=255)
                                                              tGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATa                                                                                        >  2:853142/1‑150 (MQ=255)
                                                               gcgcGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATAc                                                                                       <  2:481622/150‑1 (MQ=255)
                                                                     ccACCTGTGGCGCGTCGCGGTGAGGCGACCGGGCATTGGGCACCAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGa                                                                                 <  2:160843/150‑1 (MQ=255)
                                                                      cACCTGTGGCGTGTCGCGGTGCGGAGGCCGGGCATTGGGAAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCACCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTCGAATGCCGATACGCCGGAc                                                                                >  1:416262/1‑150 (MQ=255)
                                                                       aCCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTggac                                                                                                                                >  1:968062/1‑99 (MQ=255)
                                                                       aCCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACt                                                                               >  1:169792/1‑150 (MQ=255)
                                                                       aCCTGGGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCAATTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACt                                                                               <  2:416262/150‑1 (MQ=255)
                                                                         cTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATCGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGc                                                                             >  2:2997239/1‑150 (MQ=255)
                                                                            tgGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCGCATCCAGCCCCATGACCTCGGCCACTTCCAGCGCGCTGCGCAAGTCGCCGCTCGGGACGTGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTg                                                                          <  1:2382236/150‑1 (MQ=255)
                                                                             ggCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCGACGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTgg                                                                         <  1:2338925/150‑1 (MQ=255)
                                                                             ggCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTgg                                                                         <  1:2844224/150‑1 (MQ=255)
                                                                              gCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTg                                                                                                                                   >  1:627717/1‑91 (MQ=255)
                                                                              gCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTg                                                                                                                                   <  2:627717/91‑1 (MQ=255)
                                                                              gCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGcc                                                                            <  2:331000/146‑1 (MQ=255)
                                                                              gCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGcc                                                                            >  1:331000/1‑146 (MQ=255)
                                                                              gCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTggg                                                                        >  2:2082832/1‑150 (MQ=255)
                                                                                   tCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGGCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTg                                                                   <  1:1614579/150‑1 (MQ=255)
                                                                                     gcgGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTac                                                                                                         <  2:1405696/110‑1 (MQ=255)
                                                                                     gcgGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTa                                                                                                   <  1:2678816/116‑1 (MQ=255)
                                                                                     gcgGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTa                                                                                                   >  2:2678816/1‑116 (MQ=255)
                                                                                              cggcCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTcagc                                                        <  2:1967501/150‑1 (MQ=255)
                                                                                              cggcCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTcagc                                                        <  1:610360/150‑1 (MQ=255)
                                                                                               ggcCGGGAATTGGGCAACAGGCGTCGCTCATCCCGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTcagca                                                       <  2:1697559/150‑1 (MQ=255)
                                                                                                   gggCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGAc                                                   <  2:2713471/150‑1 (MQ=255)
                                                                                                   gggCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGAc                                                   <  2:2311564/150‑1 (MQ=255)
                                                                                                     gCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGAccc                                                 <  1:1956164/150‑1 (MQ=255)
                                                                                                       aTTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTGCGCAGCGCGCTgg                                                                                                                                          >  2:2423482/1‑60 (MQ=25)
                                                                                                       aTTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCt                                               <  1:1867753/150‑1 (MQ=255)
                                                                                                        ttGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTg                                              <  2:1238585/150‑1 (MQ=255)
                                                                                                           ggCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGcc                                           <  1:2997239/150‑1 (MQ=255)
                                                                                                            gCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCCCCTCCCGCGCGCCTCGCCAGGGGCCCATAATCCCGGGGTCGAAGGTACCCCCTAAAAAGCCCATACGCCGGACCTCCCGGGCCTTGCCCAACTCAGCCAGAACCCTTGccc                                          >  2:49959/1‑149 (MQ=255)
                                                                                                            gCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGc                                                                                             <  2:2819418/99‑1 (MQ=255)
                                                                                                            gCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGa                                                    >  1:1450699/1‑140 (MQ=255)
                                                                                                            gCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCcg                                          >  2:2738697/1‑150 (MQ=255)
                                                                                                            gCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCcg                                          >  1:1462165/1‑150 (MQ=255)
                                                                                                              aaCAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATTCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCcgcg                                        >  1:354624/1‑150 (MQ=255)
                                                                                                              aaCAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGAccc                                                 >  1:1701645/1‑141 (MQ=255)
                                                                                                              aaCAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGAccc                                                 <  2:1701645/141‑1 (MQ=255)
                                                                                                              aaCAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCcgcg                                        >  1:2563210/1‑150 (MQ=255)
                                                                                                               aCAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCgcgc                                       >  2:55213/1‑150 (MQ=255)
                                                                                                                   gCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTgg                                                                         >  2:2126423/1‑112 (MQ=255)
                                                                                                                     gTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTgg                                 >  1:917315/1‑150 (MQ=255)
                                                                                                                       aGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGcc                               >  1:463192/1‑150 (MQ=255)
                                                                                                                       aGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGACTGATTTCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGcc                               <  1:1753783/150‑1 (MQ=255)
                                                                                                                         cTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGTACCCCTGCCCGCGCAACTGGCCAt                             >  2:426355/1‑150 (MQ=255)
                                                                                                                            gtCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATCGGCCTTACTGCCTGGGCCTTGCCCAACTAAGCAGGACCCCTGGCCGCGCAACTGGCCATAGc                          >  2:2105674/2‑150 (MQ=255)
                                                                                                                               cAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGa                       <  1:49959/150‑1 (MQ=255)
                                                                                                                                aGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAAc                                                             >  1:2798519/1‑111 (MQ=255)
                                                                                                                                aGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAAc                                                             <  2:2798519/111‑1 (MQ=255)
                                                                                                                                aGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGAt                      >  1:717780/1‑150 (MQ=255)
                                                                                                                                    ccATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTaccac                                                                                                      <  2:2587067/66‑1 (MQ=255)
                                                                                                                                     cATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATcacca                 <  2:463192/150‑1 (MQ=255)
                                                                                                                                          tatCGGCCACTTGCGCAGCGCGCTGGAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATca                    <  1:2222966/140‑1 (MQ=38)
                                                                                                                                           cTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGAccc                                                 >  2:2930687/1‑112 (MQ=255)
                                                                                                                                           cTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGAGACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGAccc                                                 <  1:2930687/112‑1 (MQ=255)
                                                                                                                                               ggCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATc                     >  1:912986/1‑137 (MQ=255)
                                                                                                                                               ggCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATc                     <  2:912986/137‑1 (MQ=255)
                                                                                                                                                 ccACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGaa                                                                                                <  2:509545/60‑1 (MQ=255)
                                                                                                                                                 ccACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGaa                                                                                                >  1:509545/1‑60 (MQ=255)
                                                                                                                                                 ccACTGCCAGCGCGCTGCGCAAGTGGACGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATCACCACATACTTCTCg      <  1:822057/150‑1 (MQ=255)
                                                                                                                                                     ttCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATCACCACATACTTCTCGCAgg  >  1:2874646/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GGGCATCGACGCTCAGGCAAGCGACGCCCTGCACCGCTTCACGGACCATCGCCAGGCGATCCTGCGCGGCCACCTGTGGCGTGTCGCGGTGCGGCGGCCGGGCATTGGGCAACAGGCGTAGCTCATCCAGCCCCATGAACTCGGCCACTTCCAGCGCGCTGCGCAAGTGGCCGATATGCACGGGGTCGAAGGTACCACCTAGAATGCCGATACGCCGGACTGCCTGGGCCTTGCTCAACTCAGCAGGACCCCTGGCCGCGCAACTGGCCATCGCCGATCACCACATACTTCTCGCAGG  >  NC_002947/5473332‑5473629

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