Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,493 T→A 6.5% N98I (AAT→ATT)  PP_0076 ← choline betaine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294787,4930TA 6.5% 236.0 / 5.2 77N98I (AAT→ATT) PP_0076choline betaine‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (34/38);  new base A (2/3);  total (36/41)
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.84e-01

CGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAG  >  NC_002947/87347‑87641
                                                                                                                                                  |                                                                                                                                                    
cGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCTCCCGGGCATTCGACACATAGTCGGGAACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGGGCGACAGCTGCTGCACATCGCCATTg                                                                                                                                                    <  4:306278/149‑1 (MQ=255)
      cGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                               <  3:139256/148‑1 (MQ=255)
      cGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTAGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGTGAt                                                                                                                                               <  1:117036/148‑1 (MQ=255)
       gTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATg                                                                                                                                              >  1:316201/1‑148 (MQ=255)
        tACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGaa                                                                                                                                            <  5:115094/149‑1 (MQ=255)
        tACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGATgtg                                                                                                                                            <  4:23961/149‑3 (MQ=255)
            gTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                               <  4:136701/142‑1 (MQ=255)
            gTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGAt                                                                                                                                               >  3:136701/1‑142 (MQ=255)
            gTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTg                                                                                                                                        >  3:225878/1‑149 (MQ=255)
             ttttttGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGt                                                                                                                                       >  1:117515/1‑149 (MQ=255)
                   gTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATgg                                                                                                                                  <  2:117515/148‑1 (MQ=255)
                      gAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAt                                                                                                                                                      >  3:73214/1‑125 (MQ=255)
                      gAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAt                                                                                                                                                      <  4:73214/125‑1 (MQ=255)
                      gAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACtt                                                                                                                                         <  6:3646/138‑1 (MQ=255)
                      gAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACtt                                                                                                                                         >  5:3646/1‑138 (MQ=255)
                                gCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTgcgc                                                                                                                    >  6:72242/1‑149 (MQ=255)
                                        ttGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATc                                                                                                             <  4:225878/148‑1 (MQ=255)
                                            aGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGtt                                                                                                        >  7:22771/1‑149 (MQ=255)
                                               tCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTcc                                                                                                     <  8:22771/149‑1 (MQ=255)
                                                 aGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTcc                                                                                                     >  3:112094/1‑147 (MQ=255)
                                                     aaGTCCTTCACCCCGGCATTCCACACTTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCTTTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGaa                                                                                               <  2:7442/149‑1 (MQ=255)
                                                     aaGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGaa                                                                                               <  5:72242/149‑1 (MQ=255)
                                                      aGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                              <  2:316201/149‑1 (MQ=255)
                                                      aGTCCTTCACCCCGGCATTCCACACATAGTCGGGCAATGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                              <  8:321191/149‑1 (MQ=255)
                                                          cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGaa                                                                                                                                <  2:249308/111‑1 (MQ=255)
                                                          cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGaa                                                                                                                                >  1:249308/1‑111 (MQ=255)
                                                          cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGt                                                                                          >  1:201418/1‑149 (MQ=255)
                                                             cACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTc                                                                                         <  2:201418/147‑1 (MQ=255)
                                                             cACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCg                                                                                        <  4:112094/148‑1 (MQ=255)
                                                                cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACtt                                                                                                                                         <  6:37592/96‑1 (MQ=255)
                                                                cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCt                                                                                    >  2:255087/1‑149 (MQ=255)
                                                                cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCt                                                                                    >  6:318172/1‑149 (MQ=255)
                                                                 ccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGaa                                                                                               >  4:823/1‑137 (MQ=255)
                                                                 ccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGaa                                                                                               <  3:823/137‑1 (MQ=255)
                                                                 ccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                   <  7:331995/149‑1 (MQ=255)
                                                                  cGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTgg                                                                                  >  8:99074/1‑149 (MQ=255)
                                                                    gCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGcc                                                                                <  1:255087/149‑1 (MQ=255)
                                                                           acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTc                                                                                         <  1:147474/133‑1 (MQ=255)
                                                                           acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTc                                                                                         >  2:147474/1‑133 (MQ=255)
                                                                              caTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTgg                                                                                  <  5:252045/137‑1 (MQ=255)
                                                                              caTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTgg                                                                                  >  6:252045/1‑137 (MQ=255)
                                                                               aTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCg                                                                               >  6:126310/1‑139 (MQ=255)
                                                                               aTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCg                                                                               <  5:126310/139‑1 (MQ=255)
                                                                               aTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGgc                                                                     >  6:214391/1‑149 (MQ=255)
                                                                               aTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACACCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTg                                                                         >  4:106523/1‑145 (MQ=255)
                                                                               aTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACACCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTg                                                                         <  3:106523/145‑1 (MQ=255)
                                                                                     gggCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTtg                                                                                                                          >  7:1260/1‑90 (MQ=255)
                                                                                     gggCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTtg                                                                                                                          <  8:1260/90‑1 (MQ=255)
                                                                                      ggCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTa                                                              <  5:318172/149‑1 (MQ=255)
                                                                                              gAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGt                                                               >  7:104251/1‑140 (MQ=255)
                                                                                              gAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGt                                                               <  8:104251/140‑1 (MQ=255)
                                                                                                 ggTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTa                                                   <  5:214391/149‑1 (MQ=255)
                                                                                                    gAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGACGCCGTAGATGATCGGTACgg                                                >  5:269951/1‑149 (MQ=255)
                                                                                                           gTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCt                                         >  2:116459/1‑149 (MQ=255)
                                                                                                                cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCgg                                                     <  2:214415/132‑1 (MQ=255)
                                                                                                                cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCgg                                                     >  1:214415/1‑132 (MQ=255)
                                                                                                                 tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGAc                                                                                      >  2:114214/1‑98 (MQ=255)
                                                                                                                 tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGAc                                                                                      <  1:114214/98‑1 (MQ=255)
                                                                                                                 tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAg                                                                       <  8:265490/113‑1 (MQ=255)
                                                                                                                 tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAg                                                                       >  7:265490/1‑113 (MQ=255)
                                                                                                                 tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCa                                             >  4:218758/1‑139 (MQ=255)
                                                                                                                 tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCa                                             <  3:218758/139‑1 (MQ=255)
                                                                                                                 tCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGa                                    >  6:95842/1‑148 (MQ=255)
                                                                                                                       gtGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCAc                             <  5:153477/148‑1 (MQ=255)
                                                                                                                          cgcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTg                          <  1:208349/149‑1 (MQ=255)
                                                                                                                           gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTgg                         <  2:251939/149‑1 (MQ=255)
                                                                                                                           gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTgg                         >  6:62375/1‑149 (MQ=255)
                                                                                                                           gcgACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTgg                         >  1:251939/1‑149 (MQ=255)
                                                                                                                             gACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGt                        >  4:336912/1‑148 (MQ=255)
                                                                                                                             gACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGt                        >  4:342204/1‑148 (MQ=255)
                                                                                                                                       gttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCa                  <  4:202800/141‑1 (MQ=255)
                                                                                                                                       gttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAgagat             >  3:202800/4‑145 (MQ=255)
                                                                                                                                       gcACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCtt             >  6:69529/1‑149 (MQ=255)
                                                                                                                                        ttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATc                                  <  4:59700/125‑1 (MQ=255)
                                                                                                                                        ttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATc                                  >  3:59700/3‑127 (MQ=255)
                                                                                                                                         acaTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTc           <  5:62375/149‑1 (MQ=255)
                                                                                                                                           aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTc           >  1:343592/1‑147 (MQ=255)
                                                                                                                                           aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCa         <  1:116459/149‑1 (MQ=255)
                                                                                                                                                 aTTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACa   <  7:99074/149‑1 (MQ=255)
                                                                                                                                                  ttGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAg  <  2:343592/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAACAG  >  NC_002947/87347‑87641

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