Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,450,173 C→G 11.4% A68G (GCC→GGC)  nicA → nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,450,1730CG11.4% 265.6 / 24.1 97A68G (GCC→GGC) nicAnicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base C (39/47);  new base G (6/5);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

Only 100 of 103 total aligned reads displayed.

ATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCA  >  NC_002947/4450025‑4450322
                                                                                                                                                    |                                                                                                                                                     
aTGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTCTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGgcgcg                                                                                                                                                      <  2:1649859/150‑2 (MQ=255)
    cTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGccgccg                                                                                                                                                  <  2:2028329/150‑1 (MQ=255)
         acCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaa                                                                                                                                             <  1:1294306/150‑3 (MQ=255)
         acCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCCGCGTCATCCCCTTGGCCGGCGGCgccggccaa                                                                                                                                             <  1:993226/150‑3 (MQ=255)
          ccccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                            <  2:1847497/150‑4 (MQ=255)
          ccccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTCGGCCGGCGGCgccggccaag                                                                                                                                            <  1:1218890/150‑4 (MQ=255)
           cccGCTGCTGCTGATCCTGCGTAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCaa                                                                                                                                           >  1:1557503/1‑150 (MQ=255)
              gctgctgctGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAatat                                                                                                                                        >  1:438190/1‑150 (MQ=255)
                tgctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATca                                                                                                                                      <  2:1751508/150‑1 (MQ=255)
                     ctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCaccacc                                                                                                                                 >  2:1433502/1‑150 (MQ=255)
                      tgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCAccc                                                                                                                                >  2:145851/1‑150 (MQ=255)
                         tCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCg                                                                                                                             >  1:921954/1‑150 (MQ=255)
                           cTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                           >  1:2380917/1‑150 (MQ=255)
                            tGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAg                                                                                                                          >  1:1124379/1‑150 (MQ=255)
                              cgcAACGACCTGTGCCTGAACGTCCCCAAGTACGGCTGCGCCCTGGCCGAATGCGGTGCTTCTACGGTGATCATCGACTGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAggg                                                                                                                        <  2:1196920/150‑1 (MQ=255)
                                 aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                            >  1:919960/1‑124 (MQ=255)
                                 aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                            <  2:919960/127‑4 (MQ=255)
                                 aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCAGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                     <  1:444115/150‑1 (MQ=255)
                                     aCCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                 <  2:1899489/150‑1 (MQ=255)
                                     aCCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                 <  1:1767461/150‑1 (MQ=255)
                                      ccTGTGCCTGAACGGCCCCAAGTACGGCTTCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcag                                                                                                                >  1:1874184/1‑150 (MQ=255)
                                       cTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccgg                                                                                                                                                 >  2:2235045/1‑116 (MQ=255)
                                       cTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccgg                                                                                                                                                 <  1:2235045/116‑1 (MQ=255)
                                        tgtgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagca                                                                                                              >  2:43761/1‑150 (MQ=255)
                                        tgtgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagca                                                                                                              >  2:1383248/1‑150 (MQ=255)
                                         gtgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCaa                                                                                                             <  2:1630516/150‑1 (MQ=255)
                                              tGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccg                                                                                                        >  2:1635279/1‑150 (MQ=255)
                                                aaCGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAat                                                                                                                                          >  2:2409092/1‑114 (MQ=255)
                                                       ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCGCCACCCTCGAAGGGCTGGGCAGCAAGGCCgcaccgcacc                                                                                               >  1:22854/1‑150 (MQ=255)
                                                       ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                    <  2:2307112/129‑1 (MQ=255)
                                                       ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                    >  1:2307390/1‑129 (MQ=255)
                                                       ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                    >  1:2307112/1‑129 (MQ=255)
                                                       ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                    <  2:2307390/129‑1 (MQ=255)
                                                        cAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCAccc                                                                                              >  1:1074306/1‑150 (MQ=255)
                                                         aaGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                     <  1:381454/126‑1 (MQ=255)
                                                         aaGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                     >  2:381454/1‑126 (MQ=255)
                                                           ttACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGGCGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCTGCGCCGCCGGCCGCAATATCACCTCCCGCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGt                                                                                           <  1:797130/149‑1 (MQ=255)
                                                           gTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATCTCAACAACCACGAAGAGCTGGGCAGCAAGGCCGCACCGCACCCCGt                                                                                           >  2:2323151/1‑150 (MQ=255)
                                                           gTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                 <  1:1813647/128‑1 (MQ=255)
                                                           gTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                 >  2:1813647/1‑128 (MQ=255)
                                                                  tGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAg                                                                                    <  2:1124379/150‑1 (MQ=255)
                                                                   gCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGGGGCCGCCCGCTCCTGGATCATCCCCTTGTCCTGCGCCTCCGTCCGCAATATCTCCCTCCTCGTAGCGCTGGGCAGCAAGGCCGCAACGCACCCCGTGTAACAgg                                                                                   >  2:1025064/1‑150 (MQ=255)
                                                                    cGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCCCCACCGCCGGACGCAATTGCACCGCCCTTTTAATGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGc                                                                                  <  1:45747/150‑1 (MQ=255)
                                                                      gCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATGCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGt                                                                                <  1:168186/150‑1 (MQ=255)
                                                                      gCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGt                                                                                <  1:556735/150‑1 (MQ=255)
                                                                      gCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGt                                                                                <  2:921954/150‑1 (MQ=255)
                                                                         tGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTg                                                                                          <  1:869653/137‑1 (MQ=255)
                                                                         tGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTg                                                                                          >  2:869653/1‑137 (MQ=255)
                                                                          gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCTGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCAt                                                                            <  1:1760049/150‑1 (MQ=255)
                                                                          gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCAt                                                                            <  1:190796/150‑1 (MQ=255)
                                                                          gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCAt                                                                            <  2:1557503/150‑1 (MQ=255)
                                                                          gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGCACGGCTCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCAt                                                                            <  1:1731532/150‑1 (MQ=255)
                                                                            gCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAg                                                                                                            <  1:1379554/116‑1 (MQ=255)
                                                                            gCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAg                                                                                                            >  2:1379554/1‑116 (MQ=255)
                                                                             cGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCaa                                                                                                             >  2:211851/1‑114 (MQ=255)
                                                                             cGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCaa                                                                                                             <  1:211851/114‑1 (MQ=255)
                                                                              gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGcgccgc                                                                                                                                                    <  2:25928/74‑1 (MQ=255)
                                                                              gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGcgccgc                                                                                                                                                    >  1:25928/1‑74 (MQ=255)
                                                                              gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgac                                                                        <  1:2309479/150‑1 (MQ=255)
                                                                              gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgac                                                                        >  1:1041005/1‑150 (MQ=255)
                                                                              gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgac                                                                        <  2:1874184/150‑1 (MQ=255)
                                                                                aTGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgacga                                                                      <  2:163157/150‑1 (MQ=255)
                                                                                  gCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGc                                                                    <  2:462649/150‑1 (MQ=255)
                                                                                         gtgtACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAggcggc                                                             >  1:73788/1‑150 (MQ=255)
                                                                                                 tGATCATCGACGGCGTGGCCGCTCGCTCCTGCGTCACCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCGCCCTCGTGGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAGGAGGCGTTCATCGACGAGCAGGCGGCCCAGtgtg                                                     <  2:1768402/150‑1 (MQ=255)
                                                                                                 tGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGtgtg                                                     <  1:1433502/150‑1 (MQ=255)
                                                                                                  gATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCCCCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTgg                                                    <  2:1782940/150‑1 (MQ=255)
                                                                                                    tcatcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGtt                                                  >  2:1529794/1‑150 (MQ=255)
                                                                                                      ctcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTAc                                                <  1:288056/149‑1 (MQ=255)
                                                                                                         gACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGAGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGc                                             >  2:1491646/1‑150 (MQ=255)
                                                                                                          aCGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGcc                                            <  2:1654231/150‑2 (MQ=255)
                                                                                                           cGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCAt                                           <  1:1635279/150‑1 (MQ=255)
                                                                                                               gTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAAc                                       <  2:966954/150‑1 (MQ=255)
                                                                                                               gTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAAc                                       <  1:2323151/150‑1 (MQ=255)
                                                                                                                 ggCCGCCCTCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACgg                                     <  1:1553411/150‑1 (MQ=255)
                                                                                                                 ggCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACgg                                     >  2:1736667/1‑150 (MQ=255)
                                                                                                                  gccgccCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGc                                    >  1:1897713/1‑150 (MQ=255)
                                                                                                                  gccgccCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGc                                    <  2:1038987/150‑1 (MQ=255)
                                                                                                                   ccgcccgcTCCTCCGTCATACCTCTAGCCGGTGCCGCCGGCCGCAATAACACCACCCTCCAAGGGTTTCACATCAAGGCCGCACCGCACCCCGTGCAACAGGCGTGCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCa                                   <  1:42809/150‑1 (MQ=255)
                                                                                                                        cgcTCCTGCGTCAGCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATc                              <  2:964747/150‑1 (MQ=255)
                                                                                                                            ccTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGa                          <  2:1897713/150‑1 (MQ=255)
                                                                                                                                    aTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAg                  >  2:571238/1‑150 (MQ=255)
                                                                                                                                    aTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACAGcacag                  >  1:2042882/1‑146 (MQ=255)
                                                                                                                                       attgcggCCGGCGGCGCCGGCCACAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGcc               >  1:247607/6‑150 (MQ=255)
                                                                                                                                        ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCAt                            >  2:2021794/1‑136 (MQ=255)
                                                                                                                                        ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCAt                            <  1:2021794/136‑1 (MQ=255)
                                                                                                                                        ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGAcc                        >  1:884598/1‑140 (MQ=255)
                                                                                                                                        ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGAcc                        <  2:884598/140‑1 (MQ=255)
                                                                                                                                        ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGccc              >  2:707575/1‑150 (MQ=255)
                                                                                                                                         cTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCt             >  2:1589178/1‑150 (MQ=255)
                                                                                                                                          gcggCCGGCGGCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCCCGAAGGCCCTg            >  1:166164/3‑150 (MQ=255)
                                                                                                                                           cggCCGGCGGCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGc           >  1:1390010/2‑150 (MQ=255)
                                                                                                                                             gCCGGCGGCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTc         >  1:1433075/1‑150 (MQ=255)
                                                                                                                                             gCCGGCGGCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTc         >  2:88471/1‑150 (MQ=255)
                                                                                                                                             gCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCCCCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTc         <  1:1383248/150‑1 (MQ=255)
                                                                                                                                              ccGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCg        >  2:2246502/1‑150 (MQ=255)
                                                                                                                                                ggCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGAt      <  1:145851/150‑1 (MQ=255)
                                                                                                                                                ggCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCAAGGCCCTGCTCGATc     <  1:1736667/150‑1 (MQ=255)
                                                                                                                                                 gcgcCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATc     <  2:2380917/150‑1 (MQ=255)
                                                                                                                                                    ccgccgGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCa  <  2:438190/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCA  >  NC_002947/4450025‑4450322

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: