Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,914,751 T→C 6.1% T134A (ACC→GCC)  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,7510TC 6.1% 188.4 / 2.2 65T134A (ACC→GCC) ccmEcytochrome c‑type biogenesis protein CcmE
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

GCAACGGCACACTGGCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATG  >  NC_002947/4914608‑4914893
                                                                                                                                               |                                                                                                                                              
gCAACGGCACACTGGCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACtt                                                                                                                                           >  8:97734/1‑149 (MQ=255)
     ggCACACTGGCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGggc                                                                                                                                      <  5:334445/149‑1 (MQ=255)
       cacaCTGGCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCTACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCggg                                                                                                                                       <  5:255184/146‑1 (MQ=255)
       cacaCTGGCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCTCCGGCAACGGACTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCagg                                                                                                                                       >  6:255184/1‑146 (MQ=255)
              gCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGt                                                                                                                                                >  4:177144/1‑130 (MQ=255)
              gCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGt                                                                                                                                                <  3:177144/130‑1 (MQ=255)
              gCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTAt                                                                                                                             >  7:100791/1‑149 (MQ=255)
              gCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTAt                                                                                                                             <  7:97734/149‑1 (MQ=255)
                   tACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTtct                                                                                                                         <  8:100791/148‑1 (MQ=255)
                         ggCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTg                                                                                                                  >  4:135966/1‑149 (MQ=255)
                            aaaGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCtt                                                                                                               >  4:94445/1‑149 (MQ=255)
                             aaGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTg                                                                                                              <  3:135966/149‑1 (MQ=255)
                                  gATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCtcgtc                                                                                                                     >  4:112971/1‑137 (MQ=255)
                                  gATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCtcgtc                                                                                                                     <  3:112971/137‑1 (MQ=255)
                                         aGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCa                                                                                                  >  6:79382/1‑149 (MQ=255)
                                            atcatcGCCATCTGGCCCAGTTCGGGGACTTCCATTGCTTCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCTGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCg                                                                                               <  2:5628/149‑1 (MQ=255)
                                            atcatcGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCg                                                                                               <  3:10064/149‑1 (MQ=255)
                                               atcGCCAACTGGACCAGTTCGTGGACTAACAGTGCTGCGTTCATGGCTTCCTCTCACCGCCGGCAGCGGCCTTGCAGCTTTCCTTAAGGGCCTTGATCACTTCTGGCGGCATGTATTTATCGTCGTGCTTGACCAGCACTTCATCGGc                                                                                             >  8:287792/1‑148 (MQ=255)
                                                  gCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTacc                                                                                         <  5:79382/149‑1 (MQ=255)
                                                     aaCTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTaccacc                                                                                      >  7:54348/1‑149 (MQ=255)
                                                            ccAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACtt                                                                                                    <  6:80705/128‑1 (MQ=255)
                                                            ccAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACtt                                                                                                    >  5:80705/1‑128 (MQ=255)
                                                            ccAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGt                                                                               >  1:34023/1‑149 (MQ=255)
                                                                        cTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCt                                                                   <  7:287792/149‑1 (MQ=255)
                                                                         tACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                        <  4:95013/143‑1 (MQ=255)
                                                                         tACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCAtt                                                                        >  3:95013/1‑143 (MQ=255)
                                                                         tACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCtt                                                                  <  8:54348/149‑1 (MQ=255)
                                                                           ccAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTg                                                                 >  1:54124/1‑148 (MQ=255)
                                                                           ccAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGc                                                                >  3:39768/1‑149 (MQ=255)
                                                                                gctgcGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCGCACCGTCGGCATTGAGCTTGCCAAgg                                                           <  3:94445/149‑1 (MQ=255)
                                                                                 ctgcGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATc                                                                                                <  4:119233/111‑1 (MQ=255)
                                                                                 ctgcGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATc                                                                                                >  3:119233/1‑111 (MQ=255)
                                                                                     gTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGcgac                                                      <  3:88591/149‑1 (MQ=255)
                                                                                      ttCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGcgac                                                      >  5:267825/1‑148 (MQ=255)
                                                                                       tCATGGCTTCGTCTCCCCGCCCGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGcgacga                                                    >  4:254926/1‑149 (MQ=255)
                                                                                        cATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAggg                                                          >  7:176538/1‑142 (MQ=255)
                                                                                        cATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAggg                                                          <  8:176538/142‑1 (MQ=255)
                                                                                           ggCTTCGTCCCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgcc                                                >  1:303979/1‑149 (MQ=255)
                                                                                                cGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccctgc                                            <  6:267825/148‑1 (MQ=255)
                                                                                                 gTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccctgccc                                          <  2:112281/149‑1 (MQ=255)
                                                                                                     ccccGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTcg                                      >  8:174826/1‑149 (MQ=255)
                                                                                                       ccgccgGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGCGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccctgcc                                           <  5:199910/142‑1 (MQ=255)
                                                                                                       ccgccgGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGCGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccctgcc                                           >  6:199910/1‑142 (MQ=255)
                                                                                                       ccgccgGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGCGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccctgcc                                           <  1:322123/142‑1 (MQ=255)
                                                                                                       ccgccgGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGCGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGAtgccctgcc                                           >  2:322123/1‑142 (MQ=255)
                                                                                                        cgccgGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCgg                                   <  2:303979/149‑1 (MQ=255)
                                                                                                         gccgGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCgg                                   >  3:339407/1‑148 (MQ=255)
                                                                                                           ggaaaGCGGCCAGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTTGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTcgcg                                    <  7:11772/141‑1 (MQ=255)
                                                                                                           ggaaaGCGGCCAGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTTGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTcgcg                                    >  8:11772/5‑145 (MQ=255)
                                                                                                            ggCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGaa                                 >  5:335599/1‑147 (MQ=255)
                                                                                                                     cTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGgga                      <  6:335599/149‑1 (MQ=255)
                                                                                                                        gCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggagg                    <  7:174826/148‑1 (MQ=255)
                                                                                                                        gCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggag                     <  4:39768/147‑1 (MQ=255)
                                                                                                                          cGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATg                 >  1:265287/1‑149 (MQ=255)
                                                                                                                               tCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCcgcg            >  4:76772/1‑149 (MQ=255)
                                                                                                                                 cTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCgg           <  2:265287/148‑1 (MQ=255)
                                                                                                                                 cTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCgg           <  2:54124/148‑1 (MQ=255)
                                                                                                                                 cTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGt          >  8:99387/1‑149 (MQ=255)
                                                                                                                                     aGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggagg                    <  8:194392/135‑1 (MQ=255)
                                                                                                                                     aGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGggagg                    >  7:194392/1‑135 (MQ=255)
                                                                                                                                       ggacTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCATCACTTCATCGGCTACCACCACACCCTCGGCAATGAGCTTACCAAGGGCTACGATGCCCTGCCATTAGCGTAACAGGTCCGGGAGGATGCCGAGGTAGGTgc    >  8:290529/4‑148 (MQ=255)
                                                                                                                                         ccTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATg  <  4:339407/149‑1 (MQ=255)
                                                                                                                                          cTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACAcc                                                                                 <  3:56616/69‑1 (MQ=255)
                                                                                                                                          cTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACAcc                                                                                 >  4:56616/1‑69 (MQ=255)
                                                                                                                                           ttGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAg        <  5:210464/141‑1 (MQ=255)
                                                                                                                                           ttGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAg        >  6:210464/1‑141 (MQ=255)
                                                                                                                                             ggTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCaa                                                             >  8:70304/1‑86 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GCAACGGCACACTGGCCTGTACGACGGCAAAGCAGATGGCCAGAATCATCGCCAACTGGCCCAGTTCGGGGACTACCAGTGCTGCGTTCATGGCTTCGTCTCCCCGCCGGCAGCGGCCTGGCCGCTTTCCTTCAGGGCCTTGGTGACTTCGGGCGGCATGTATTTCTCGTCGTGCTTGGCCAGCACTTCATCGGCTACCACCACACCGTCGGCATTGAGCTTGCCAAGGGCGACGATGCCCTGCCCTTCGCGGAACAGGTCGGGGAGGATGCCGCGGTAGGTGATG  >  NC_002947/4914608‑4914893

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: