Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 165,932 A→C 6.3% D325E (GAT→GAG)  pntAA ← pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947165,9320AC 6.3% 228.3 / 4.9 76D325E (GAT→GAG) pntAApyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base C (3/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCTGCGCCA  >  NC_002947/165784‑166080
                                                                                                                                                    |                                                                                                                                                     
cTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCa                                                                                                                                                       >  4:324272/1‑149 (MQ=255)
 ttAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCAt                                                                                                                                                      >  8:72472/1‑149 (MQ=255)
       gTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCAc                                                                                                                                                <  8:216065/149‑1 (MQ=255)
        tttttGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACt                                                                                                                                               <  8:236450/149‑1 (MQ=255)
        tttttGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCAAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACt                                                                                                                                               >  7:136687/1‑149 (MQ=255)
            tGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCAc                                                                                                                                            <  8:148064/148‑1 (MQ=255)
              cGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTg                                                                                                                                         <  6:399228/149‑1 (MQ=255)
               ggACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTgg                                                                                                                                        >  3:115854/1‑149 (MQ=255)
               ggACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTgg                                                                                                                                        >  3:49843/1‑149 (MQ=255)
                  cgacCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCTAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGc                                                                                                                                      <  8:136687/148‑1 (MQ=255)
                     ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                                   >  7:19616/1‑148 (MQ=255)
                     ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                                   >  8:47296/1‑148 (MQ=255)
                     ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTAGTCGAACAGCACCTTCATGAAGTCCAGCAGCTTGCGTGCGGACAGTGCCGATGCATCGCCACTCACCTGGGCCgc                                                                                                                                   >  8:329541/1‑147 (MQ=255)
                      cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                                   >  4:38281/1‑147 (MQ=255)
                         aCCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGt                                                                                                                              <  6:94463/149‑1 (MQ=255)
                           cATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTg                                                                                                                            >  2:184493/1‑149 (MQ=255)
                            aTCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTgg                                                                                                                           >  6:51700/1‑149 (MQ=255)
                               aCGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCg                                                                                                                        >  5:221345/1‑149 (MQ=255)
                                 ggCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCggg                                                                                                                      <  6:275059/149‑1 (MQ=255)
                                  gCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGc                                                                                                                     <  6:342596/149‑1 (MQ=255)
                                   cacaTCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGcc                                                                                                                    >  2:246515/1‑149 (MQ=255)
                                       tCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCcgac                                                                                                                 >  3:168823/1‑148 (MQ=255)
                                         aggcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgat                                                                                                              <  3:324272/149‑1 (MQ=255)
                                               gcCGCGACGATGTCGTCTTCGAGGTTGATGATCACCGCTCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAc                                                                                                        <  8:91026/149‑1 (MQ=255)
                                                ccGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑g                                                                                                      <  7:280222/149‑1 (MQ=255)
                                                 cgcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGGAGGTTGGTCGGGCCGACGATGATCAC‑Gc                                                                                                     <  1:220184/149‑1 (MQ=255)
                                                 cgcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑Gc                                                                                                     <  5:248529/149‑1 (MQ=255)
                                                  gctacgaAGTCGTCTTCGCGGTTCATGATCATCGCGCCATCCTTTTCGAACAGCAGCTTCATTAAGTCCAGCAGGTTGCTTGCGTACATTGCCTAGTAATAGGAACTCACCTGGTCCGGCAGTTTGGTCGGGCCGACGATGATCAC‑Gcc                                                                                                    <  2:153089/146‑1 (MQ=255)
                                                         tGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTcc                                                                                             <  5:51700/149‑1 (MQ=255)
                                                           tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTg                                                                                           <  2:120256/149‑1 (MQ=255)
                                                             gtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGc                                                                                          <  4:298449/148‑1 (MQ=255)
                                                              tcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGcac                                                                                        >  4:76272/1‑149 (MQ=255)
                                                                ttCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGcacca                                                                                      >  5:236656/1‑149 (MQ=255)
                                                                    aGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTg                                                                                  <  3:38281/149‑1 (MQ=255)
                                                                     ggTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCt                                                                                            >  7:327015/1‑138 (MQ=255)
                                                                     ggTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCt                                                                                            <  8:327015/138‑1 (MQ=255)
                                                                         gatgatCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCt                                                                              <  4:49843/148‑1 (MQ=255)
                                                                           tgatCAGCGCGCCATCCTTGTCGAACAGCAGCTTCAAGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGCTTGGTCGGCCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGCTCTGc                                                                            >  1:62478/1‑148 (MQ=255)
                                                                                    cgcCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGc                                                                  <  6:221345/149‑1 (MQ=255)
                                                                                       cATCCTTGTCGAACAGCAGCTTCATGAAGTGCAGCAGGTGGCGTGCGTACAGTGCCGAGGCTTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCggg                                                               <  7:329541/149‑1 (MQ=255)
                                                                                       cATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCggg                                                               <  7:47296/149‑1 (MQ=255)
                                                                                        aTCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTCGTCTGCGACGGTCAGCggg                                                               >  7:310130/1‑148 (MQ=255)
                                                                                         tCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCa                                                             <  8:214571/149‑1 (MQ=255)
                                                                                         tCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCa                                                             <  1:283212/149‑1 (MQ=255)
                                                                                           cTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACTGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCa                                                             <  1:76077/148‑1 (MQ=255)
                                                                                           cTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGc                                                                  >  1:202033/1‑142 (MQ=255)
                                                                                           cTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGc                                                                  <  2:202033/142‑1 (MQ=255)
                                                                                            ttGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGtt                                                          >  8:369416/1‑149 (MQ=255)
                                                                                               tCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTgcc                                                       >  3:112329/1‑149 (MQ=255)
                                                                                                     agcagcTTCATGAAGTCCAGCTGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTg                                                 <  8:129867/149‑1 (MQ=255)
                                                                                                     agcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTg                                                 >  6:125113/1‑149 (MQ=255)
                                                                                                     agcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTg                                                 <  1:184493/149‑1 (MQ=255)
                                                                                                      gcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                                 <  1:88602/69‑1 (MQ=255)
                                                                                                      gcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                                 >  2:88602/1‑69 (MQ=255)
                                                                                                       cagcTTCATGAAGTCCAGCAGGTTGCGTCCCTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATAAC‑GCCGTTCTCCTGCACCACCAGGTCTGCCACGGTCAGCGGGCAGTTGCCGCCTTggg                                               >  6:169598/1‑149 (MQ=255)
                                                                                                             cATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCa                                          <  7:369416/148‑1 (MQ=255)
                                                                                                             cATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCa                                          <  3:76272/148‑1 (MQ=255)
                                                                                                             cATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCa                                          <  4:168823/148‑1 (MQ=255)
                                                                                                             cATCAAGTCCTGCAGGTTGGGTGCGTACAGTGCCGTGGCTTCGGCAGTCACCTTGTCCGGCAGGTTGGTCGGTCCGTCGATGTTCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGAGGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCa                                          <  5:169598/148‑1 (MQ=255)
                                                                                                                 aaGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAg                                     >  3:50704/1‑149 (MQ=255)
                                                                                                                  aGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAgg                                    <  3:367583/149‑1 (MQ=255)
                                                                                                                       agcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGAt                               >  8:95096/1‑149 (MQ=255)
                                                                                                                       agcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGAt                               <  1:246515/149‑1 (MQ=255)
                                                                                                                       agcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGAt                               <  7:72472/149‑1 (MQ=255)
                                                                                                                       agcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGAt                               >  7:349170/1‑149 (MQ=255)
                                                                                                                        gcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGACCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGAt                               >  4:59870/1‑148 (MQ=255)
                                                                                                                                       aCAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGc               >  3:324516/1‑149 (MQ=255)
                                                                                                                                        gaGTGCCGATGCCTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCgg                 <  2:267101/145‑1 (MQ=255)
                                                                                                                                        gaGTGCCGATGCCTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCgg                 >  1:267101/2‑146 (MQ=255)
                                                                                                                                         aGTGCCGATGCCTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCt              >  4:260214/1‑148 (MQ=255)
                                                                                                                                         aGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCtt             >  7:179202/1‑149 (MQ=255)
                                                                                                                                         aGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCtt             >  1:51068/1‑149 (MQ=255)
                                                                                                                                          gTGCCGATGCCTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCtt             >  6:370704/1‑148 (MQ=255)
                                                                                                                                             ccGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATc         >  3:196852/1‑149 (MQ=255)
                                                                                                                                              cGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCt        >  6:62763/1‑149 (MQ=255)
                                                                                                                                              cGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCt        >  3:39200/1‑149 (MQ=255)
                                                                                                                                               gAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACAGAGCCGGGCTTCATCTg       >  6:77342/1‑149 (MQ=255)
                                                                                                                                                    aTCGGCACTCACATGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTATCCTGCACAACCTGGTCTGCGACCGTCCGCGGGCAGTTGCCGCCTTGGGCTGCCGCCAGGTCGATGACCACCGAGCAGGGCTTCATCTGCGCCa  >  2:86458/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCTGCGCCA  >  NC_002947/165784‑166080

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑