Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 354,876 A→C 6.8% F100C (TTC→TGC)  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,8760AC 6.8% 222.5 / 5.6 82F100C (TTC→TGC) cbcWcholine/betaine/carnitine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (42/34);  new base C (3/3);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGGGTCAGGGTAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGG  >  NC_002947/354729‑355025
                                                                                                                                                   |                                                                                                                                                     
aaCAGGGTCAGGGTAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAg                                                                                                                                                     <  1:1405473/150‑1 (MQ=255)
aaCAGGGTCAGGGTAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAg                                                                                                                                                     <  2:1807307/150‑1 (MQ=255)
      gtcagggtAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGAcc                                                                                                                                               >  1:224755/1‑150 (MQ=255)
      gtcagggtAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGAcc                                                                                                                                               >  1:1453171/1‑150 (MQ=255)
       gcagcgtAGGGATCATGTAGACCAAGGTGTGCACCGGCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGGGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGCAGGTCACCt                                                                                                                                              <  2:991505/149‑1 (MQ=255)
               gggATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                      <  1:1462105/150‑1 (MQ=255)
               gggATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAgg                                                                                                                                      <  2:224755/150‑1 (MQ=255)
                 gATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGaa                                                                                                                                                      >  2:1631315/1‑132 (MQ=255)
                 gATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGaa                                                                                                                                                      <  1:1631315/132‑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                                                                                                                  <  1:929661/150‑1 (MQ=255)
                                   gggCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCa                                                                                                                  >  2:429411/1‑150 (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:123690/1‑150 (MQ=255)
                                                 tCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGa                                                                                                    <  1:2391707/150‑1 (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                                                                                  <  1:1994830/136‑1 (MQ=255)
                                                                                 aaCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACa                                                                                  >  2:1994830/1‑136 (MQ=255)
                                                                                 aaCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                    <  2:423053/150‑1 (MQ=255)
                                                                                 aaCATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCg                                                                    <  2:1453171/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                                                           >  2:476335/1‑131 (MQ=255)
                                                                                                             ccAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGcc                                                           <  1:476335/131‑1 (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                                                      >  2:1058223/1‑124 (MQ=255)
                                                                                                                         cGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACtt                                                      <  1:1058223/124‑1 (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                                            >  2:169033/1‑117 (MQ=255)
                                                                                                                                          accGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATg                                            <  1:169033/117‑1 (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)
                                                                                                                                                   |                                                                                                                                                     
AACAGGGTCAGGGTAGGGATCAGGTAGACGAAGGTGGGCACCGTCTGCATCAGGTCGAGCACCGGGCGCATGGCGGTATAGAACATCGGCTTGTGCGCGGCGAGGATGCCCAGTGGCACGCCGATCACCACACATACCACCGTGGCGAAGGTGACCTGCGCCAGGGTCTCCATGGTTTCCTGCCAGTAACCCAGGTTGAGGATCAGCAGGAACGACAGGGCGACAAATGCGGTCAGTGCCCACTTGCGCTGGATGAAATGCGCCAGGGCGGCGAATAGGGCGATAAGAACGAACGGG  >  NC_002947/354729‑355025

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: