Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,587 A→T 6.1% K98* (AAG→TAG)  petA → ubiquinol‑cytochrome c reductase iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,504,5870AT 6.1% 198.9 / 2.5 66K98* (AAG→TAG) petAubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base T (2/2);  total (35/31)
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.09e-01

GGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGG  >  NC_002947/1504448‑1504734
                                                                                                                                           |                                                                                                                                                   
ggTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCgg                                                                                                                                            <  6:159973/149‑1 (MQ=255)
    aaCATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAca                                                                                                                                        <  2:147800/149‑1 (MQ=255)
            cAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCCGTCCGACCCCGAGTCCAAg                                                                                                                                                   <  3:148318/130‑1 (MQ=255)
            cAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCCGTCCGACCCCGAGTCCAAg                                                                                                                                                   >  4:148318/1‑130 (MQ=255)
                 tCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTc                                                                                                                           <  5:134180/149‑1 (MQ=255)
                      gcagggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGAccc                                                                                                                      >  4:172299/1‑149 (MQ=255)
                      gcagggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGAcc                                                                                                                       <  8:142244/148‑1 (MQ=255)
                        agggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTa                                                                                                                               <  5:249850/138‑1 (MQ=255)
                        agggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTa                                                                                                                               >  6:249850/1‑138 (MQ=255)
                        agggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTg                                                                                                                    <  2:135494/149‑1 (MQ=255)
                          ggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGAc                                                                                                                        >  8:176006/1‑143 (MQ=255)
                          ggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGAc                                                                                                                        <  7:176006/143‑1 (MQ=255)
                          ggcagCAGATGGTGGCCGAGTGGCGGGGCCAGCCCGTATTCAACGAGCGGCGAACGGAGGAATTCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGCCCGACCCCGACGCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAg                                                                                                                  <  6:182430/149‑1 (MQ=255)
                               cagATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCg                                                                                                             >  5:163332/1‑149 (MQ=255)
                                      tggCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAcagcca                                                                                                                                    <  5:117062/119‑1 (MQ=255)
                                      tggCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAcagcca                                                                                                                                    >  6:117062/1‑119 (MQ=255)
                                      tggCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATc                                                                                                      <  6:163332/149‑1 (MQ=255)
                                         ccGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTc                                                                                                          >  6:120373/1‑142 (MQ=255)
                                         ccGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTc                                                                                                          <  5:120373/142‑1 (MQ=255)
                                          cGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCaa                                                                                                    <  3:286836/147‑1 (MQ=255)
                                          cGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCaa                                                                                                    >  4:286836/1‑147 (MQ=255)
                                           gAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGcc                                                                                                 >  6:178547/1‑149 (MQ=255)
                                                gcgtggcCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCAAGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGaaa                                                                                            >  3:308502/1‑149 (MQ=255)
                                                  gtggcCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAat                                                                                           >  6:46512/1‑148 (MQ=255)
                                                      ccAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAattctca                                                                                      >  8:10422/1‑149 (MQ=255)
                                                      ccAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACATCGACCCTGAGGTTCGCTCGATCAAGCCGGAAattctca                                                                                      <  2:197578/149‑1 (MQ=255)
                                                       cAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAattctcat                                                                                     <  5:46512/149‑1 (MQ=255)
                                                       cAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAattctcat                                                                                     >  2:12211/1‑149 (MQ=255)
                                                          cgcGTATTCTTCGTGCGGCGAACGGAGCAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAattctcattc                                                                                   <  8:57268/146‑1 (MQ=255)
                                                               aTTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCtcgtcg                                                                             >  4:335524/1‑149 (MQ=255)
                                                                   aTCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCt                                                                         >  8:78976/1‑149 (MQ=255)
                                                                   aTCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCt                                                                         >  3:95078/1‑149 (MQ=255)
                                                                      gTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCtgtg                                                                      >  2:71637/1‑149 (MQ=255)
                                                                          ggcgAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTgcacg                                                                  >  3:108835/1‑149 (MQ=255)
                                                                           gcgAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTgcacgc                                                                 <  7:86816/149‑1 (MQ=255)
                                                                           gcgAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTgcacgc                                                                 >  7:63772/1‑149 (MQ=255)
                                                                               aCGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCtgt                                                                       >  3:312640/1‑139 (MQ=255)
                                                                               aCGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCtgt                                                                       <  4:312640/139‑1 (MQ=255)
                                                                               aCGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCt                                                             <  5:178547/149‑1 (MQ=255)
                                                                                  gaggaAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTggg                                                          >  6:137256/1‑149 (MQ=255)
                                                                                   aggaAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGgc                                                         >  4:299013/1‑149 (MQ=255)
                                                                                     gaAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGgctg                                                       <  3:335524/149‑1 (MQ=255)
                                                                                        aTCCTTGGCAATCTGAAAAAAATCTCGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTc                                                    <  7:10422/149‑1 (MQ=255)
                                                                                           cTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGc                                                  <  5:137256/148‑1 (MQ=255)
                                                                                               gCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCAcc                                             >  4:333139/1‑149 (MQ=255)
                                                                                                     tGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCgg                                       >  1:34596/1‑149 (MQ=255)
                                                                                                       aaaaaaaTCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCg                                        <  4:108835/146‑1 (MQ=255)
                                                                                                       aaaaaaaTCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGc                                      >  6:164519/1‑148 (MQ=255)
                                                                                                          aaaaTCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAg                                   >  3:239719/1‑148 (MQ=255)
                                                                                                                aCGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTgg                             <  3:172299/148‑1 (MQ=255)
                                                                                                                     tGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCtcgt                                                                               >  3:93420/1‑93 (MQ=255)
                                                                                                                     tGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCtcgt                                                                               <  4:93420/93‑1 (MQ=255)
                                                                                                                           gTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGAcc                 >  5:193118/1‑149 (MQ=255)
                                                                                                                           gTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGACGGAAATTCTCATTCTCGTCGGACTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCAAGAAGTGGCGCCCGCCGAcc                 >  1:245769/1‑149 (MQ=255)
                                                                                                                            tCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCt                >  6:294687/1‑149 (MQ=255)
                                                                                                                            tCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCGCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGAc                  <  1:71637/147‑1 (MQ=255)
                                                                                                                              cGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTg              >  5:258268/1‑149 (MQ=255)
                                                                                                                                aCCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTggg            <  5:164519/149‑1 (MQ=255)
                                                                                                                                   ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                            >  3:326743/1‑130 (MQ=255)
                                                                                                                                   ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                            <  6:362494/130‑1 (MQ=255)
                                                                                                                                   ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                            >  5:362494/1‑130 (MQ=255)
                                                                                                                                   ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                            <  4:326743/130‑1 (MQ=255)
                                                                                                                                       gTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCAGCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAAt     >  5:313164/1‑149 (MQ=255)
                                                                                                                                       gTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAAt     <  4:239719/149‑1 (MQ=255)
                                                                                                                                          cAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATgg   <  3:299013/148‑1 (MQ=255)
                                                                                                                                          cAAGCCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAAATCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCACAAGTGGCGCCCGCCGACCTTGGTCCGAAATggg  >  6:146831/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGG  >  NC_002947/1504448‑1504734

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