Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,591 C→A 6.0% A99D (GCC→GAC)  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,5910CA 6.0% 217.8 / 3.6 67A99D (GCC→GAC) petAubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGGTAG  >  NC_002947/1504448‑1504737
                                                                                                                                               |                                                                                                                                                  
ggTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCgg                                                                                                                                               <  6:159973/149‑1 (MQ=255)
    aaCATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTAca                                                                                                                                           <  2:147800/149‑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                                                                                                                                  >  6:249850/1‑138 (MQ=255)
                        agggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTa                                                                                                                                  <  5:249850/138‑1 (MQ=255)
                        agggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTg                                                                                                                       <  2:135494/149‑1 (MQ=255)
                          ggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGAc                                                                                                                           <  7:176006/143‑1 (MQ=255)
                          ggcagCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGAc                                                                                                                           >  8:176006/1‑143 (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                                                                                        >  2:12211/1‑149 (MQ=255)
                                                       cAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAattctcat                                                                                        <  5:46512/149‑1 (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:63772/1‑149 (MQ=255)
                                                                           gcgAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTgcacgc                                                                    <  7:86816/149‑1 (MQ=255)
                                                                               aCGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCtgt                                                                          <  4:312640/139‑1 (MQ=255)
                                                                               aCGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCtgt                                                                          >  3:312640/1‑139 (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                               >  5:362494/1‑130 (MQ=255)
                                                                                                                                   ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                               <  6:362494/130‑1 (MQ=255)
                                                                                                                                   ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                               <  4:326743/130‑1 (MQ=255)
                                                                                                                                   ccGAGGCCTTGGACTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                               >  3:326743/1‑130 (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)
                                                                                                                                             ggCCTCGGTACAGGCTGCGTACGTGGACCCTGAGGTTCGCTCTATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCGGCCGACCTTGGGCCGAAATGGGTAg  <  4:95078/149‑1 (MQ=255)
                                                                                                                                             ggCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGGTAg  <  5:294687/149‑1 (MQ=255)
                                                                                                                                              gCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGgc                               <  8:4908/119‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGTCAACATCGCCAAGGTCGAGGCAGGGCAGCAGATGGTGGCCGAGTGGCGTGGCCAGCCCGTATTCATCGTGCGGCGAACGGAGGAAATCCTTGGCAATCTGAAAAAAATCACGGGTGACCTGTCCGACCCCGAGTCCAAGGCCTCGGTACAGCCAGCGTACGTCGACCCTGAGGTTCGCTCGATCAAGCCGGAAATTCTCATTCTCGTCGGCCTGTGCACGCACCTGGGCTGCTCGCCCACCTTCCGGCCAGAAGTGGCGCCCGCCGACCTTGGGCCGAAATGGGTAG  >  NC_002947/1504448‑1504737

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