Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,620,295 T→G 8.1% S20R (AGC→CGC)  PP_3190 ← ornithine cyclodeaminase/mu‑crystallin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,620,2950TG 8.1% 228.1 / 10.1 74S20R (AGC→CGC) PP_3190ornithine cyclodeaminase/mu‑crystallin family protein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base G (4/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGCGTTCGGCGAGCCT  >  NC_002947/3620148‑3620441
                                                                                                                                                     |                                                                                                                                                   
aaaGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCt                                                                                                                                                     >  4:286659/1‑148 (MQ=255)
       ttgcgtgaagtcggcctacatcgaaggcagcccgaccttcacGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTg                                                                                                                                                    <  6:329555‑M2/100‑1 (MQ=255)
       ttgcgtgaagtcggcctacatcgaaggcagcccgaccttcacGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTg                                                                                                                                                    >  5:329555‑M2/43‑142 (MQ=255)
        tgcgtgaagtcggcctacatcgaaggcagcccgaccttcacGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCaa                                                                                                                                            >  3:122169‑M2/42‑149 (MQ=255)
          ggTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACt                                                                                                                                          >  2:173874/1‑149 (MQ=255)
           GTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCGGTCAGCCAActg                                                                                                                                         >  8:23383‑M1/1‑148 (MQ=255)
           GTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGGCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCGGTCAGCCAActg                                                                                                                                         >  6:92668‑M1/1‑148 (MQ=255)
            TGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCATTTGGCTGACGGCCAACtcg                                                                                                                                        >  2:136764‑M1/1‑148 (MQ=255)
             GAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCg                                                                                                                                       <  4:319768‑M1/149‑1 (MQ=255)
                        tGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATgg                                                                                                                            >  6:51413/1‑149 (MQ=255)
                         gCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGCCACCAGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTAGGCTGACCGCCAACTCAGCTTTGAAATggc                                                                                                                           >  2:231518/1‑148 (MQ=255)
                          ccTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTTCACCGGCGGCAACTTCACCTTTCCTTCGGAATAGGCGGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATggg                                                                                                                           <  1:85513/149‑1 (MQ=255)
                          ccTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGaaa                                                                                                                               >  1:159452/1‑144 (MQ=255)
                          ccTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGaaa                                                                                                                               <  2:159452/144‑1 (MQ=255)
                            ttCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATggg                                                                                                                           >  7:298019/1‑146 (MQ=255)
                            ttCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATggg                                                                                                                           <  8:298019/146‑1 (MQ=255)
                            ttCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGcga                                                                                                                        >  2:79137/1‑149 (MQ=255)
                            ttCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGc                                                                                                                          <  7:57641/147‑1 (MQ=255)
                            ttCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGc                                                                                                                          >  8:57641/1‑147 (MQ=255)
                            ttCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCAGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGcga                                                                                                                        >  5:204459/1‑149 (MQ=255)
                                 tGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTg                                                                                                                                  >  5:246990/1‑134 (MQ=255)
                                 tGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTg                                                                                                                                  <  6:246990/134‑1 (MQ=255)
                                  gTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCt                                                                                                                  >  6:113909/1‑149 (MQ=255)
                                   tAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTg                                                                                                                <  1:173874/149‑1 (MQ=255)
                                    aGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCGGTCAGCCAACtg                                                                                                                                         <  1:196666/124‑2 (MQ=255)
                                    aGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCGGTCAGCCAACtg                                                                                                                                         >  2:196666/1‑123 (MQ=255)
                                     ggCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTgcg                                                                                                              >  6:92935/1‑149 (MQ=255)
                                      gCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTgcgc                                                                                                             >  2:294455/1‑149 (MQ=255)
                                      gCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTgcgc                                                                                                             <  3:319768/149‑1 (MQ=255)
                                       ccGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGcc                                                                                                            >  4:254180/1‑149 (MQ=255)
                                         gACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCtt                                                                                                          <  2:128/149‑1 (MQ=255)
                                          aCTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTg                                                                                                         >  3:167741/1‑149 (MQ=255)
                                               acGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGaa                                                                                                    >  6:317213‑M2/3‑149 (MQ=255)
                                               acGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGaa                                                                                                    <  5:51413‑M2/147‑1 (MQ=255)
                                                          aCCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCAtt                                                                                          <  4:167741/148‑1 (MQ=255)
                                                          aCCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCAtt                                                                                          <  5:113909/148‑1 (MQ=255)
                                                           ccGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTg                                                                                         >  3:1754/1‑148 (MQ=255)
                                                            cGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGaa                                                                                       >  6:83542/1‑149 (MQ=255)
                                                             gTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAAc                                                                                      <  3:286659/149‑1 (MQ=255)
                                                              ttCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCACCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACg                                                                                     >  3:188792/1‑149 (MQ=255)
                                                                   gcCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCCTTCCTCCAGTTGGCTGACCGCCAACTCCggg                                                                                                                                     <  1:218660/98‑1 (MQ=255)
                                                                   gcCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCggg                                                                                                                                     >  2:218660/1‑97 (MQ=255)
                                                                     ccTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGa                                                                              <  3:254180/149‑1 (MQ=255)
                                                                     ccTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCACCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGa                                                                              <  4:188792/149‑1 (MQ=255)
                                                                      cTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAg                                                                             >  5:126609/1‑149 (MQ=255)
                                                                         gAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAgg                                                                            <  2:48751/147‑1 (MQ=255)
                                                                           aTGCGAAGCCCTTCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCGAACTCCGGGTTGAAATGGGCGAGGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAAc                                                                         <  1:231518/148‑1 (MQ=255)
                                                                               gAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTc                                                                    <  5:289303/149‑1 (MQ=255)
                                                                                aaGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTcc                                                                   <  4:122169/149‑1 (MQ=255)
                                                                                 aGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCca                                                                  <  5:317213/149‑1 (MQ=255)
                                                                                  gCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTAGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCca                                                                  <  1:79137/149‑1 (MQ=255)
                                                                                         aCGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGaa                                                                          <  6:113646/133‑1 (MQ=255)
                                                                                         aCGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGaa                                                                          >  5:113646/1‑133 (MQ=255)
                                                                                         aCGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCcagca                                                               >  8:143708/1‑144 (MQ=255)
                                                                                         aCGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCcagca                                                               <  7:143708/144‑1 (MQ=255)
                                                                                            ggcggcACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTa                                                       <  8:38871/149‑1 (MQ=255)
                                                                                                gcacctgcacctTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGa                                                                                                     >  5:52187/1‑99 (MQ=255)
                                                                                                gcacctgcacctTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGa                                                                                                     <  6:52187/99‑1 (MQ=255)
                                                                                                 cacctgcacctTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAg                                                                             >  5:94803/1‑122 (MQ=255)
                                                                                                 cacctgcacctTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAg                                                                             <  6:94803/122‑1 (MQ=255)
                                                                                                       cacctTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAg                                            <  6:204459/149‑1 (MQ=255)
                                                                                                       cacctTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAg                                            <  4:1754/149‑1 (MQ=255)
                                                                                                        acctTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAg                                            >  4:293008/1‑148 (MQ=255)
                                                                                                                          gCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCca                                                                  <  6:136716/108‑1 (MQ=255)
                                                                                                                          gCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCca                                                                  >  5:136716/1‑108 (MQ=255)
                                                                                                                               tGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGc                                   <  2:54779/135‑1 (MQ=255)
                                                                                                                               tGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGc                                   >  1:54779/1‑135 (MQ=255)
                                                                                                                                  aGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGAt                  <  3:293008/149‑1 (MQ=255)
                                                                                                                                      cTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGc               >  3:24122/1‑149 (MQ=255)
                                                                                                                                      cTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGc               <  4:24122/149‑1 (MQ=255)
                                                                                                                                         cccggaGTTGGCGGTCAGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGCGtt            >  1:174084/6‑149 (MQ=255)
                                                                                                                                           tccAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCCCACAGGATCACGGCGACACGCTTGCCCTGATAGCGTTCg          >  8:51000/1‑149 (MQ=255)
                                                                                                                                               gTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACACTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGCGTTCGGc        >  6:197046/1‑147 (MQ=255)
                                                                                                                                               gTTGGCGGTCAGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGCGTTCGGCGa      <  2:174084/149‑1 (MQ=255)
                                                                                                                                                   gCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGCGTTCGGCGAGCct  >  3:74403/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGCACGGGCGGCACCTGCACCTTGCCTTCGGAATAGGCGATGAAGCCTTCCTCCAGTTGGCTGACCGCCAACTCCGGGTTGAAATGGGCGACGATCTGCGCCTTGCTGAAGACCCTCATTGAACGGCCTCGAGGAACTTTTCCAGCAGGATGTTACGGCCACACAGGATCACGGCGACACGCTTGCCCTGATAGCGTTCGGCGAGCCT  >  NC_002947/3620148‑3620441

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: