Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 86,245 A→C 7.2% 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 7.2% 216.1 / 7.4 69M158R (ATG→AGG) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base C (3/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

GCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAG  >  NC_002947/86097‑86391
                                                                                                                                                    |                                                                                                                                                   
gCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCa                                                                                                                                                     >  5:258847/1‑149 (MQ=255)
       cAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcag                                                                                                                                              <  1:189445/149‑1 (MQ=255)
        aGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagc                                                                                                                                             <  7:171360/149‑1 (MQ=255)
           gTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcagccag                                                                                                                                          <  2:216771/149‑4 (MQ=255)
            tGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCGTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGGCATGCCCAGCAGGt                                                                                                                                         <  4:123902/149‑1 (MQ=255)
                cgcCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTggg                                                                                                                                     <  6:34824/149‑1 (MQ=255)
                       gTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTgg                                                                                                                              <  6:258847/149‑1 (MQ=255)
                        tttGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGc                                                                                                                             >  8:7023/1‑149 (MQ=255)
                             aGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                                        <  7:7023/149‑1 (MQ=255)
                               gCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAgg                                                                                                                       >  7:186737/1‑148 (MQ=255)
                                 cGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCaa                                                                                                                    >  7:176430/1‑149 (MQ=255)
                                  gCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGTCCAGGTCCTGCACGGACATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑t                                                                                                                  >  8:64391/1‑149 (MQ=255)
                                     ccGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGAc                                                                                                               <  4:162663/149‑1 (MQ=255)
                                           acgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATggc                                                                                                         >  6:203141/1‑149 (MQ=255)
                                            cgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATggcg                                                                                                        <  8:32707/149‑1 (MQ=255)
                                             gacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATggcgg                                                                                                       >  3:214023/1‑149 (MQ=255)
                                                ggTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCATGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCAc                                                                                                    <  2:137460/149‑1 (MQ=255)
                                                ggTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCAc                                                                                                    <  2:201091/149‑1 (MQ=255)
                                                  tGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCGCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCAcc                                                                                                   >  5:60285/1‑148 (MQ=255)
                                                    gCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAg                                                                                                >  5:98751/1‑149 (MQ=255)
                                                     ccAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGc                                                                                               <  8:166623/149‑1 (MQ=255)
                                                     ccAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGc                                                                                               >  8:197830/1‑149 (MQ=255)
                                                         gCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGa                                                                                           <  5:89342/149‑1 (MQ=255)
                                                          cccAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCAGCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCGCGCCGCGTg                                                                                            >  8:109272/1‑147 (MQ=255)
                                                           ccAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAg                                                                                         <  6:98751/149‑1 (MQ=255)
                                                           ccAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAg                                                                                         <  3:135100/149‑1 (MQ=255)
                                                                atcaATGAAGGTTGGTCGACCTCGCCTGCGTTGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                                     <  8:176430/146‑1 (MQ=255)
                                                                agcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGTCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                                     <  6:72967/148‑1 (MQ=255)
                                                                    aggaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCGGAGcacc                                                                                <  4:214023/149‑1 (MQ=255)
                                                                    aggaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGACCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGcac                                                                                 >  1:166547/1‑148 (MQ=255)
                                                                          gtaggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGaa                                                                          >  5:62313/1‑149 (MQ=255)
                                                                             ggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATgg                                                                       <  5:203141/149‑1 (MQ=255)
                                                                              gtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATgg                                                                       <  7:64391/148‑1 (MQ=255)
                                                                                       gCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGt                                                             >  6:42332/1‑149 (MQ=255)
                                                                                        ccTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                                   <  4:26028/142‑1 (MQ=255)
                                                                                        ccTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                                   >  3:26028/1‑142 (MQ=255)
                                                                                        ccTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGtt                                                            <  7:135390/149‑1 (MQ=255)
                                                                                           ctgggcctggcCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGAc                                                         <  8:169257/143‑1 (MQ=255)
                                                                                            cGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGTGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGAc                                                         >  5:124593/1‑148 (MQ=255)
                                                                                             gTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCa                                                       <  7:197830/149‑1 (MQ=255)
                                                                                                      cAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCCCGGAATGGGAGCCGAAGTTGACCAGCGCGCCg                                               >  4:166836/1‑148 (MQ=255)
                                                                                                      cAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCg                                               >  5:163362/1‑148 (MQ=255)
                                                                                                            cTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcgc                                         <  6:163362/148‑1 (MQ=255)
                                                                                                              gCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCaa                                       >  2:93055/1‑148 (MQ=255)
                                                                                                               cACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATa                                     <  4:81227/149‑1 (MQ=255)
                                                                                                                aCGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAg                                    >  1:177388/1‑149 (MQ=255)
                                                                                                                      atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                                     >  1:42218/1‑110 (MQ=255)
                                                                                                                      atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                                     <  2:42218/110‑1 (MQ=255)
                                                                                                                      atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGag                              <  7:109272/149‑1 (MQ=255)
                                                                                                                      atcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGag                              <  8:186737/149‑1 (MQ=255)
                                                                                                                       tcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGag                              >  7:176784/1‑148 (MQ=255)
                                                                                                                        caGCGCGGTGGCCTGGCTGGGAAGGGCCATGCCCAGCAGCTTGGGTAGCTGGCCCAACGCAA‑TGACGATGGCTGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTAGACAAGCGCGCCGAAGCGCAAAAGACCGacca                            >  8:220704/1‑149 (MQ=255)
                                                                                                                         aGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGa                                                                                           >  4:208731/1‑85 (MQ=255)
                                                                                                                         aGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGa                                                                                           <  3:208731/85‑1 (MQ=255)
                                                                                                                          gcgcGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcag                           >  8:131662/1‑148 (MQ=255)
                                                                                                                                    ctggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                                     >  1:216763/1‑80 (MQ=255)
                                                                                                                                    ctggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                                     <  2:216763/80‑1 (MQ=255)
                                                                                                                                    ctggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGaaaa                >  7:82146/1‑149 (MQ=255)
                                                                                                                                      tgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATg              >  6:115521/2‑149 (MQ=255)
                                                                                                                                       gctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGa                               <  1:39136/131‑1 (MQ=255)
                                                                                                                                       gctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGa                               >  2:39136/1‑131 (MQ=255)
                                                                                                                                       gctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgc             <  8:176784/149‑1 (MQ=255)
                                                                                                                                           ggCAGGGCCATGCCCAGGAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgccgcc         <  8:82146/149‑1 (MQ=255)
                                                                                                                                           ggCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATTGACGATGGCGGCACCCAGCGTGAAGCCGAGCCCCAGGGAATGGGAGACGAAGTTGACCACCGCGCCGACGCGCAATAGCCCCAGCAGCTGATGGAAAATgccgc          >  1:56117/1‑149 (MQ=255)
                                                                                                                                            gCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCa                                                       >  8:84196/1‑102 (MQ=255)
                                                                                                                                            gCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCa                                                       <  7:84196/102‑1 (MQ=255)
                                                                                                                                            gCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgccgcc         >  6:261796/1‑148 (MQ=255)
                                                                                                                                                ggCCCTGCCCAGCCGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAgg     >  8:16172/1‑148 (MQ=255)
                                                                                                                                                  ccATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAg  <  7:131662/149‑1 (MQ=255)
                                                                                                                                                  ccATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAg  <  5:264968/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAA‑TGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAG  >  NC_002947/86097‑86391

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