Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,633,178 C→T 6.4% L203L (CTG→TTG)  PP_3201 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,633,1780CT 6.4% 326.1 / 7.4 109L203L (CTG→TTG) PP_3201BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (51/51);  new base T (4/3);  total (55/54)
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

Only 100 of 116 total aligned reads displayed.

AGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGC  >  NC_002947/3633029‑3633325
                                                                                                                                                     |                                                                                                                                                    
aGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGccc                                                                                                                                                      <  2:1077165/150‑1 (MQ=255)
  ccgatcttgGCCCTGAATTTCACCGATGCCCGCATCGGCCCGGTCGTCGGTTCCTACGGGGTGGCCATGCGCACCGACGTCGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACCGGTGTGGACTGCACCTGTATGCCCTg                                                                                                                                                    <  2:901189/144‑1 (MQ=255)
  ccGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTg                                                                                                                                                    <  2:2777185/150‑1 (MQ=255)
  ccGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTg                                                                                                                                                    >  2:320579/1‑150 (MQ=255)
   cGTTGTTGGCCCTGAATTTCACAGATGCCCGCAACGCCCTCGTCGTAGGTGCCTACGGGGTTGCTATCCGCCCCGTCGACGCAGGGCGCAGCTCGCAGTCAGAGATCGGCCGACTGACGCCCGGGCGTGGACTGCACCTGTATGCACTgg                                                                                                                                                   >  2:228735/1‑150 (MQ=255)
          ggCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCATTCGGAGATCGGCCGACTGCCGAACGGGCGTGGAGTGCACCTGTATGCCCTggagatcgg                                                                                                                                            >  1:901189/1‑143 (MQ=255)
           gCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCaa                                                                                                                                           <  2:534611/150‑1 (MQ=255)
            cccTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAg                                                                                                                                          <  1:1405349/150‑1 (MQ=255)
               tGAATTTCCCCGATACGGTGAAGGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCa                                                                                                                                       >  2:128908/1‑150 (MQ=255)
                   tttCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCa                                                                                                                                   >  2:418036/1‑150 (MQ=255)
                    ttCCCCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAg                                                                                                                                  <  1:1943624/150‑1 (MQ=255)
                    ttCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCat                                                                                                                                  <  2:2109557/150‑2 (MQ=255)
                      cACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgt                                                                                                                                <  1:530488/150‑1 (MQ=255)
                       aCCGATGCCCGCAACTTCCTGGTAGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGAGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgtg                                                                                                                               >  2:1996735/1‑150 (MQ=255)
                       aCCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgtg                                                                                                                               >  1:2870612/1‑150 (MQ=255)
                        ccGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgtgt                                                                                                                              <  1:418036/150‑1 (MQ=255)
                          gATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGCGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAGCCAGgtgggg                                                                                                                            >  1:2821549/1‑147 (MQ=255)
                            tGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGCCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACAt                                                                                                                          >  2:481723/1‑150 (MQ=255)
                                   aaCGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTg                                                                                                                                                    <  1:1518849/117‑1 (MQ=255)
                                   aaCGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGTTGTACATCGCCGGc                                                                                                                  >  2:351127/1‑150 (MQ=255)
                                      ggCCTGGTCGTCGGGGCCTACGGGGTGGCGAGGCGCACCGACGACGGCGGGCGCATTTTGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAg                                                                                                               <  2:416242/150‑1 (MQ=255)
                                       gCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAg                                                                                                                                  <  1:178134/131‑1 (MQ=255)
                                       gCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAg                                                                                                                                  >  2:178134/1‑131 (MQ=255)
                                       gCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGc                                                                                                              >  1:2900139/1‑150 (MQ=255)
                                        ccTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGAAGACGGCAGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCa                                                                                                             >  1:714640/1‑150 (MQ=255)
                                         cTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAg                                                                                                            >  1:1879982/1‑150 (MQ=255)
                                                tcgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCctgc                                                                                                     >  2:108307/1‑150 (MQ=255)
                                                 cgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGTTGTACATCGCCGGCGAGCAGGGCctgct                                                                                                    <  1:351127/150‑1 (MQ=255)
                                                  ggTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCctgctg                                                                                                   >  1:103640/1‑150 (MQ=255)
                                                   gTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATACCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc                                                                                                         <  2:1277948/143‑1 (MQ=255)
                                                   gTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATACCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc                                                                                                         >  1:1277948/1‑143 (MQ=255)
                                                          cGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtcgctc                                                                                           <  1:1213114/150‑1 (MQ=255)
                                                             ggTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTcgcg                                                                                        >  2:2023190/1‑150 (MQ=255)
                                                                  ccATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgtg                                                                                                                               >  2:1459699/1‑107 (MQ=255)
                                                                  ccATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgtg                                                                                                                               <  1:1459699/107‑1 (MQ=255)
                                                                  ccATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATg                                                                                   >  1:1229816/1‑150 (MQ=255)
                                                                     tGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCg                                                                                <  2:1879982/150‑1 (MQ=255)
                                                                      gcgcACGGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGGGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCTGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCgg                                                                               <  2:886257/150‑1 (MQ=255)
                                                                      gcgcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCgg                                                                               >  1:666443/1‑150 (MQ=255)
                                                                       cgcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGt                                                                              <  1:2023190/150‑1 (MQ=255)
                                                                        gcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtc                                                                                               >  2:1436549/1‑132 (MQ=255)
                                                                        gcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtc                                                                                               <  1:1436549/132‑1 (MQ=255)
                                                                         cACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGACCGGCCGAGTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGa                                                                            <  1:2475691/150‑1 (MQ=255)
                                                                                 gcggcggGCGCAGCTGGCAGTCAGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACg                                                                    <  1:228735/149‑1 (MQ=255)
                                                                                  cggcggGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGa                                                                   <  2:2623566/150‑1 (MQ=255)
                                                                                  cggcggGCGCAGCTGGCAGTCGGAGATCGCCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACCTCAGCCGGGGCTAGCGGATCCTCCTCCGCACTCGTTATGGCGGTGAGCATTACGa                                                                   <  1:386088/150‑1 (MQ=255)
                                                                                   ggcgggcgCAGCTGGCAGTCGGAGATCGGCCGGCTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATCACGAg                                                                  >  1:1514012/1‑150 (MQ=255)
                                                                                     cgggcgCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGgc                                                                >  2:1079013/1‑150 (MQ=255)
                                                                                      gggcgCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGgcg                                                               >  1:2524663/1‑150 (MQ=255)
                                                                                         cgcAGCTGGCAGTCGGAGATCGGCCGGCTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATCACGAGGCGCTg                                                            <  2:1514012/150‑1 (MQ=255)
                                                                                          gcAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGc                                                           >  1:76810/1‑150 (MQ=255)
                                                                                                ggCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGc                                                     <  2:1928341/150‑1 (MQ=255)
                                                                                                ggCAGTCGGAGACCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAACGCAACCAGGTGTACTTCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGc                                                     <  2:713813/150‑1 (MQ=255)
                                                                                                 gCAGTCGGAGATCGGCCGACTGCCGAACGGTCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGcc                                                    <  2:714640/150‑1 (MQ=255)
                                                                                                 gCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGcc                                                    >  1:1740270/1‑150 (MQ=255)
                                                                                                     tCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTAc                                                <  1:1249462/150‑1 (MQ=255)
                                                                                                       ggAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtc                                                                                               <  2:258582/101‑1 (MQ=255)
                                                                                                        gagaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGAt                                             <  1:2187911/150‑1 (MQ=255)
                                                                                                          gaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATCTCGGTTAGCATTACGAGGCGCTGCACGGGCCCTACGATgg                                           >  1:2039245/1‑150 (MQ=255)
                                                                                                            tCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                         <  1:1079013/150‑1 (MQ=255)
                                                                                                               gCCGACTGCCGAACGGTGGTGGACGTCACCTGGATGCCCGGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCt                                      >  2:41063/1‑150 (MQ=255)
                                                                                                               gCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCt                                      >  2:2360363/1‑150 (MQ=255)
                                                                                                                 cGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                         <  2:799002/145‑1 (MQ=255)
                                                                                                                 cGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                         >  1:799002/1‑145 (MQ=255)
                                                                                                                 cGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCttc                                    <  1:108307/150‑3 (MQ=255)
                                                                                                                  gACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGt                                                                                       >  1:2810076/1‑98 (MQ=255)
                                                                                                                  gACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGt                                                                                       <  2:2810076/98‑1 (MQ=255)
                                                                                                                     tGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATg                                                                                   <  1:222353/99‑1 (MQ=255)
                                                                                                                     tGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGt                                                                              <  1:408236/104‑1 (MQ=255)
                                                                                                                     tGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGt                                                                              >  2:408236/1‑104 (MQ=255)
                                                                                                                       ccGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTcgcc                              >  2:2750841/1‑150 (MQ=255)
                                                                                                                                 gTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCgtgt                    >  1:1736858/1‑150 (MQ=255)
                                                                                                                                    gACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGc                 <  2:1241499/150‑1 (MQ=255)
                                                                                                                                    gACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGAGTGATGGCGGTGAGCATTACTAGGCTCTGTACTGGCCGTACGATGGCAGCTACTTCGCCGCCACCGTGTTGc                 <  1:58342/150‑1 (MQ=255)
                                                                                                                                       ttcACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGcc                                                                                                        <  2:1952240/58‑1 (MQ=255)
                                                                                                                                       tGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGcc                                                                                                        >  1:1952240/1‑60 (MQ=255)
                                                                                                                                         cACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGAt            >  1:266224/1‑150 (MQ=255)
                                                                                                                                         cACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGAt            <  2:76810/150‑1 (MQ=255)
                                                                                                                                          aCCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATg           <  1:1983687/150‑1 (MQ=255)
                                                                                                                                          aCCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATg           >  2:1218062/1‑150 (MQ=255)
                                                                                                                                           cctggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATgg          >  2:14738/7‑150 (MQ=255)
                                                                                                                                           ccTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATgg          <  2:266224/150‑1 (MQ=255)
                                                                                                                                            ctggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                         >  2:2295662/6‑118 (MQ=38)
                                                                                                                                            ctggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                         <  1:2295662/113‑1 (MQ=38)
                                                                                                                                            cTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGa                                                                            >  2:2165300/1‑83 (MQ=255)
                                                                                                                                             tggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                         <  1:323473/113‑1 (MQ=38)
                                                                                                                                             tggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                         >  2:323473/5‑117 (MQ=38)
                                                                                                                                             tGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGcc        <  2:1736858/150‑1 (MQ=255)
                                                                                                                                               gttGCCTTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCa                          <  2:2197709/128‑1 (MQ=255)
                                                                                                                                               gttGCCTTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCa                          >  1:2197709/3‑130 (MQ=255)
                                                                                                                                                aTGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAg                                                                           >  2:2107984/1‑80 (MQ=255)
                                                                                                                                                aTGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAg                                                                           <  1:2107984/80‑1 (MQ=255)
                                                                                                                                                aTGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTc     >  1:1484192/1‑150 (MQ=255)
                                                                                                                                                 tGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTct    >  2:1952496/1‑150 (MQ=255)
                                                                                                                                                   cccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGa                                                                            <  1:1422429/76‑1 (MQ=255)
                                                                                                                                                   cccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGa                                                                            >  2:1422429/1‑76 (MQ=255)
                                                                                                                                                   cccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTctgc  >  1:2965670/1‑150 (MQ=255)
                                                                                                                                                   cccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTctgc  >  2:2493325/1‑150 (MQ=255)
                                                                                                                                                    ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAgg                                                                 >  1:1384980/1‑86 (MQ=255)
                                                                                                                                                    ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAgg                                                                 <  2:1384980/86‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTCTGC  >  NC_002947/3633029‑3633325

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