Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,091,694 T→G 5.8% L40R (CTG→CGG)  PP_5343 → putative Transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,091,6940TG 5.8% 220.7 / 4.2 69L40R (CTG→CGG) PP_5343putative Transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base G (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

TGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCC  >  NC_002947/6091550‑6091834
                                                                                                                                                |                                                                                                                                            
tGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCg                                                                                                                                          <  6:198212/149‑1 (MQ=255)
   aaCCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTg                                                                                                                                       >  4:5879/1‑149 (MQ=255)
   aaCCATTCGACAAGGATCCAGCCATGTACAACAGTAAGCCGCATCAGGTACACGACCTCACACGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCTGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTg                                                                                                                                       <  6:176018/149‑1 (MQ=255)
        ttCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCGGCAGGTGcgtgg                                                                                                                                  <  2:223072/149‑5 (MQ=255)
        ttCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCA‑GCACCGGCAGGTGcgtggc                                                                                                                                 <  7:134160/149‑6 (MQ=255)
              aaGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATAc                                                                                                                            >  3:68133/1‑149 (MQ=255)
                ggATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATAccg                                                                                                                          >  3:264187/1‑149 (MQ=255)
                 gATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATAccgc                                                                                                                         <  1:267897/149‑1 (MQ=255)
                      cGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAgg                                                                                                                    <  4:68133/149‑1 (MQ=255)
                      cGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAgg                                                                                                                    >  3:127082/1‑149 (MQ=255)
                         cATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGc                                                                                                                  >  6:47242/1‑148 (MQ=255)
                         cATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCg                                                                                                                 <  3:5879/149‑1 (MQ=255)
                                aacaGTAAGCCGCATCAGTAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCgg                                                                                                          <  6:129419/149‑1 (MQ=255)
                                aacaGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCTGTGCTGGTCGATACCGCCGATGGCGATTTCgg                                                                                                          >  6:215797/1‑149 (MQ=255)
                                  caGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCa                                                                                                        >  8:105388/1‑149 (MQ=255)
                                    gTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATTCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACt                                                                                                      <  2:125768/149‑1 (MQ=255)
                                        gCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACtgtg                                                                                                   >  3:118865/1‑148 (MQ=255)
                                           gCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATg                                                                                               <  4:264187/149‑1 (MQ=255)
                                            cATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCgg                                                                                                          <  8:48678/137‑1 (MQ=255)
                                            cATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCgg                                                                                                          >  7:48678/1‑137 (MQ=255)
                                               cAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCa                                                                                           >  2:195058/1‑149 (MQ=255)
                                                          ccTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCggg                                                                                <  3:242932/149‑1 (MQ=255)
                                                           cTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGc                                                                               >  7:1768/1‑149 (MQ=255)
                                                               cccGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATAccgcc                                                                                                                        >  3:186980/1‑104 (MQ=255)
                                                               cccGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATAccgcc                                                                                                                        <  4:186980/104‑1 (MQ=255)
                                                                  gCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGc                                                                        <  4:118865/149‑1 (MQ=255)
                                                                    ttGCACCAGCACATGCTAGACACCCGCCTGGCTGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAAcc                                                                          >  5:99786/1‑145 (MQ=255)
                                                                    ttGCACCAGCACATGCTAGACACCCGCCTGGCTGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAAcc                                                                          <  6:99786/145‑1 (MQ=255)
                                                                        aCCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATgg                                                                  >  5:259838/1‑149 (MQ=255)
                                                                          cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGcc                                                                              <  8:103805/135‑1 (MQ=255)
                                                                          cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGcc                                                                              >  7:103805/1‑135 (MQ=255)
                                                                          cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGc                                                                 >  5:280437/1‑148 (MQ=255)
                                                                          cAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCa                                                                <  5:60510/149‑1 (MQ=255)
                                                                                   cTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTgca                                                       <  8:1768/149‑1 (MQ=255)
                                                                                      gACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTgcagca                                                    >  7:232901/1‑149 (MQ=255)
                                                                                       acacCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGTGCGATTTCTGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAc                                                   <  5:215797/149‑1 (MQ=255)
                                                                                           ccGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcggcg                                               >  8:243460/1‑149 (MQ=255)
                                                                                           ccGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcggcg                                               >  1:308015/1‑149 (MQ=255)
                                                                                               cTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcgg                                                 <  6:7559/143‑1 (MQ=255)
                                                                                               cTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcgg                                                 >  5:7559/1‑143 (MQ=255)
                                                                                                  gCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGccc                                                                                   <  4:166465/106‑1 (MQ=255)
                                                                                                  gCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGccc                                                                                   >  3:166465/1‑106 (MQ=255)
                                                                                                    cGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATTCCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGccc                                                                                   >  1:252406/1‑104 (MQ=255)
                                                                                                    cGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGccc                                                                                   <  2:252406/104‑1 (MQ=255)
                                                                                                    cGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCAc                                      >  3:39820/1‑149 (MQ=255)
                                                                                                      ttCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCActg                                    >  2:89415/1‑149 (MQ=255)
                                                                                                           gTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCActgct                                  <  6:280437/146‑1 (MQ=255)
                                                                                                           gTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTgg                                >  3:128624/1‑148 (MQ=255)
                                                                                                               gCCAGGGCGAACCGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTc                           <  3:77919/149‑1 (MQ=255)
                                                                                                               gCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTc                           <  4:39820/149‑1 (MQ=255)
                                                                                                                   gggCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCg                       <  6:259838/149‑1 (MQ=255)
                                                                                                                         aCAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTg                                                         >  4:157497/1‑109 (MQ=255)
                                                                                                                         aCAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTg                                                         <  3:157497/109‑1 (MQ=255)
                                                                                                                          cAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGa                <  2:308015/149‑1 (MQ=255)
                                                                                                                               ccTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCt           <  7:105388/149‑1 (MQ=255)
                                                                                                                                 ttCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCa                                       >  4:67090/1‑119 (MQ=255)
                                                                                                                                 ttCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCa                                       <  3:67090/119‑1 (MQ=255)
                                                                                                                                   cTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGTCCAACCGCCAATGGCAAGACCTGCAGCACGGAGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTAt         <  8:232901/147‑1 (MQ=255)
                                                                                                                                   cTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATg        <  1:89415/148‑1 (MQ=255)
                                                                                                                                    tGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGt       >  7:136492/1‑148 (MQ=255)
                                                                                                                                    tGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGt       >  1:1992/1‑148 (MQ=255)
                                                                                                                                     ggCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGc             <  6:111772/141‑1 (MQ=255)
                                                                                                                                     ggCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGc             >  5:111772/1‑141 (MQ=255)
                                                                                                                                      gaccagCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCat    >  7:65450/6‑148 (MQ=255)
                                                                                                                                       ccACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGc   <  1:229092/149‑1 (MQ=255)
                                                                                                                                        ccagCACCGGCAGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGcc  >  8:217802/4‑149 (MQ=255)
                                                                                                                                        cACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGcc  >  2:55971/1‑149 (MQ=255)
                                                                                                                                           gCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGc                                                                                              >  4:71367/1‑54 (MQ=255)
                                                                                                                                           gCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGc                                                                                              <  3:71367/54‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCC  >  NC_002947/6091550‑6091834

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