Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,912,822 C→G 7.1% G626A (GGC→GCC)  ccmF ← holocytochrome c synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,912,8220CG 7.1% 224.8 / 5.9 70G626A (GGC→GCC) ccmFholocytochrome c synthetase
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base G (3/2);  total (38/32)
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.69e-01

CTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCACGGTGT  >  NC_002947/4912680‑4912969
                                                                                                                                              |                                                                                                                                                   
cTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCccagc                                                                                                                                                >  2:21335/1‑148 (MQ=255)
 tGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCccagc                                                                                                                                                >  8:118432/1‑147 (MQ=255)
 tGCAAGGGGCACTACCAAGATCCAACGCTTCATGCAGCTGCTCCAGACACGCACAGGACATCCCGCACCCGGGTCTTGACCTAGCCGCGGTAGCGCCGGTCGAGGGCCGCCAGCCACCCGCCCAGGCCGGTCCGCAGACCGCccagcca                                                                                                                                              >  1:350695/1‑149 (MQ=255)
    aaGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCccagc                                                                                                                                                >  2:342019/1‑144 (MQ=255)
      ggggCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                                          >  6:222287/1‑148 (MQ=255)
           aCTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                                          >  4:128377/1‑143 (MQ=255)
           aCTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                                          <  3:128377/143‑1 (MQ=255)
              aCCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACg                                                                                                                                 >  1:126674/1‑149 (MQ=255)
               ccATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACg                                                                                                                                 <  1:180074/148‑1 (MQ=255)
               ccATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGt                                                                                                                                <  6:86039/149‑1 (MQ=255)
               ccATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGt                                                                                                                                <  2:222942/149‑1 (MQ=255)
                                 tGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTAcgc                                                                                                              >  4:341947/1‑149 (MQ=255)
                                        gctCCAGACACGCCTAGGACATCAGGCCCCCGGGGCTTGACCTTTACGCGGTATCGCCGGTCGAGGGCCGCCAGCAACCCCCCCAGCCGGGGCAGCTGACCGCCCAGCCAGATCCACCGGACGTAAGGCTTAATATGTACGCGCACAGc                                                                                                       <  2:350695/149‑1 (MQ=255)
                                         ctCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTAc                                                                                                                >  6:220030/1‑139 (MQ=255)
                                         ctCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTAc                                                                                                                <  5:220030/139‑1 (MQ=255)
                                            cAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAg                                                                                                   >  3:2750/1‑149 (MQ=255)
                                               acacGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGgcg                                                                                                >  1:72016/1‑149 (MQ=255)
                                                cacGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGgcg                                                                                                <  4:2750/148‑1 (MQ=255)
                                                 acGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGgcgc                                                                                               >  7:206786/1‑148 (MQ=255)
                                                   gCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACCGCCCAGGCGCCAt                                                                                            >  6:9716/1‑149 (MQ=255)
                                                       aGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATtctc                                                                                        >  3:277235/1‑149 (MQ=255)
                                                          acctCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGGCAGCAGACCGCCCAGCCCGATCCAGCGGACGTAAGGCTTAGTATGTACGCGCACAGCCCAGGCGCCATTCTCCAg                                                                                     <  8:238831/146‑1 (MQ=255)
                                                             tCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTgg                                                                                  >  1:41449/1‑149 (MQ=255)
                                                              cACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCTCCCAGCCAGATc                                                                                                                                          <  3:3187/92‑1 (MQ=255)
                                                              cACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                                          >  4:3187/1‑92 (MQ=255)
                                                                cGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTc                                                                               >  6:341017/1‑149 (MQ=255)
                                                                cGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTc                                                                               <  5:222287/149‑1 (MQ=255)
                                                                  cACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGc                                                                             <  1:342019/149‑1 (MQ=255)
                                                                   aCCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGc                                                                             >  8:114634/1‑148 (MQ=255)
                                                                   aCCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGc                                                                             >  2:325430/1‑148 (MQ=255)
                                                                      cGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAg                                                                         <  2:72016/149‑1 (MQ=255)
                                                                             tgaccttgacGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGATCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACg                                                                  >  8:36006/1‑149 (MQ=255)
                                                                                ccttgacGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTAc                                                               <  8:206786/149‑1 (MQ=255)
                                                                                         ggTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                                      <  1:21335/149‑1 (MQ=255)
                                                                                         ggTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                                      <  5:47777/149‑1 (MQ=255)
                                                                                         ggTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGATCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                                      <  7:36006/149‑1 (MQ=255)
                                                                                          gTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTg                                                                                   <  8:86759/119‑1 (MQ=255)
                                                                                          gTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTg                                                                                   >  7:86759/1‑119 (MQ=255)
                                                                                           tAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGt                                                    <  1:325430/149‑1 (MQ=255)
                                                                                             gcgcCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACCTGTACGTcgcgcttg                                                   <  2:126674/148‑5 (MQ=255)
                                                                                               gcCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGCCGTACa                                                              >  2:252734/1‑135 (MQ=255)
                                                                                               gcCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                              <  1:252734/135‑1 (MQ=255)
                                                                                               gcCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                              >  2:247658/1‑135 (MQ=255)
                                                                                               gcCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                              <  1:247658/135‑1 (MQ=255)
                                                                                                ccGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGc                                               >  5:250674/1‑149 (MQ=255)
                                                                                                 cGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGcc                                              <  3:341947/149‑1 (MQ=255)
                                                                                                   gTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAgc                                                                        <  4:17546/121‑1 (MQ=255)
                                                                                                   gTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAgc                                                                        >  3:17546/1‑121 (MQ=255)
                                                                                                    tCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTcgcg                                                       <  5:121579/137‑1 (MQ=255)
                                                                                                    tCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTcgcg                                                       >  6:121579/1‑137 (MQ=255)
                                                                                                        gggCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGt                                         <  1:349646/147‑1 (MQ=255)
                                                                                                        gggCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGt                                         >  2:349646/1‑147 (MQ=255)
                                                                                                         ggCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGa                                      >  7:197769/1‑149 (MQ=255)
                                                                                                              ccAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGccc                                 <  7:114634/149‑1 (MQ=255)
                                                                                                               cAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCg                                >  3:84749/1‑149 (MQ=255)
                                                                                                                    aCCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCg                           >  6:206114/1‑149 (MQ=255)
                                                                                                                        gCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACCGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtca                       >  7:251948/1‑149 (MQ=255)
                                                                                                                           ctgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatca                    >  8:43734/6‑149 (MQ=255)
                                                                                                                            tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc                     >  2:92948/5‑147 (MQ=255)
                                                                                                                            tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc                     <  1:97566/143‑1 (MQ=255)
                                                                                                                            tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc                     >  2:97566/5‑147 (MQ=255)
                                                                                                                            tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc                     <  1:92948/143‑1 (MQ=255)
                                                                                                                                cGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGAc               <  4:277235/149‑1 (MQ=255)
                                                                                                                                      gcagACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTg          >  4:104128/1‑148 (MQ=255)
                                                                                                                                      gcagACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGc         >  4:99196/1‑149 (MQ=255)
                                                                                                                                       cagACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCatcatc                  <  2:346812/139‑1 (MQ=255)
                                                                                                                                       cagACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCatcatc                  >  1:346812/1‑139 (MQ=255)
                                                                                                                                        agACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCAc       >  8:269642/1‑149 (MQ=255)
                                                                                                                                             gCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCACGgtgt  <  4:84749/149‑1 (MQ=255)
                                                                                                                                             gCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCACGgtgt  <  7:118432/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCACGGTGT  >  NC_002947/4912680‑4912969

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