Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 354,886 G→T 5.4% Q97K (CAG→AAG)  cbcW ← choline/betaine/carnitine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947354,8860GT 5.4% 285.3 / 2.0 89Q97K (CAG→AAG) cbcWcholine/betaine/carnitine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (46/38);  new base T (3/2);  total (49/40)
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.95e-01

GGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAA  >  NC_002947/354744‑355034
                                                                                                                                              |                                                                                                                                                    
gggATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                               <  2:224755/150‑1 (MQ=255)
gggATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                               <  1:1462105/150‑1 (MQ=255)
  gATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGCAGGTCACCTTCGCCa                                                                                                                                                 <  2:819215/146‑1 (MQ=255)
  gATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGCAGGTCACCTTCGCCa                                                                                                                                                 >  1:819215/1‑146 (MQ=255)
   aTCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGtc                                                                                                                                            >  2:993085/1‑150 (MQ=255)
   aTCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCCTGGCGCAGGTCACCTTCGCCACGGtg                                                                                                                                            <  2:226281/150‑2 (MQ=255)
       ggTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCa                                                                                                                                        <  1:1400567/150‑1 (MQ=255)
        gTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCAt                                                                                                                                       <  1:993085/150‑1 (MQ=255)
           gACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATGCTACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGt                                                                                                                                    >  1:2202105/1‑150 (MQ=255)
              gAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTc                                                                                                                                 >  2:1364779/1‑150 (MQ=255)
                    gggCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCa                                                                                                                           >  2:429411/1‑150 (MQ=255)
                    gggCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCa                                                                                                                           <  1:929661/150‑1 (MQ=255)
                       cACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTa                                                                                                                        >  1:1713135/1‑150 (MQ=255)
                             cTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAg                                                                                                                  <  1:264603/150‑1 (MQ=255)
                              tGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAgg                                                                                                                 <  1:177044/150‑1 (MQ=255)
                                  tCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGa                                                                                                             >  1:521604/1‑150 (MQ=255)
                                  tCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGa                                                                                                             <  1:2391707/150‑1 (MQ=255)
                                  tCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGa                                                                                                             >  1:123690/1‑150 (MQ=255)
                                     ggTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                               >  1:19720/1‑113 (MQ=255)
                                     ggTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                               <  2:19720/113‑1 (MQ=255)
                                      gTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGAt                                                                                                         >  2:1957827/1‑150 (MQ=255)
                                         gAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGtctc                                                                                                                                          >  2:415737/1‑114 (MQ=255)
                                         gAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGtctc                                                                                                                                          <  1:415737/114‑1 (MQ=255)
                                            cACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATcagcag                                                                                                   >  1:2271733/1‑150 (MQ=255)
                                                           ggTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGAc                                                                                    >  2:866134/1‑150 (MQ=255)
                                                           ggTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGAc                                                                                    >  2:685078/1‑150 (MQ=255)
                                                             tataGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACaa                                                                                  >  2:928822/1‑150 (MQ=255)
                                                                aGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATg                                                                               <  1:685078/150‑1 (MQ=255)
                                                                  aaCCTCGGCTTGTGCGCGGCGAGGAGGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCTAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                             <  2:810494/150‑1 (MQ=255)
                                                                  aaCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACa                                                                                           >  2:1994830/1‑136 (MQ=255)
                                                                  aaCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACa                                                                                           <  1:1994830/136‑1 (MQ=255)
                                                                  aaCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                             <  2:1453171/150‑1 (MQ=255)
                                                                  aaCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                             <  2:423053/150‑1 (MQ=255)
                                                                    cATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGt                                                                           >  2:1906844/1‑150 (MQ=255)
                                                                     aTCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAgg                                                                                         <  1:2191929/135‑1 (MQ=255)
                                                                     aTCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAgg                                                                                         >  2:2191929/1‑135 (MQ=255)
                                                                      tCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCa                                                                         <  1:429411/150‑1 (MQ=255)
                                                                         gCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTg                                                                      >  1:228313/1‑150 (MQ=255)
                                                                         gCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTg                                                                      <  1:928822/150‑1 (MQ=255)
                                                                          cTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGc                                                                     <  1:1957827/150‑1 (MQ=255)
                                                                            tgtgCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGccc                                                                   >  1:2178959/1‑150 (MQ=255)
                                                                                              ccAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                             <  2:610633/122‑1 (MQ=255)
                                                                                              ccAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                             >  1:610633/1‑122 (MQ=255)
                                                                                              ccAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGcc                                                                    <  1:476335/131‑1 (MQ=255)
                                                                                              ccAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGcc                                                                    >  2:476335/1‑131 (MQ=255)
                                                                                               cAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCa                                                                                                                   <  1:168403/83‑1 (MQ=255)
                                                                                               cAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCa                                                                                                                   >  2:168403/1‑83 (MQ=255)
                                                                                                   ggCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGcc                                            <  1:866134/150‑1 (MQ=255)
                                                                                                      aCGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGCACGGATGACCAACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAgg                                         >  2:1269980/1‑150 (MQ=255)
                                                                                                      aCGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGt                                                                       <  2:648120/120‑1 (MQ=255)
                                                                                                      aCGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGt                                                                       >  1:648120/1‑120 (MQ=255)
                                                                                                      aCGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAgg                                         >  2:2343889/1‑150 (MQ=255)
                                                                                                        gCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGTGCCAAGGTCTCTATCGCTTCCTTCCAGTAACCCAGCCTGAGGATCAGCACGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTCGCGCTGGATGAACTGCGCCAGggc                                       >  2:428373/1‑150 (MQ=255)
                                                                                                        gCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGggc                                       >  1:2004647/1‑150 (MQ=255)
                                                                                                          cGATCACCACACATACCACCGTGGGGAATGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGggcgg                                     >  1:1531789/1‑150 (MQ=255)
                                                                                                          cGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACtt                                                               <  1:1058223/124‑1 (MQ=255)
                                                                                                          cGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACtt                                                               >  2:1058223/1‑124 (MQ=255)
                                                                                                                ccacACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAGTGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAggg                                        >  2:758751/1‑141 (MQ=255)
                                                                                                                ccacACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGACTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAGTGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAggg                                        <  1:758751/141‑1 (MQ=255)
                                                                                                                  acacaTACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAgg                             >  1:1334977/1‑150 (MQ=255)
                                                                                                                   cacaTCCCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAggg                            <  1:511255/150‑1 (MQ=255)
                                                                                                                      aTACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGaaa                                                  <  2:1292573/125‑1 (MQ=255)
                                                                                                                      aTACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGaaa                                                  >  1:1292573/1‑125 (MQ=255)
                                                                                                                      aTACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGATGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGa                         >  2:189232/1‑150 (MQ=255)
                                                                                                                          caccGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                             >  1:1158080/1‑94 (MQ=255)
                                                                                                                          caccGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                             <  2:1158080/94‑1 (MQ=255)
                                                                                                                          caccGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAg                     >  1:2012931/1‑150 (MQ=255)
                                                                                                                           accGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATg                                                     <  1:169033/117‑1 (MQ=255)
                                                                                                                           accGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATg                                                     >  2:169033/1‑117 (MQ=255)
                                                                                                                            ccGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATc                                                                                                        <  2:920778/65‑1 (MQ=255)
                                                                                                                            ccGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATc                                                                                                        >  1:920778/1‑65 (MQ=255)
                                                                                                                             cctGGCGCAGGTCACCTTCGCCAGGGTCTCCATGGTTTCCTGCGAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAgaac                  >  1:1162391/3‑150 (MQ=255)
                                                                                                                             cctGGCGCAGGTCACCTTCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACACAATGGGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGCCGAATAGGGCGATAAgaac                  >  2:2393580/3‑150 (MQ=255)
                                                                                                                                  cGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACgg            >  1:319846/1‑150 (MQ=255)
                                                                                                                                   gAAGGTTACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACgg            <  2:2178959/150‑1 (MQ=255)
                                                                                                                                   gAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACgg            >  1:1583981/1‑150 (MQ=255)
                                                                                                                                   gAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACgg            >  1:1226099/1‑150 (MQ=255)
                                                                                                                                    aaGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACggg           <  2:1713135/150‑1 (MQ=255)
                                                                                                                                      ggTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGtt         >  1:2029541/1‑150 (MQ=255)
                                                                                                                                      ggTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACACAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGtt         >  1:2103642/1‑150 (MQ=255)
                                                                                                                                       gTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAggg                                        <  1:2114169/118‑1 (MQ=255)
                                                                                                                                       gTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAggg                                        >  2:2114169/1‑118 (MQ=255)
                                                                                                                                         gACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGt                                                                       <  1:878591/85‑1 (MQ=255)
                                                                                                                                         gACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGt                                                                       >  2:878591/1‑85 (MQ=255)
                                                                                                                                         gACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACTGGTTgaa      <  2:123690/150‑1 (MQ=255)
                                                                                                                                           ccTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTgaaga    <  2:1531789/150‑1 (MQ=255)
                                                                                                                                            cTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTgaagaa   >  1:1704073/1‑150 (MQ=255)
                                                                                                                                             tGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGaaa  <  2:2197164/150‑1 (MQ=255)
                                                                                                                                             tGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGaaa  >  1:164638/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGGTTGAAGAAA  >  NC_002947/354744‑355034

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: