Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,759,938 C→T 6.2% L335L (CTG→CTA)  PP_2413 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,759,9380CT 6.2% 189.8 / 2.7 64L335L (CTG→CTA) PP_2413GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGCACCAAGGGCTGCTGCGCTACCACGCCCCAG  >  NC_002947/2759791‑2760085
                                                                                                                                                   |                                                                                                                                                   
ccggccgATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCcaca                                                                                                                                                    <  5:58150/149‑1 (MQ=255)
 cggccgATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCATAg                                                                                                                                                   >  4:189499/1‑149 (MQ=255)
     cgATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCAGCGGCGAGTTTCCACAgcggg                                                                                                                                               >  7:343642/1‑149 (MQ=255)
      gATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAgcgggc                                                                                                                                              >  3:307240/1‑149 (MQ=255)
          cGCTGTTCCAGCAAGCTCCGGCTGAACACCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGcgggcggg                                                                                                                                           <  8:233504/148‑1 (MQ=255)
           gCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCgggcgggc                                                                                                                                          >  5:83357/1‑148 (MQ=255)
           gCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCgggcgggc                                                                                                                                          >  1:221836/1‑148 (MQ=255)
           gCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCgggcgggc                                                                                                                                          >  5:141568/1‑148 (MQ=255)
               ttGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGAt                                                                                                                                       <  8:234374/147‑1 (MQ=255)
                 gcagcaAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGt                                                                                                                                    <  7:350682/148‑1 (MQ=255)
                 gcagcaAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGt                                                                                                                                    <  4:322487/148‑1 (MQ=255)
                     caAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACggg                                                                                                                               <  3:228632/149‑1 (MQ=255)
                       aGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGgc                                                                                                                              >  5:74281/1‑148 (MQ=255)
                        gTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGgcca                                                                                                                            <  5:2546/149‑1 (MQ=255)
                         ttCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGgccag                                                                                                                           >  8:84271/1‑149 (MQ=255)
                            aGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGgccagcca                                                                                                                        >  2:49471/1‑149 (MQ=255)
                             ggCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGccac                                                                                                                       <  3:122624/149‑1 (MQ=255)
                                tGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGccacca                                                                                                                     >  2:360935/1‑148 (MQ=255)
                                  aaCAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCAcag                                                                                                                  <  6:74281/149‑1 (MQ=255)
                                   aCAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCAGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCAcagc                                                                                                                 <  8:343642/149‑1 (MQ=255)
                                    cAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCAcag                                                                                                                  >  4:310617/1‑147 (MQ=255)
                                           gcgcgcTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCt                                                                                                         >  8:157879/1‑149 (MQ=255)
                                                 ttGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAAc                                                                                                   >  7:122749/1‑149 (MQ=255)
                                                 ttGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAAc                                                                                                   >  2:305865/1‑149 (MQ=255)
                                                         ggCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAgcgc                                                                                             <  7:260161/147‑1 (MQ=255)
                                                         ggCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGcc                                                                                            >  8:206436/1‑148 (MQ=255)
                                                            cgccTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAgg                                                                                        >  3:336833/1‑149 (MQ=255)
                                                            cgccTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAgg                                                                                        <  2:117341/149‑1 (MQ=255)
                                                            cgccTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAgg                                                                                        <  3:101641/149‑1 (MQ=255)
                                                            cgccTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCCGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAgg                                                                                        <  1:305865/149‑1 (MQ=255)
                                                            cgccTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGAGAGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTCCCAACCAGCGCCAAgg                                                                                        >  6:129841/1‑149 (MQ=255)
                                                              ccTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGgc                                                                                      <  7:310079/149‑1 (MQ=255)
                                                                tCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCGCCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCGACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGgcgc                                                                                    >  2:234362/1‑149 (MQ=255)
                                                                  gAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGcc                                                                                   <  6:141568/148‑1 (MQ=255)
                                                                     gTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCATAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAAc                                                                               >  2:315543/1‑149 (MQ=255)
                                                                        ccAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGccacca                                                                                                                     >  7:299982/1‑108 (MQ=255)
                                                                        ccAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGccacca                                                                                                                     <  8:299982/108‑1 (MQ=255)
                                                                        ccAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAgg                                                                                        >  5:360924/1‑137 (MQ=255)
                                                                        ccAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAgg                                                                                        <  6:360924/137‑1 (MQ=255)
                                                                        ccAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTg                                                                             <  8:347264/148‑1 (MQ=255)
                                                                         cAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGc                                                                           >  3:333008/1‑149 (MQ=255)
                                                                          aGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGa                                                                                                                                        <  7:342176/87‑1 (MQ=255)
                                                                          aGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGa                                                                                                                                        >  8:342176/1‑87 (MQ=255)
                                                                           ggCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGa                                                                         <  2:293616/149‑1 (MQ=255)
                                                                           ggCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGa                                                                         <  4:297695/149‑1 (MQ=255)
                                                                           ggCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGa                                                                         >  7:35994/1‑149 (MQ=255)
                                                                                     tGATGCAATTGCTCGGCGGTACGGGCTGGGTCCTTCGTGCGGGCACCGGCGGCGAGTTTCCATAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAg                                                               <  3:189499/149‑1 (MQ=255)
                                                                                     tGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAg                                                               >  3:369116/1‑149 (MQ=255)
                                                                                           aaaTGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCATAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAg                                                         <  1:315543/149‑1 (MQ=255)
                                                                                            aaTGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAgg                                                        >  8:197483/1‑149 (MQ=255)
                                                                                                 tCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCg                                                                          <  3:31368/126‑1 (MQ=255)
                                                                                                 tCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCg                                                                          >  4:31368/1‑126 (MQ=255)
                                                                                                   ggcggTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGgtggtg                                                   <  4:336833/147‑1 (MQ=255)
                                                                                                        tACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGgc                                            <  5:200454/149‑1 (MQ=255)
                                                                                                          ccGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGgcg                                           <  8:21944/148‑1 (MQ=255)
                                                                                                              cgcgGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAggtgg                                                     <  1:184084/134‑1 (MQ=255)
                                                                                                              cgcgGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAggtgg                                                     >  2:184084/1‑134 (MQ=255)
                                                                                                                     cATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGc                                <  7:157879/148‑1 (MQ=255)
                                                                                                                     cATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGCa                               <  3:310617/149‑1 (MQ=255)
                                                                                                                       tCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACTAGTTTTTGCAACGGCGCCACGGTCTGCAcc                             <  8:35994/149‑1 (MQ=255)
                                                                                                                          tGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGCACCAAg                          <  6:38015/149‑1 (MQ=255)
                                                                                                                               gCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGgc                                                                                      >  4:266172/1‑84 (MQ=255)
                                                                                                                               gCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGgc                                                                                      <  3:266172/84‑1 (MQ=255)
                                                                                                                                          gAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGCACCAAGGGCTGCTGCGCTACCa          >  1:197939/1‑149 (MQ=255)
                                                                                                                                           aGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGCACCAAGGGCTGCTGCGCTACCAc         <  8:122749/149‑1 (MQ=255)
                                                                                                                                                  acaGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGCACCAAGGGCTGCTGCGCTACCACGCCCCAg  <  1:49471/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGGCCGATACGCTCTTGCAGCAAGTTCAGGCTGAACAGCAACGCGCGCTTGAGCTCGGCCGCCTCGAAGTACCAGGCCCGCACCTGATGCAAATGCTCGGCGGTACCGGCGCGGTCCATCGTGCGGGCACCGGCGGCGAGTTTCCACAGCGGGCGGGCGATGGTACGGGCCAGCCACCACAGCAGCAGGCTGCCAACCAGCGCCAAGGGCGCCGAACTGGCGATCACGTTCAGCACCAGGTGGTGCAACGGCGCCACCGTCTGCACCAAGGGCTGCTGCGCTACCACGCCCCAG  >  NC_002947/2759791‑2760085

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: ‑