Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 866,735 A→C 5.8% E63A (GAA→GCA) ‡ PP_0748 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947866,7350AC 5.8% 314.3 / 7.1 99E63A (GAA→GCA) ‡PP_0748hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/48);  new base C (3/3);  total (47/52)
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.05e-01

GGCCGAAAAGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTTGGTGGTGCTGG  >  NC_002947/866592‑866881
                                                                                                                                               |                                                                                                                                                  
ggCCGAAAAGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCt                                                                                                                                               <  7:324675/149‑1 (MQ=255)
  ccGAAAAGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCtt                                                                                                                                              >  1:18284/1‑148 (MQ=255)
   cGAAAAGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGa                                                                                                                                            <  6:348625/149‑1 (MQ=255)
    gAAAAGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGa                                                                                                                                            >  5:107461/1‑148 (MQ=255)
    gAAAAGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGAt                                                                                                                                           <  5:183201/149‑1 (MQ=255)
     aaaaGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGa                                                                                                                                            >  5:81900/1‑147 (MQ=255)
         ggCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAgcgc                                                                                                                                      <  2:130701/149‑1 (MQ=255)
            cGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTc                                                                                                                                    >  3:91998/1‑148 (MQ=255)
             gCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCaa                                                                                                                                  <  5:289752/149‑1 (MQ=255)
                gCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAggg                                                                                                                               >  1:404251/1‑149 (MQ=255)
                gCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAggg                                                                                                                               <  8:158695/149‑1 (MQ=255)
                 ccAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTg                                                                                                                                             <  3:209815/134‑1 (MQ=255)
                 ccAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTg                                                                                                                                             >  4:209815/1‑134 (MQ=255)
                  cAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGc                                                                                                                             >  6:86045/1‑149 (MQ=255)
                  cAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGc                                                                                                                             >  2:77659/1‑149 (MQ=255)
                   aaCCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCt                                                                                                                                               <  8:205368/130‑1 (MQ=255)
                   aaCCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCg                                                                                                                            >  4:173217/1‑149 (MQ=255)
                   aaCCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCCTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCCCTTGAGCGCTATCAAGGTTTCAGCCTTGTCGGCGCTCaaatggc                                                                                                                             >  4:384215/1‑143 (MQ=21)
                    aCCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTCAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGt                                                                                                                           <  6:81900/149‑1 (MQ=255)
                     cccGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGtt                                                                                                                          >  4:392767/1‑149 (MQ=255)
                      ccGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGttt                                                                                                                         >  5:309946/1‑149 (MQ=255)
                       cgcgcgGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCCCTTGAGCGCTATCAAGGTTTCAGCCTTGTCGGCGCTCAAATGGCgatgc                                                                                                                        <  3:384215/149‑5 (MQ=21)
                       cgcgcgGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGtttt                                                                                                                        <  5:86045/149‑1 (MQ=255)
                          gcgggcggCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGc                                                                                                                     <  6:107461/149‑1 (MQ=255)
                            gggcggCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGGTAGCGCTCAAGGGGCGTTTTGGCCt                                                                                                                   >  5:41810/1‑149 (MQ=255)
                            gggcggCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCt                                                                                                                   <  4:188714/149‑1 (MQ=255)
                                ggcAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGtt                                                                                                                >  5:49371/1‑148 (MQ=255)
                                ggcAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGtt                                                                                                                >  8:31353/1‑148 (MQ=255)
                                 gcagcaCCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGa                                                                                                              <  8:298845/149‑1 (MQ=255)
                                 gcagcaCCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCATGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGa                                                                                                              >  1:179466/1‑149 (MQ=255)
                                  cagcaCCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGAt                                                                                                             <  4:91998/149‑1 (MQ=255)
                                   agcaCCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATc                                                                                                            <  3:392767/149‑1 (MQ=255)
                                   agcaCCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATc                                                                                                            >  6:93957/1‑149 (MQ=255)
                                      aCCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTc                                                                                                         >  2:39692/1‑149 (MQ=255)
                                       ccATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGGTAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCg                                                                                                        <  6:41810/149‑1 (MQ=255)
                                        cATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGTAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGa                                                                                                       <  8:370384/149‑1 (MQ=255)
                                          tAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGACGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATg                                                                                                     <  5:351053/149‑1 (MQ=255)
                                              tttCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCTATCAAGGTTTCAGCCTTGTCGGCGCTCAAATGGCg                                                                                                                            >  8:112238/1‑122 (MQ=21)
                                              tttCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAAc                                                                                                  >  7:320493/1‑148 (MQ=255)
                                              tttCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGGCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCTATCAAGGTTTCAGCCTTGTCGGCGCTCAAATGGCg                                                                                                                            <  7:112238/122‑1 (MQ=21)
                                                tCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGTGGCGTTTTGGCCTGTTGATCTTCGATGAACAAg                                                                                               >  4:366865/1‑149 (MQ=255)
                                                   cGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCTATCAAGGTTTCAGCCTTGTCGGCGCTCaaa                                                                                                                                 <  8:359492/112‑2 (MQ=25)
                                                   cGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCTATCAAGGTTTCAGCCTTGTCGGCGCTCaaa                                                                                                                                 >  7:359492/1‑111 (MQ=25)
                                                   cGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGgcg                                                                                            <  7:31353/149‑1 (MQ=255)
                                                        tGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCgg                                                                                         <  8:320493/147‑1 (MQ=255)
                                                          gAAGAGCCCTGCCACCGCAGGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGtgac                                                                                     <  4:78741/149‑1 (MQ=255)
                                                           aaGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGtgact                                                                                    >  4:387485/1‑149 (MQ=255)
                                                              agCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGtgactgac                                                                                 <  6:68308/149‑1 (MQ=255)
                                                                  cTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGc                                                                             >  6:14479/1‑149 (MQ=255)
                                                                     ccACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCggg                                                                          <  2:179466/149‑1 (MQ=255)
                                                                            aTGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGtgctg                                                                     >  7:116263/1‑147 (MQ=255)
                                                                             tGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAg                                                                  <  1:131854/149‑1 (MQ=255)
                                                                               gCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGc                                                                <  4:351117/149‑1 (MQ=255)
                                                                                   tCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGCCAAGCGGGTGCTGCAGGCCGGt                                                            >  5:248333/1‑149 (MQ=255)
                                                                                         tGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGAc                                                      >  3:100530/1‑149 (MQ=255)
                                                                                         tGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGGGACTGACAAGCGGGTGCTGCAGGACGGTGGCGAc                                                      >  5:415148/1‑149 (MQ=255)
                                                                                         tGCAACCAAGCACATACATTCCGCCGGATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGAc                                                      <  8:297465/149‑1 (MQ=255)
                                                                                            aaCCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGt                                                                                        >  1:329074/1‑112 (MQ=255)
                                                                                            aaCCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGt                                                                                        <  2:329074/112‑1 (MQ=255)
                                                                                             aCCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTc                                                  >  2:209744/1‑149 (MQ=255)
                                                                                                aaGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTgg                                               >  4:82382/1‑149 (MQ=255)
                                                                                                 aGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTggg                                              <  8:300230/149‑1 (MQ=255)
                                                                                                  gCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGgt                                             >  6:227403/1‑149 (MQ=255)
                                                                                                    aCCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGgtgt                                           <  6:309946/149‑1 (MQ=255)
                                                                                                    aCCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGgtgt                                           >  7:20345/1‑149 (MQ=255)
                                                                                                     ccTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTgg                                                          <  2:290396/133‑1 (MQ=255)
                                                                                                     ccTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGtgtg                                          <  1:77659/149‑1 (MQ=255)
                                                                                                      cTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGtgtg                                          >  8:261867/1‑148 (MQ=255)
                                                                                                       tACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGa                                        >  2:351635/1‑149 (MQ=255)
                                                                                                        aCATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGAt                                       >  1:295932/1‑149 (MQ=255)
                                                                                                         cATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATc                                      >  3:163167/1‑149 (MQ=255)
                                                                                                         cATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCACGCCGGTGGCGACTGTCTGGGTGTGGATc                                      <  1:231529/149‑1 (MQ=255)
                                                                                                          aTTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCt                                     <  7:261867/149‑1 (MQ=255)
                                                                                                          aTTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCt                                     <  4:203979/149‑1 (MQ=255)
                                                                                                              cgccgcATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGgc                                  <  5:14479/148‑1 (MQ=255)
                                                                                                                        ccATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGgt                                             <  1:239826/127‑1 (MQ=255)
                                                                                                                        ccATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGgt                                             >  2:239826/1‑127 (MQ=255)
                                                                                                                        ccATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGtacc                        >  8:291703/1‑148 (MQ=255)
                                                                                                                        ccATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGtacc                        >  7:322733/1‑148 (MQ=255)
                                                                                                                          aTTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGtacctac                     <  8:20345/149‑1 (MQ=255)
                                                                                                                              gAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCAt                  <  3:173217/148‑1 (MQ=255)
                                                                                                                               aGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCAta                 <  1:209744/148‑2 (MQ=255)
                                                                                                                               aGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATgg                >  3:251323/1‑149 (MQ=255)
                                                                                                                                gcgcCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGc               <  6:49371/149‑1 (MQ=255)
                                                                                                                                gcgcCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGc               >  2:243257/1‑149 (MQ=255)
                                                                                                                                gcgcCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGc               <  3:387485/149‑1 (MQ=255)
                                                                                                                                    cGACAAGGCTGAAACCTTGTTTGCGCTCAAGGGGCGTTTTGGCCTGTTGAGCTTCGATGAACAATGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTtgg           <  8:376001/149‑1 (MQ=255)
                                                                                                                                     gACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTg                             >  4:386799/1‑130 (MQ=255)
                                                                                                                                      aCAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTtggtg         <  2:404251/149‑1 (MQ=255)
                                                                                                                                      aCAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTtggtg         <  3:82382/149‑1 (MQ=255)
                                                                                                                                       catggCTGAAACCTTTATAGCGCTCATGGGCCGTTTTGGCCTGTTGATCTTCTATGAACAAGGCGCGGTGACTGACTCGGGGGTGCTGCAGGCCGGTGGGGACTGTGTGGGTGTGGTTCTGGCGCGTGGAACCTACCATGGCTtggggg        <  8:116263/146‑1 (MQ=255)
                                                                                                                                           gCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTtggtg         <  4:389830/144‑1 (MQ=255)
                                                                                                                                           gCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTtggtg         >  3:389830/1‑144 (MQ=255)
                                                                                                                                            cTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTTGGTGGTGCTg   <  2:18284/149‑1 (MQ=255)
                                                                                                                                             tGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCt                                     >  7:29063/1‑114 (MQ=255)
                                                                                                                                             tGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCt                                     <  8:29063/114‑1 (MQ=255)
                                                                                                                                              gAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTgg                                                          <  4:31828/92‑1 (MQ=255)
                                                                                                                                              gAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTgg                                                          >  3:31828/1‑92 (MQ=255)
                                                                                                                                              gAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTTGGTGGTGCTg   >  3:298618/1‑147 (MQ=255)
                                                                                                                                               aaaCCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTTGGTGGTGCTgg  <  5:227403/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GGCCGAAAAGGCCGCTGCCAACCCGCGCGGGCGGCAGCACCATAACTTTCACGCTATGGAAGAGCCCTGCCACCGCATGGCCGTCGGGCTGCAACCAAGCACCTACATTCCGCCGCATCGCCATTTGAGCGCCGACAAGGCTGAAACCTTGATAGCGCTCAAGGGGCGTTTTGGCCTGTTGATCTTCGATGAACAAGGCGCGGTGACTGACAAGCGGGTGCTGCAGGCCGGTGGCGACTGTCTGGGTGTGGATCTGGCGCCTGGTACCTACCATGGCTTGGTGGTGCTGG  >  NC_002947/866592‑866881

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: