Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,914,788 C→A 5.1% L121L (CTG→CTT ccmE ← cytochrome c‑type biogenesis protein CcmE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,914,7880CA 5.1% 254.9 / 3.3 78L121L (CTG→CTTccmEcytochrome c‑type biogenesis protein CcmE
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAACGCACGTC  >  NC_002947/4914645‑4914934
                                                                                                                                               |                                                                                                                                                  
ggCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGAAc                                                                                                                                               >  8:307200/1‑149 (MQ=255)
  ccAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACtt                                                                                                                                             >  5:307849/1‑149 (MQ=255)
   cAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTGGCCAAGCAcgac                                                                                                                                            <  6:165020/149‑4 (MQ=255)
           tcGACAACTGGCACATTTAGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCAGCAGGCAGCGGACTGGCCGCTTTCCTTCAGGGCCGTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCAATTCATCGGCTa                                                                                                                                    >  2:235576/1‑149 (MQ=255)
             gCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTac                                                                                                                                   <  2:155735/148‑1 (MQ=255)
               cAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTaccac                                                                                                                                <  6:307849/149‑1 (MQ=255)
                   tGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGTATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATc                                                                                                                                         <  8:303553/136‑1 (MQ=255)
                   tGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGGATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATc                                                                                                                                         >  7:303553/1‑136 (MQ=255)
                     gCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACAcc                                                                                                                          >  4:174395/1‑149 (MQ=255)
                       ccAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGt                                                                                                                        >  1:34611/1‑149 (MQ=255)
                       ccAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTAGGGCGGCATGTATTTATAGTAGTGCTTGGCCAGCACATCAAAGGCTACCACCAAACCg                                                                                                                         >  2:156604/1‑148 (MQ=255)
                             cGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAt                                                                                                                  >  2:259855/1‑149 (MQ=255)
                              ggggACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                                                                 <  6:323760/149‑1 (MQ=255)
                              ggggACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                                                                 <  7:8037/149‑1 (MQ=255)
                               gggACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                                                                 >  1:257991/1‑148 (MQ=255)
                               gggACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                                                                 >  7:5112/1‑148 (MQ=255)
                                 gACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAg                                                                                                              <  4:150022/149‑1 (MQ=255)
                                      ccAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGc                                                                                                         >  5:302612/1‑149 (MQ=255)
                                      ccAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGc                                                                                                         >  8:306253/1‑149 (MQ=255)
                                       cAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGcc                                                                                                        <  3:174395/149‑1 (MQ=255)
                                            ctgcGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAggg                                                                                                   <  3:253645/149‑1 (MQ=255)
                                                gTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGcgac                                                                                               <  7:306253/149‑1 (MQ=255)
                                                gTTCATGGCTTCGTCTCCCCCCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCGCTTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGcgac                                                                                               <  1:235576/149‑1 (MQ=255)
                                                   cATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTAAATCTCGTCGTGCTTAGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAt                                                                                            >  1:278861/1‑149 (MQ=255)
                                                    aTGGCTTCGTCTCCCCTCCTGCAGCGGCCTTGCAGCGTTGCTTCAGGGCCTTGGTGTCTTCTGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCAACCACCACACCGTAGGCATTGAGCTTGCCAAGGGCGACGAtg                                                                                           <  4:45896/149‑1 (MQ=255)
                                                     tGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgc                                                                                          >  8:68855/1‑149 (MQ=255)
                                                      ggCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgcc                                                                                         >  5:30829/1‑149 (MQ=255)
                                                       gCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccc                                                                                        >  6:246061/1‑149 (MQ=255)
                                                       gCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgcc                                                                                         <  2:34611/148‑1 (MQ=255)
                                                        cTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccct                                                                                       <  6:30829/149‑1 (MQ=255)
                                                            gTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccctgccc                                                                                   >  8:241159/1‑149 (MQ=255)
                                                             tctcCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATgccctgccct                                                                                  >  2:245816/1‑149 (MQ=255)
                                                               tcCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTc                                                                                >  8:2512/1‑149 (MQ=255)
                                                                 cccGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTcg                                                                               <  8:208795/148‑1 (MQ=255)
                                                                 cccGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACGTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGGATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTgcgc                                                                              >  3:33038/1‑149 (MQ=255)
                                                                  ccgccgGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTcgcg                                                                             >  5:94340/1‑149 (MQ=255)
                                                                    aaggaaaGCGGCCAGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCgg                                                                            >  4:227267/7‑148 (MQ=255)
                                                                      ggaaaGCGGCCAGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAggg                                                                                                   >  6:128561/5‑123 (MQ=255)
                                                                      ggaaaGCGGCCAGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAggg                                                                                                   <  5:128561/119‑1 (MQ=255)
                                                                        gCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAg                                                                       <  8:5112/149‑1 (MQ=255)
                                                                              gacTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCggg                                                                 >  6:347624/3‑149 (MQ=255)
                                                                                 tGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggag                                                              >  3:329175/1‑149 (MQ=255)
                                                                                   gCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggagga                                                            <  1:259855/149‑1 (MQ=255)
                                                                                   gCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggagga                                                            <  6:302612/149‑1 (MQ=255)
                                                                                   gCCGCGTTCCTTCCGGGCCTTGGGGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggag                                                              <  4:33038/147‑1 (MQ=255)
                                                                                    ccGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGAt                                                           >  3:266358/1‑149 (MQ=255)
                                                                                        tttCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCcg                                                       >  1:75070/1‑149 (MQ=255)
                                                                                            cTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCgg                                                    <  7:307200/148‑1 (MQ=255)
                                                                                             ttCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCgg                                                    >  2:122163/1‑147 (MQ=255)
                                                                                             ttCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCgg                                                    >  6:133944/1‑147 (MQ=255)
                                                                                             ttCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTa                                                  >  8:109966/1‑149 (MQ=255)
                                                                                                   gCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCTCGGAGGATGCCGCGGTAGGTGAt                                            <  7:50630/149‑1 (MQ=255)
                                                                                                   gCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGGCGCGGTAGGTGAt                                            <  5:347624/149‑1 (MQ=255)
                                                                                                     cTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATgg                                          <  7:109966/149‑1 (MQ=255)
                                                                                                           gACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCCGGGAGGATGCCGCGGTAGGTGATGGGTAcc                                     >  6:228134/1‑148 (MQ=255)
                                                                                                            aCTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGaa                                                                          <  4:32003/110‑1 (MQ=255)
                                                                                                            aCTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGaa                                                                          >  3:32003/1‑110 (MQ=255)
                                                                                                                 gggCGGAATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGTGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACttg                               <  2:278861/148‑1 (MQ=255)
                                                                                                                     ggcATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGaa                          <  2:75070/149‑1 (MQ=255)
                                                                                                                     ggcATGTATTTATCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGaa                          <  1:156604/149‑1 (MQ=255)
                                                                                                                      gcATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACttgtt                             <  8:34448/145‑1 (MQ=255)
                                                                                                                      gcATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACttgtt                             >  7:34448/1‑145 (MQ=255)
                                                                                                                       cATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGt                        <  5:17313/149‑1 (MQ=255)
                                                                                                                         tGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCg                      <  7:68855/149‑1 (MQ=255)
                                                                                                                           tATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGt                    >  2:272275/1‑149 (MQ=255)
                                                                                                                            aTTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTc                       <  1:160207/145‑1 (MQ=255)
                                                                                                                            aTTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTc                       >  2:160207/1‑145 (MQ=255)
                                                                                                                            aTTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTa                  >  2:227275/1‑149 (MQ=255)
                                                                                                                               tctcGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTAc                 <  4:329175/148‑1 (MQ=255)
                                                                                                                                  tgaagtGCTGGCCAAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCAc              >  6:61407/5‑148 (MQ=255)
                                                                                                                                   gtcgtGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCa               <  5:133944/146‑1 (MQ=255)
                                                                                                                                   gtcgtGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACg             >  4:203221/1‑148 (MQ=255)
                                                                                                                                   gaagtGCTGGCCAAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggag                                                              <  2:45552/96‑1 (MQ=255)
                                                                                                                                   gaagtGCTGGCCAAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggag                                                              >  1:45552/4‑99 (MQ=255)
                                                                                                                                    tcgtGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGa            >  5:298200/1‑148 (MQ=255)
                                                                                                                                      gtGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAAc         <  4:266358/149‑1 (MQ=255)
                                                                                                                                            ggCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAACGCAc     <  6:94340/147‑1 (MQ=255)
                                                                                                                                            ggCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAACGCAc     <  3:227267/147‑1 (MQ=255)
                                                                                                                                             cccAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAACGCACGTc  <  5:228134/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAACGCACGTC  >  NC_002947/4914645‑4914934

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: