Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,755,634 G→A 7.0% G359R (GGG→AGG)  PP_3319 → GGDEF domain‑containing sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,755,6340GA 7.0% 251.5 / 7.1 86G359R (GGG→AGG) PP_3319GGDEF domain‑containing sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (43/37);  new base A (3/3);  total (46/40)
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.18e-01

TGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGAAGACGCTGTCC  >  NC_002947/3755492‑3755778
                                                                                                                                              |                                                                                                                                                 
tGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCact                                                                                                                                             <  3:320589/149‑1 (MQ=255)
  aCGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCactac                                                                                                                                           >  5:136565/1‑149 (MQ=255)
      gCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCg                                                                                                                                      >  5:364260/1‑149 (MQ=255)
            gggTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGAt                                                                                                                                 >  6:118066/1‑148 (MQ=255)
            gggTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATc                                                                                                                                <  7:293812/149‑1 (MQ=255)
              gTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCaa                                                                                                                              <  3:176443/149‑1 (MQ=255)
                    ttAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCt                                                                                                                         >  3:98119/1‑148 (MQ=255)
                    ttAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCt                                                                                                                         >  3:263240/1‑148 (MQ=255)
                    ttAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCt                                                                                                                         >  3:258308/1‑148 (MQ=255)
                    ttAACCGCCTGTCTCTGGAAGACCGGCTTCTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCt                                                                                                                         >  8:125768/1‑148 (MQ=255)
                        ccGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATc                                                                                                                     >  8:210690/1‑148 (MQ=255)
                               tctcTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTgg                                                                                                              >  5:263698/1‑148 (MQ=255)
                                    ggAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgc                                                                                                        >  2:294613/1‑149 (MQ=255)
                                    ggAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGATGATCGAGGTGGCACAgc                                                                                                        >  2:325833/1‑149 (MQ=255)
                                     tctgACCGGCTTGTCCAGGCGTGTCGAAGCACCCGTAAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGTATGTATTAAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgc                                                                                                        <  2:184478/145‑1 (MQ=255)
                                     gAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgcg                                                                                                       >  7:205979/1‑149 (MQ=255)
                                     gAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgcg                                                                                                       >  1:173248/1‑149 (MQ=255)
                                     gAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgc                                                                                                        >  7:63038/1‑148 (MQ=255)
                                     gAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGAAGATCTAGGTCGaacacc                                                                                                        >  4:360275/1‑146 (MQ=255)
                                       aGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTc                                                                                                                                       >  8:48071/1‑115 (MQ=255)
                                       aGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTc                                                                                                                                       <  7:48071/115‑1 (MQ=255)
                                         aCCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgc                                                                                                        >  1:184478/1‑144 (MQ=255)
                                          ccGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGa                                                                                                  <  4:108685/149‑1 (MQ=255)
                                           cGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATc                                                                                                                     >  5:135929/1‑129 (MQ=255)
                                           cGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATc                                                                                                                     <  6:135929/129‑1 (MQ=255)
                                           cGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTAAACTGGCGGTAATGTACAAAAACATGGATGCATTCAAATCAATCAACAAATCAAACGGGCAATACTAAGGCGATAAATTGAATATCCAGGTGGCAAAgccattaac                                                                                                 >  2:331368/1‑142 (MQ=255)
                                              cTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGAcc                                                                                                >  5:201219/1‑147 (MQ=255)
                                                 gTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcagc                                                                                            <  7:210690/148‑1 (MQ=255)
                                                 gTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcagc                                                                                            <  6:201219/148‑1 (MQ=255)
                                                 gTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcag                                                                                             >  5:354019/1‑147 (MQ=255)
                                                  tCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcagca                                                                                           <  6:136565/148‑1 (MQ=255)
                                                     aGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACCTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGgtc                                                                                       >  7:198189/1‑149 (MQ=255)
                                                      ggCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGgtcc                                                                                      <  5:118066/149‑1 (MQ=255)
                                                              gcagcaGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGa                                                                                                  <  5:128329/129‑1 (MQ=255)
                                                              gcagcaGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGa                                                                                                  >  6:128329/1‑129 (MQ=255)
                                                              gcagcaGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCt                                                                              >  2:126892/1‑149 (MQ=255)
                                                                agcagcCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGa                                                                            >  6:360034/1‑149 (MQ=255)
                                                                     ccttCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACcg                                                                       <  5:92602/146‑1 (MQ=255)
                                                                     ccGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACcg                                                                       >  2:140165/1‑149 (MQ=255)
                                                                       gTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACcgtc                                                                     <  1:294613/149‑1 (MQ=255)
                                                                              tAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTa                                                              <  5:360034/149‑1 (MQ=255)
                                                                              tAAAGTTCAACTGGAGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCGATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTa                                                              <  4:98119/149‑1 (MQ=255)
                                                                               aaaGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcagc                                                                                            <  1:53287/118‑1 (MQ=255)
                                                                               aaaGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcagc                                                                                            >  2:53287/1‑118 (MQ=255)
                                                                                aaGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATg                                                            <  1:126892/149‑1 (MQ=255)
                                                                                  gTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATggg                                                          >  2:73045/1‑149 (MQ=255)
                                                                                    tCAACTTGCGGTAATGTACATAGACCTTGATGGATTCAAATCCATCAACGACTCCATCGTGCACTACTTCGGCGATCAAGTGCTGATCGATGTGGCACAGCGATTGACACAGCAGGTCCTTCCTGCTGATACCGTCGTTCGTATGggc                                                         <  3:327843/148‑1 (MQ=255)
                                                                                      aaCTGGCTGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGggcggc                                                      <  6:354019/149‑1 (MQ=255)
                                                                                      aaCTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGggcgg                                                       <  1:73045/148‑1 (MQ=255)
                                                                                           gCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGa                                                 <  1:140165/149‑1 (MQ=255)
                                                                                           gCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGa                                                 <  4:263240/149‑1 (MQ=255)
                                                                                           gCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGa                                                 <  4:258308/149‑1 (MQ=255)
                                                                                           gCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACCTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGa                                                 <  8:198189/149‑1 (MQ=255)
                                                                                             ggTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCAGGCACTACGATCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAg                                                >  1:106535/1‑149 (MQ=255)
                                                                                             ggTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCAGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGt                                               >  7:268539/1‑149 (MQ=255)
                                                                                             ggTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCAGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGt                                               >  7:263574/1‑149 (MQ=255)
                                                                                              gTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGtt                                              >  7:226288/1‑149 (MQ=255)
                                                                                               tCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGttt                                             <  2:173248/149‑1 (MQ=255)
                                                                                               tCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGttt                                             <  7:125768/149‑1 (MQ=255)
                                                                                                 aTGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGt                                           <  8:63038/149‑1 (MQ=255)
                                                                                                 aTGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGt                                           <  1:331368/149‑1 (MQ=255)
                                                                                                     aCATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGtt                                       >  7:44185/1‑149 (MQ=255)
                                                                                                     aCATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCATGCGTATGGGCGGCGATGAGTTTGTCGtt                                       >  3:366773/1‑149 (MQ=255)
                                                                                                          gACCTGGATGGATTCAAATCGATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTAGGGGCGGCGATGAGTTTGTCGTTctgct                                  <  3:360275/149‑1 (MQ=255)
                                                                                                          gACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTctgc                                   <  6:364260/148‑1 (MQ=255)
                                                                                                                gatggatTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACg                              <  8:205979/147‑1 (MQ=255)
                                                                                                                gatggatTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCATGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACgg                             <  4:366773/148‑1 (MQ=255)
                                                                                                                 atggatTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGAc                           >  2:186348/1‑149 (MQ=255)
                                                                                                                     atTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTg                                            <  7:185225/128‑1 (MQ=255)
                                                                                                                     atTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTg                                            >  8:185225/1‑128 (MQ=255)
                                                                                                                      ttCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCt                      >  7:50715/1‑149 (MQ=255)
                                                                                                                      ttCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGAAGATCTAGGTGGAACACCTATAAAAAAAAAAAATAAGACCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACATGCt                      >  8:211270/1‑149 (MQ=255)
                                                                                                                       tCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCtt                     >  5:197138/1‑149 (MQ=255)
                                                                                                                        cAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCtt                     <  1:325833/148‑1 (MQ=255)
                                                                                                                         aaaTCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGAAGATCTAGGTGGAACAGCGATTGACACATCAAGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGTCCTGCTTCg                   >  6:349740/1‑149 (MQ=255)
                                                                                                                          aaTCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGt                  >  5:125257/1‑149 (MQ=255)
                                                                                                                            tCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGa                >  1:21933/1‑149 (MQ=255)
                                                                                                                               aTCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACg              >  6:294413/1‑148 (MQ=255)
                                                                                                                                tCAACGACTCCATCAGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGaa            <  2:106535/149‑1 (MQ=255)
                                                                                                                                tCAACGACTCCATCAGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGaa            <  8:263574/149‑1 (MQ=255)
                                                                                                                                tCAACGACTCCATCAGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGaa            <  8:268539/149‑1 (MQ=255)
                                                                                                                                   aCGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGAAGAc         >  6:192596/1‑149 (MQ=255)
                                                                                                                                     gACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGAAGACGc       <  8:50715/149‑1 (MQ=255)
                                                                                                                                      aCTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGAAGACGCt      <  6:263698/149‑1 (MQ=255)
                                                                                                                                        tCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGAAGACGCTGt    >  5:238065/1‑149 (MQ=255)
                                                                                                                                          cATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGAAGACGCTGTcc  >  4:48508/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTCGTGACGAAGACGCTGTCC  >  NC_002947/3755492‑3755778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: