Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,091,697 C→A 5.8% P41Q (CCG→CAG)  PP_5343 → putative Transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,091,6970CA 5.8% 223.1 / 3.3 69P41Q (CCG→CAG) PP_5343putative Transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base A (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 = 2.23e-01

TGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCCCCGGCTATTA  >  NC_002947/6091550‑6091844
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                    >  7:48678/1‑137 (MQ=255)
                                            cATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCgg                                                                                                                    <  8:48678/137‑1 (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                                                           >  5:7559/1‑143 (MQ=255)
                                                                                               cTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCAcgg                                                           <  6:7559/143‑1 (MQ=255)
                                                                                                  gCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGccc                                                                                             >  3:166465/1‑106 (MQ=255)
                                                                                                  gCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGccc                                                                                             <  4:166465/106‑1 (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                                                                   <  3:157497/109‑1 (MQ=255)
                                                                                                                         aCAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTg                                                                   >  4:157497/1‑109 (MQ=255)
                                                                                                                          cAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGa                          <  2:308015/149‑1 (MQ=255)
                                                                                                                               ccTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCt                     <  7:105388/149‑1 (MQ=255)
                                                                                                                                 ttCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCa                                                 <  3:67090/119‑1 (MQ=255)
                                                                                                                                 ttCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCa                                                 >  4:67090/1‑119 (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                       >  5:111772/1‑141 (MQ=255)
                                                                                                                                     ggCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGc                       <  6:111772/141‑1 (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                                                                                                        <  3:71367/54‑1 (MQ=255)
                                                                                                                                           gCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGc                                                                                                        >  4:71367/1‑54 (MQ=255)
                                                                                                                                                  ccGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCCCCGGCTAtt   <  4:128624/148‑1 (MQ=255)
                                                                                                                                                  ccGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCCCCGGCTATTa  <  1:195058/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGAACCATTCGACAAGGATCCCGCCATGTACAACAGTAAGCCGCATCAGGAACACGACCTCACCCGCTTGCACCAGCACATGCTAGACACCCGCCTGGCCGTTCTGGTCAGCCAGGGCGAACAGGGCCTTCTGGCCACGCACCTGCCGGTGCTGGTCGATACCGCCGAGGGCGATTTCGGCACTGTGTATGCCCACCTGGCCCGGGCCAACCGCCAATGGCAAGACCTGCAGCACGGCGGCGAAGCACTGCTGGTATTCCCCGGCGCCGATGCCTATGTCAGCCCCGGCTATTA  >  NC_002947/6091550‑6091844

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