Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,755,655 C→T 8.2% F351F (TTC→TTT 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,6550CT 8.2% 196.9 / 5.3 70F351F (TTC→TTTPP_3319GGDEF domain‑containing sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base T (3/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTC  >  NC_002947/3755512‑3755803
                                                                                                                                               |                                                                                                                                                    
cAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATcc                                                                                                                                                >  1:185312/1‑150 (MQ=255)
cAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATcc                                                                                                                                                >  2:2360203/1‑150 (MQ=255)
  gTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCAt                                                                                                                                              <  1:1802076/150‑1 (MQ=255)
  gTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCAt                                                                                                                                              <  2:2120101/150‑1 (MQ=255)
    tGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTTGAATCCATcc                                                                                                                                            <  1:2392569/150‑2 (MQ=255)
    tGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCACAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAatcaatca                                                                                                                                            >  1:762716/1‑150 (MQ=255)
      aGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAAc                                                                                                                                          <  1:2274386/150‑1 (MQ=255)
       gTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTTGATGGATTTGAATCCATccagg                                                                                                                                         <  1:833414/150‑5 (MQ=255)
          tCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACt                                                                                                                                      >  2:118715/1‑150 (MQ=255)
                 cagccaTGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCgg                                                                                                                               >  1:1696894/1‑150 (MQ=255)
                           gggCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACtt                                                                                                                      <  1:430238/149‑1 (MQ=255)
                           gggCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACtt                                                                                                                      >  2:430238/1‑149 (MQ=255)
                             gCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAAc                                                                                                                                          <  1:2362591/127‑1 (MQ=255)
                             gCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAAc                                                                                                                                          >  2:2362591/1‑127 (MQ=255)
                               tGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCg                                                                                                                 >  1:1578598/1‑150 (MQ=255)
                                gACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCg                                                                                                                                <  1:187971/134‑1 (MQ=255)
                                gACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCg                                                                                                                                >  2:187971/1‑134 (MQ=255)
                                  cGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATc                                                                                                              <  2:79602/150‑1 (MQ=255)
                                           ttAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAAc                                                                                                                                          <  1:667797/113‑1 (MQ=255)
                                           ttAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAAc                                                                                                                                          >  2:667797/1‑113 (MQ=255)
                                           ttAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGa                                                                                                     >  2:1007872/1‑150 (MQ=255)
                                             aaCCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATc                                                                                                   >  2:1425561/1‑150 (MQ=255)
                                                   cTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTg                                                                                             <  2:1696894/150‑1 (MQ=255)
                                                          tGGAGGAGCGGGTTGTCTAGGCGTGGCGCGGCGGCCGTCAACGTAAAGTTCGACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgc                                                                                      >  2:1298058/1‑150 (MQ=255)
                                                           ggAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAANGTAAAGNTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAgcg                                                                                     >  1:1649713/1‑150 (MQ=255)
                                                            gAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGt                                                                                              >  1:238334/1‑140 (MQ=255)
                                                            gAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGt                                                                                              <  2:238334/140‑1 (MQ=255)
                                                                    gCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATATACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCca                                                                            >  2:766805/1‑150 (MQ=255)
                                                                         tCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCa                                                                                                                                               >  1:941583/1‑78 (MQ=255)
                                                                         tCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCa                                                                                                                                               <  2:941583/78‑1 (MQ=255)
                                                                         tCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGGCCCAGCAgg                                                                       >  2:1923342/1‑150 (MQ=255)
                                                                         tCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcagc                                                                          <  2:1494200/147‑1 (MQ=255)
                                                                         tCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCcagc                                                                          >  1:1494200/1‑147 (MQ=255)
                                                                         tCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCAATAGTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAgg                                                                       >  1:147955/1‑150 (MQ=255)
                                                                           cAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGgtc                                                                     >  1:993235/1‑150 (MQ=255)
                                                                            aGGCGAGTAGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGgtcc                                                                    >  2:1956402/1‑150 (MQ=255)
                                                                              gCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGgtccgt                                                                  >  1:433777/1‑150 (MQ=255)
                                                                                  gtCGCAGCAGCTCTCACCGTTAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGCTTCAAATCCATCAACGACTCCATCGGGTCCTACTTCGGCGATCAACTGATGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCgg                                                              <  2:164894/150‑1 (MQ=255)
                                                                                   tCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGGCCCAGCAGGTCCGTCCGGc                                                             <  1:1923342/150‑1 (MQ=255)
                                                                                   tCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGc                                                             >  1:477038/1‑150 (MQ=255)
                                                                                       agcagcCGTCAACGTAAAGGTCAACTGGCGGTCATGTACGTCGACCTGGCTGGATTCAAATCCACCAACGACTCCATCGGGCAGTACTTCGGCGATCAAGTGCTGATGGAGGTGGCACAGCGCTTGACCCAACAGGTCCGTCCGCCTGAt                                                         >  2:985344/1‑150 (MQ=255)
                                                                                             cGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGACGGCTTCCGCTCCATCAACGACTCCATCGGGCACTACTTCGGCGATCACGTGCTGATCNAGGTGGNACAGCGCTTGACNCAGCAGGTCCGTCCGGCTGATACcgtc                                                   <  1:1955278/150‑1 (MQ=255)
                                                                                                       aaGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCgg                                                                                                                               >  1:73467/1‑64 (MQ=255)
                                                                                                       aaGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCgg                                                                                                                               <  2:73467/64‑1 (MQ=255)
                                                                                                       aaGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATgg                                         <  2:762716/150‑1 (MQ=255)
                                                                                                         gTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCtt                                                                                  >  2:1442500/1‑107 (MQ=255)
                                                                                                         gTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCtt                                                                                  <  1:1442500/107‑1 (MQ=255)
                                                                                                          ttCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGggcg                                      <  2:1578598/150‑1 (MQ=255)
                                                                                                             aaCTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGgcggcg                                   >  2:2305171/1‑150 (MQ=255)
                                                                                                                tGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATg                                <  1:2360203/150‑1 (MQ=255)
                                                                                                                 ggCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGa                               <  1:118715/150‑1 (MQ=255)
                                                                                                                      tCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTg                          >  1:956049/1‑150 (MQ=255)
                                                                                                                        aTGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTc                        <  2:477038/150‑1 (MQ=255)
                                                                                                                          gTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAgg                                                                       <  1:1120179/101‑1 (MQ=255)
                                                                                                                          gTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAgg                                                                       >  2:1120179/1‑101 (MQ=255)
                                                                                                                          gTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGt                      >  2:1558384/1‑150 (MQ=255)
                                                                                                                                cGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTctgct                <  2:433777/150‑1 (MQ=255)
                                                                                                                                 gACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTctgctg               >  2:1808219/1‑150 (MQ=255)
                                                                                                                                  aCCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTAt                                           <  2:2363429/121‑1 (MQ=255)
                                                                                                                                  aCCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTAt                                           >  1:2363429/1‑121 (MQ=255)
                                                                                                                                    gttgatggatTTGAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGttt                           <  2:226209/132‑1 (MQ=255)
                                                                                                                                    gttgatggatTTGAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGttt                           >  1:226209/4‑135 (MQ=255)
                                                                                                                                    cTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACg            <  1:1298058/150‑1 (MQ=255)
                                                                                                                                     ttgatggatTTGAATCCATCAACGACTCCATCGGGCACTACGTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTGGGGGCGGCGAGGAGTTTGTCGTTCTGCTGACgg           >  1:858194/3‑150 (MQ=255)
                                                                                                                                      tgatggatTTGAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGa          >  2:133096/2‑150 (MQ=255)
                                                                                                                                             ttCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCGAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGAc         <  1:1936178/144‑1 (MQ=255)
                                                                                                                                             ttCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCGAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGAc         >  2:1936178/1‑144 (MQ=255)
                                                                                                                                             ttCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCCCGATACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCtt   <  1:1956402/150‑1 (MQ=255)
                                                                                                                                              tCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTc  <  2:1649713/150‑1 (MQ=255)
                                                                                                                                              tCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTc  <  2:185312/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CAGTTGAGTATCAATGCCAGCCATGACGGGCTGACGGGTTTGCTTAACCGCCTGTCTCTGGAAGACCGGCTTGTCCAGGCGTGTCGCAGCAGCCGTCAACGTAAAGTTCAACTGGCGGTCATGTACATCGACCTGGATGGATTCAAATCCATCAACGACTCCATCGGGCACTACTTCGGCGATCAAGTGCTGATCGAGGTGGCACAGCGCTTGACCCAGCAGGTCCGTCCGGCTGATACCGTCGTGCGTATGGGCGGCGATGAGTTTGTCGTTCTGCTGACGGACCTGCTTC  >  NC_002947/3755512‑3755803

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