Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,914,784 T→G 5.1% K123Q (AAG→CAG)  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,7840TG 5.1% 254.3 / 4.0 79K123Q (AAG→CAG) ccmEcytochrome c‑type biogenesis protein CcmE
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base G (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 = 6.14e-01

AAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAAC  >  NC_002947/4914636‑4914927
                                                                                                                                                    |                                                                                                                                               
aaaGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCtt                                                                                                                                                 >  1:155735/1‑149 (MQ=255)
   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                                                                                                          <  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                                                                                            <  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                                                     <  6:302612/149‑1 (MQ=255)
                                                                                            gCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggagga                                                     <  1:259855/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                      >  7:34448/1‑145 (MQ=255)
                                                                                                                               gcATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACttgtt                      <  8:34448/145‑1 (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)
                                                                                                                                                    |                                                                                                                                               
AAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATGGGTACCGACTTGTTGAAGTCGGTTACCACGAAAC  >  NC_002947/4914636‑4914927

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: