Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,584 T→G 6.7% S97A (TCC→GCC)  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,5840TG 6.7% 213.6 / 4.5 71S97A (TCC→GCC) petAubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (35/30);  new base G (2/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAAT  >  NC_002947/1504437‑1504731
                                                                                                                                                   |                                                                                                                                                   
gCACCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTc                                                                                                                                                    <  5:303768/149‑1 (MQ=255)
  aCCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTcc                                                                                                                                                   >  7:142244/1‑148 (MQ=255)
   ccGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCATCGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGAcc                                                                                                                                                   >  7:57268/1‑147 (MQ=255)
           ggTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCgg                                                                                                                                         <  6:159973/149‑1 (MQ=255)
               aaCATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAca                                                                                                                                     <  2:147800/149‑1 (MQ=255)
                       cAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCCGTCCGACCCCGAGTCCAAg                                                                                                                                                >  4:148318/1‑130 (MQ=255)
                       cAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCCGTCCGACCCCGAGTCCAAg                                                                                                                                                <  3:148318/130‑1 (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                                                                                                                                 >  6:117062/1‑119 (MQ=255)
                                                 tggCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAcagcca                                                                                                                                 <  5:117062/119‑1 (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)
                                                                          aTTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTcc                                                                                                                                                   >  7:43172/1‑76 (MQ=255)
                                                                          aTTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTcc                                                                                                                                                   <  2:326395/76‑1 (MQ=255)
                                                                          aTTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTcc                                                                                                                                                   <  8:43172/76‑1 (MQ=255)
                                                                          aTTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTcc                                                                                                                                                   >  1:326395/1‑76 (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:63772/1‑149 (MQ=255)
                                                                                      gcgAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTgcacgc                                                              <  7:86816/149‑1 (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                                                                            <  4:93420/93‑1 (MQ=255)
                                                                                                                                tGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCtcgt                                                                            >  3:93420/1‑93 (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                         <  4:326743/130‑1 (MQ=255)
                                                                                                                                              ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                         >  5:362494/1‑130 (MQ=255)
                                                                                                                                                  gTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCAGCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAAt  >  5:313164/1‑149 (MQ=255)
                                                                                                                                                  gTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAAt  <  4:239719/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCACCGGTGAAGGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAAT  >  NC_002947/1504437‑1504731

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