Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,912,816 A→T 11.8% L628Q (CTG→CAG)  ccmF ← holocytochrome c synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,912,8160AT11.8% 172.3 / 10.4 66L628Q (CTG→CAG) ccmFholocytochrome c synthetase
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base T (4/4);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAACACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCAC  >  NC_002947/4912673‑4912964
                                                                                                                                                |                                                                                                                                                    
agGAACACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCg                                                                                                                                                 >  5:86039/1‑149 (MQ=255)
   aaCACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCACGCAACCCGCCCAGGGCGGTCTGCTGAcc                                                                                                                                                  <  1:202042/146‑1 (MQ=255)
   aaCACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGGCGGTCTGCTGAcc                                                                                                                                                  >  2:202042/1‑145 (MQ=255)
       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                                                                                                                                     <  3:128377/143‑1 (MQ=255)
                  aCTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                                     >  4:128377/1‑143 (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)
                                   cTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTcagca                                                                                                                                                      >  4:144780/1‑109 (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                                                                        >  2:325430/1‑148 (MQ=255)
                                                                          aCCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGc                                                                        >  8:114634/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                                                 <  5:47777/149‑1 (MQ=255)
                                                                                                ggTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                                 <  1:21335/149‑1 (MQ=255)
                                                                                                ggTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGATCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                                 <  7:36006/149‑1 (MQ=255)
                                                                                                 gTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTg                                                                              >  7:86759/1‑119 (MQ=255)
                                                                                                 gTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTg                                                                              <  8:86759/119‑1 (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:247658/135‑1 (MQ=255)
                                                                                                      gcCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                         <  1:252734/135‑1 (MQ=255)
                                                                                                      gcCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                         >  2:247658/1‑135 (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                                                  >  6:121579/1‑137 (MQ=255)
                                                                                                           tCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTcgcg                                                  <  5:121579/137‑1 (MQ=255)
                                                                                                               gggCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGt                                    >  2:349646/1‑147 (MQ=255)
                                                                                                               gggCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGt                                    <  1:349646/147‑1 (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:97566/5‑147 (MQ=255)
                                                                                                                                    tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc                >  2:92948/5‑147 (MQ=255)
                                                                                                                                    tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc                <  1:97566/143‑1 (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             >  1:346812/1‑139 (MQ=255)
                                                                                                                                               cagACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCatcatc             <  2:346812/139‑1 (MQ=255)
                                                                                                                                                agACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCAc  >  8:269642/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGGAACACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCAGCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGACTGTTGCAC  >  NC_002947/4912673‑4912964

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: