Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 86,245 A→C 5.9% 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.9% 200.1 / 2.7 67M158R (ATG→AGG) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base C (2/2);  total (34/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 = 9.34e-01

GGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTC  >  NC_002947/86096‑86394
                                                                                                                                                     |                                                                                                                                                      
ggCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATAACGATTACCGCAAGGGGCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCa                                                                                                                                                        <  2:30349/150‑1 (MQ=255)
 gCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCAt                                                                                                                                                       >  1:626668/1‑150 (MQ=255)
 gCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGCAATGCACGGTCTTCAGCGCGGTGGCCTGGCTGGGCAGGGCCAt                                                                                                                                                       >  2:533099/1‑150 (MQ=255)
  cTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATg                                                                                                                                                      <  2:817514/150‑1 (MQ=255)
           cGTGGGCGCCAGAGTTTGAAGGCCGCCCCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTGCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcagccag                                                                                                                                            <  1:112461/150‑4 (MQ=255)
           cGTGGGCGCCAGAGTTGGAAGGCCGCGCCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcagccag                                                                                                                                            <  2:265755/150‑4 (MQ=255)
            gTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGt                                                                                                                                           <  1:1086633/150‑1 (MQ=255)
            gTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGt                                                                                                                                           <  2:1099408/150‑1 (MQ=255)
                gcgcCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTggg                                                                                                                                       <  1:179068/150‑1 (MQ=255)
                    cAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGc                                                                                                                                   >  2:814980/1‑150 (MQ=255)
                     agagTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCt                                                                                                                                  >  2:649597/1‑150 (MQ=255)
                             aaGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                                          <  1:1034686/150‑1 (MQ=255)
                                gCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCa                                                                                                                       >  1:575603/1‑150 (MQ=255)
                                 ccGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCaa                                                                                                                      >  2:449028/1‑150 (MQ=255)
                                   gCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCca                                                                                                                                                  >  2:1006695/1‑121 (MQ=255)
                                   gCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCca                                                                                                                                                  <  1:1006695/121‑1 (MQ=255)
                                    cACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGa                                                                                                                   >  1:719473/1‑150 (MQ=255)
                                       cGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGa                                                                                                                <  1:649597/150‑1 (MQ=255)
                                       cGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCCGGGGTGGCGGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGa                                                                                                                <  1:577536/150‑1 (MQ=255)
                                        gATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGAt                                                                                                               <  1:281530/150‑1 (MQ=255)
                                          tCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATgg                                                                                                             >  1:1095126/1‑150 (MQ=255)
                                           cACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggc                                                                                                            <  1:1072464/150‑1 (MQ=255)
                                            acgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcg                                                                                                           <  1:1095086/150‑1 (MQ=255)
                                                     gCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTg                                                                                                                                         >  2:501334/1‑111 (MQ=255)
                                                     gCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTg                                                                                                                                         <  1:501334/111‑1 (MQ=255)
                                                      ccAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCg                                                                                                 <  1:814980/150‑1 (MQ=255)
                                                             cagcagCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGcc                                                                                          >  2:536916/1‑150 (MQ=255)
                                                             cagcagCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGcc                                                                                          <  1:102247/150‑1 (MQ=255)
                                                               gcagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCTAGCGTGAAGCCGa                                                                                        <  1:275476/150‑1 (MQ=255)
                                                               gcagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                                        <  2:39230/150‑1 (MQ=255)
                                                                       gaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGcc                                                                                                                                                    <  1:385283/83‑1 (MQ=255)
                                                                       gaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGcc                                                                                                                                                    >  2:385283/1‑83 (MQ=255)
                                                                       gaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACg                                                                                >  1:508569/1‑150 (MQ=255)
                                                                          ggtaggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGaa                                                                             <  2:1019297/150‑1 (MQ=255)
                                                                             aggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATgg                                                                          <  1:449028/150‑1 (MQ=255)
                                                                               gtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTggg                                                                                                                                       >  2:670408/1‑87 (MQ=255)
                                                                               gtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTggg                                                                                                                                       <  1:670408/87‑1 (MQ=255)
                                                                                  gACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGACCTGGCTGGGCAGGTCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCCCGGAATGGGAGAc                                                                     >  1:382282/1‑150 (MQ=255)
                                                                                    ccTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGa                                                                   >  2:177251/1‑150 (MQ=255)
                                                                                     cTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGaa                                                                  >  2:239663/1‑150 (MQ=255)
                                                                                      tCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAg                                                                                            >  1:1105036/1‑123 (MQ=255)
                                                                                      tCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAg                                                                                            <  2:1105036/123‑1 (MQ=255)
                                                                                      tCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAg                                                                 >  1:355474/1‑150 (MQ=255)
                                                                                          cTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGa                                                             >  2:71361/1‑150 (MQ=255)
                                                                                          cTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGa                                                             >  2:273368/1‑150 (MQ=255)
                                                                                                   tggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGGCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAgcgcgc                                                    >  1:868863/1‑150 (MQ=255)
                                                                                                   tggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAgcgcgc                                                    >  1:242518/1‑150 (MQ=255)
                                                                                                   tggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAgcgcgc                                                    >  1:86564/1‑150 (MQ=255)
                                                                                                       cAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGTAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGaa                                                <  1:1090559/150‑1 (MQ=255)
                                                                                                        aGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAg                                               >  2:1065708/1‑150 (MQ=255)
                                                                                                              tGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAAt                                         >  1:581143/1‑150 (MQ=255)
                                                                                                               gCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATa                                        <  1:1054064/150‑1 (MQ=255)
                                                                                                                 aCGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGc                                      >  2:198071/1‑150 (MQ=255)
                                                                                                                    gTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCg                                   >  1:374219/1‑150 (MQ=255)
                                                                                                                       atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGaa                                                <  2:1001825/134‑1 (MQ=255)
                                                                                                                       atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGaa                                                >  1:1001825/1‑134 (MQ=255)
                                                                                                                        tcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGCGACGAAGTTGACCAGCGCACCGAAGCGCAATAGCCCGagca                               >  2:559452/1‑150 (MQ=255)
                                                                                                                           gcgcGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcagca                            <  1:239663/150‑1 (MQ=255)
                                                                                                                                 ccaacctgctggGCATGGCCCTGCCCAGCCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTgg                       >  2:198765/8‑150 (MQ=255)
                                                                                                                                    ccTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGaaaa                   <  1:536916/150‑1 (MQ=255)
                                                                                                                                     ctggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAAt                  <  1:533099/150‑1 (MQ=255)
                                                                                                                                      tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                                      >  1:657270/1‑97 (MQ=255)
                                                                                                                                      tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                                      <  2:657270/97‑1 (MQ=255)
                                                                                                                                      tggctggGCAGGACCATGCACAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATg                 >  2:786359/1‑150 (MQ=255)
                                                                                                                                      ctgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGAAAGCCGAATTATACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATg                 >  2:309258/3‑150 (MQ=255)
                                                                                                                                       ggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgc                <  2:86564/150‑1 (MQ=255)
                                                                                                                                         ctggGCAGGGCCATGCCCAGCAGGTTGGGTACCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAACTGGAAAATgccg              <  2:581143/150‑1 (MQ=255)
                                                                                                                                                  gCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAg     <  2:355474/150‑1 (MQ=255)
                                                                                                                                                     aTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTc  <  1:177251/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTC  >  NC_002947/86096‑86394

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: