Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,914,786 G→C 5.1% A122G (GCC→GGC)  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,7860GC 5.1% 259.3 / 3.6 79A122G (GCC→GGC) ccmEcytochrome c‑type biogenesis protein CcmE
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base C (2/2);  total (41/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 = 3.47e-01

GCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAACGCAC  >  NC_002947/4914639‑4914931
                                                                                                                                                   |                                                                                                                                                 
gCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACGAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGc                                                                                                                                                  <  5:88066/149‑1 (MQ=255)
    aTGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGc                                                                                                                                                  <  7:322704/145‑1 (MQ=255)
    aTGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGc                                                                                                                                                  >  8:322704/1‑145 (MQ=255)
      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                                                                                                              <  7:8037/149‑1 (MQ=255)
                                    ggggACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                                                              <  6:323760/149‑1 (MQ=255)
                                     gggACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                                                              >  7:5112/1‑148 (MQ=255)
                                     gggACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                                                              >  1:257991/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                                                                                                <  5:128561/119‑1 (MQ=255)
                                                                            ggaaaGCGGCCAGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAggg                                                                                                >  6:128561/5‑123 (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                                                 >  6:133944/1‑147 (MQ=255)
                                                                                                   ttCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCgg                                                 >  2:122163/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                    >  2:160207/1‑145 (MQ=255)
                                                                                                                                  aTTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTc                    <  1:160207/145‑1 (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)
                                                                                                                                                   |                                                                                                                                                 
GCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAACGCAC  >  NC_002947/4914639‑4914931

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: