Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 165,929 G→T 7.1% A326A (GCC→GCA)  pntAA ← pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947165,9290GT 7.1% 238.1 / 4.5 77A326A (GCC→GCA) pntAApyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base T (4/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CTCCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCTG  >  NC_002947/165781‑166075
                                                                                                                                                    |                                                                                                                                                   
cTCCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAg                                                                                                                                                     <  3:166552/149‑1 (MQ=255)
 tCCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGACCAGCAGGTTGCGTCCGTACAGTCCCGAgg                                                                                                                                                    >  7:129867/1‑149 (MQ=255)
   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:49843/1‑149 (MQ=255)
                  ggACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTgg                                                                                                                                   >  3:115854/1‑149 (MQ=255)
                     cgacCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCTAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGc                                                                                                                                 <  8:136687/148‑1 (MQ=255)
                        ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                              >  8:47296/1‑148 (MQ=255)
                        ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                              >  7:19616/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                                                                                       <  8:327015/138‑1 (MQ=255)
                                                                        ggTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCt                                                                                       >  7:327015/1‑138 (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)
                                                                                                                                                    |                                                                                                                                                   
CTCCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAC‑GCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCTG  >  NC_002947/165781‑166075

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: ‑