Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 86,245 A→C 5.1% M158R (ATG→AGG)  PP_0075 ← choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294786,2450AC 5.1% 371.4 / 7.0 126M158R (ATG→AGG) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (61/58);  new base C (4/3);  total (65/61)
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.99e-01

Only 100 of 127 total aligned reads displayed.

GCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTC  >  NC_002947/86097‑86394
                                                                                                                                                     |                                                                                                                                                     
gCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGCAAGGTAGTTGACCTCGCCTGCGTGACTGGCCGGGGATCGCACGAGCAGCAGCGCGGTGGCCTGGCTGGGCAGGGCCAt                                                                                                                                                      >  1:1141744/1‑150 (MQ=255)
      ccAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcag                                                                                                                                                <  2:1055120/150‑1 (MQ=255)
             gggCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTg                                                                                                                                         >  1:1028983/1‑150 (MQ=255)
             gggCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTg                                                                                                                                         <  1:591121/150‑1 (MQ=255)
              ggCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTgg                                                                                                                                        <  2:1773936/150‑1 (MQ=255)
                 gcCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTa                                                                                                                                     >  2:269706/1‑150 (MQ=255)
                  ccAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAg                                                                                                                                    >  2:134425/1‑150 (MQ=255)
                   cAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGACGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGc                                                                                                                                   >  2:42372/1‑150 (MQ=255)
                   cAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGc                                                                                                                                   >  1:1765137/1‑150 (MQ=255)
                    ggggTTTGAAGGCCGCACCGATCACGCCGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCt                                                                                                                                  <  2:849589/147‑1 (MQ=255)
                    agagTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagca                                                                                                                                              >  1:2159974/1‑138 (MQ=255)
                    agagTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagca                                                                                                                                              <  2:2159974/138‑1 (MQ=255)
                        tttGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTTGCCTTGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGcc                                                                                                                              <  1:46727/150‑1 (MQ=255)
                        tttGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGcc                                                                                                                              >  1:1037066/1‑150 (MQ=255)
                          tGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAgg                                                                                                                                            >  2:417585/1‑134 (MQ=255)
                          tGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAgg                                                                                                                                            <  1:417585/134‑1 (MQ=255)
                           gAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCaa                                                                                                                           >  2:919560/1‑150 (MQ=255)
                            aaGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                                          <  2:1796154/150‑1 (MQ=255)
                            aaGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                                          >  1:52451/1‑150 (MQ=255)
                            aaGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                                          <  1:269706/150‑1 (MQ=255)
                                 cGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcag                                                                                                                                                >  2:1811837/1‑123 (MQ=255)
                                 cGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAAt                                                                                                                     <  1:919560/150‑1 (MQ=255)
                                   cACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGa                                                                                                                   <  1:1723788/150‑1 (MQ=255)
                                    aCCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGa                                                                                                                   <  1:124740/149‑1 (MQ=255)
                                    aCCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGa                                                                                                                   >  2:124740/1‑149 (MQ=255)
                                       gATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGAt                                                                                                               >  1:289440/1‑150 (MQ=255)
                                       gATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGAt                                                                                                               <  2:1141744/150‑1 (MQ=255)
                                         tCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATgg                                                                                                             <  2:1151877/150‑1 (MQ=255)
                                           acgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcg                                                                                                           >  1:1345777/1‑150 (MQ=255)
                                           acgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcg                                                                                                           >  2:1727133/1‑150 (MQ=255)
                                             gacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGGGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGcc                                                                                                                              <  1:1134572/129‑1 (MQ=255)
                                             gacgGGGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTTGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcggc                                                                                                         >  1:1757266/1‑150 (MQ=255)
                                              acgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCa                                                                                                        >  2:1020919/1‑150 (MQ=255)
                                                     ccAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCg                                                                                                 >  2:1086106/1‑150 (MQ=255)
                                                      cAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGt                                                                                                >  2:1927166/1‑150 (MQ=255)
                                                       aGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTg                                                                                               >  1:2270326/1‑150 (MQ=255)
                                                         gCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGaa                                                                                             <  2:52106/150‑1 (MQ=255)
                                                            cagcagCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGcc                                                                                          <  1:1305193/150‑1 (MQ=255)
                                                             agcagcAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGaa                                                                                             <  1:1557928/146‑1 (MQ=255)
                                                             agcagcAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGaa                                                                                             >  2:1557928/1‑146 (MQ=255)
                                                              gcagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                                        <  2:1037066/150‑1 (MQ=255)
                                                              gcagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCAAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                                        <  2:96146/150‑1 (MQ=255)
                                                                   aaggaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGcacc                                                                                   >  2:1290470/1‑150 (MQ=255)
                                                                   aaggaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGcacc                                                                                   >  1:102269/1‑150 (MQ=255)
                                                                          ggtaggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGGAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGaa                                                                             <  1:182723/150‑1 (MQ=255)
                                                                          ggtaggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGaa                                                                             >  2:2315904/1‑150 (MQ=255)
                                                                           gtaggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAAt                                                                            <  2:293802/150‑1 (MQ=255)
                                                                           gtaggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAAt                                                                            <  1:134425/150‑1 (MQ=255)
                                                                             aggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATgg                                                                          <  2:2270326/150‑1 (MQ=255)
                                                                               gtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGGCTGGCTGGTCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                                        >  2:798525/1‑150 (MQ=255)
                                                                               gtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                                        <  2:558389/150‑1 (MQ=255)
                                                                               gtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                                        >  2:2271655/1‑150 (MQ=255)
                                                                                  gACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGACTGGCTGGGAAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGCTTGGgagag                                                                     >  1:103210/1‑149 (MQ=255)
                                                                                      tCGCCTGCGTGGCTGGCCAGGTCCTTCACGGTCATCAGCGCTGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAg                                                                 <  1:1290470/150‑1 (MQ=255)
                                                                                        gCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGtt                                                               <  2:1757266/150‑1 (MQ=255)
                                                                                         ccTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                                        >  1:2201638/1‑140 (MQ=255)
                                                                                           tGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCca                                                                                                                                                 <  1:524691/65‑1 (MQ=255)
                                                                                           tGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCca                                                                                                                                                 >  2:524691/1‑65 (MQ=255)
                                                                                              gTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAg                                                         >  1:1207059/1‑150 (MQ=255)
                                                                                               tggctggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAgc                                                        <  2:2146471/150‑1 (MQ=255)
                                                                                                 gctggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGCCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAgcgc                                                      >  1:1587206/1‑150 (MQ=255)
                                                                                                   tggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                                      >  1:1988930/1‑132 (MQ=255)
                                                                                                   tggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                                      <  2:1988930/132‑1 (MQ=255)
                                                                                                   tggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAg                                                                 >  2:120714/1‑137 (MQ=255)
                                                                                                      ccAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGa                                                 >  2:1350777/1‑150 (MQ=255)
                                                                                                      ccAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGa                                                 >  2:866205/1‑150 (MQ=255)
                                                                                                       cAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGaa                                                >  1:2017385/1‑150 (MQ=255)
                                                                                                       cAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGaa                                                >  2:2222392/1‑150 (MQ=255)
                                                                                                       cAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGaa                                                >  1:535603/1‑150 (MQ=255)
                                                                                                         ggTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTTAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgc                                              >  1:2331764/1‑150 (MQ=255)
                                                                                                         ggTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgc                                              <  1:1350777/150‑1 (MQ=255)
                                                                                                           ctcTGCACGGTCATCAGCGCGGTGGCCGGGCTGGGCAGGGCCATGACCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcgc                                            <  2:1949285/148‑1 (MQ=255)
                                                                                                             cggCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCa                                           >  2:2154561/3‑149 (MQ=255)
                                                                                                             cTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcgc                                            <  2:456046/148‑1 (MQ=255)
                                                                                                             cTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcgc                                            >  1:456046/1‑148 (MQ=255)
                                                                                                             cTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCaa                                          >  2:1325165/1‑150 (MQ=255)
                                                                                                             cTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCa                                           <  1:2154561/149‑1 (MQ=255)
                                                                                                              tGCACGGTCATCAGCGCGGTGGCCTGGCGGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCaag                                         >  1:1949285/1‑149 (MQ=255)
                                                                                                                aaCGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAg                                       >  2:12023/2‑150 (MQ=255)
                                                                                                                   ggTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGccc                                    <  2:620106/150‑1 (MQ=255)
                                                                                                                   ggTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGccc                                    <  2:1345777/150‑1 (MQ=255)
                                                                                                                        tcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATgg                                                                          >  1:2312546/1‑107 (MQ=255)
                                                                                                                        tcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagca                               >  2:565971/1‑150 (MQ=255)
                                                                                                                              cgGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCt                         <  2:289440/150‑1 (MQ=255)
                                                                                                                                  caacctgctggGCATGGCCCTGCCCAGCCGGTGGTGTAGCTGGCCCAAGGCAATGCCGATGGCGGCACCCAGCGTGAAGCCGAGCACCCCGGACTGGGAGACGAACCAGACCCCCACGCGGACGCGCAATAGCCCGAGCAGCAGCCGGaa                     >  2:1955254/7‑150 (MQ=255)
                                                                                                                                   aacctgctggGCATGGCCCTGCCCAGCCGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGaaa                    >  1:1634291/6‑150 (MQ=255)
                                                                                                                                    ccTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGAAGTGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGaaaa                   <  1:531728/150‑1 (MQ=255)
                                                                                                                                     ctggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAAt                  >  1:1186987/1‑150 (MQ=255)
                                                                                                                                      tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATg                 >  1:1306173/1‑150 (MQ=255)
                                                                                                                                      tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATg                 >  1:16470/1‑150 (MQ=255)
                                                                                                                                      tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAATTGACCAGCGCGCCGAACCGCAATAGCCCGAGCAGCCGCCGGcaagtc                 >  2:2202265/1‑147 (MQ=255)
                                                                                                                                      ctgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATg                 >  2:874969/3‑150 (MQ=255)
                                                                                                                                       ggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgc                <  1:1727133/150‑1 (MQ=255)
                                                                                                                                             gCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTTACATCAGCCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGa          <  1:2222392/150‑1 (MQ=255)
                                                                                                                                              catggCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgc                                              <  1:1863435/110‑1 (MQ=255)
                                                                                                                                              catggCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgc                                              >  2:1863435/4‑113 (MQ=255)
                                                                                                                                                 ggCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGaa      <  2:1186987/150‑1 (MQ=255)
                                                                                                                                                 ggCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGAAGATGGAGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGaa      >  1:658661/1‑150 (MQ=255)
                                                                                                                                                  gCCCTGCCCAGCCGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAg     <  2:1634291/150‑1 (MQ=255)
                                                                                                                                                     aTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTc  >  1:1913947/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTC  >  NC_002947/86097‑86394

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: