Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 284,586 T→A 6.2% intergenic (+33/+30) betT‑II → / ← atsK choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947284,5860TA 6.2% 210.8 / 2.7 75intergenic (+33/+30)betT‑II/atsKcholine/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base A (3/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

ACCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCT  >  NC_002947/284438‑284734
                                                                                                                                                    |                                                                                                                                                            
aCCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACGTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCt                                                                                                                                                              <  5:144245/149‑1 (MQ=255)
 ccTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCtt                                                                                                                                                             >  2:242068/1‑149 (MQ=255)
    gAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGcc                                                                                                                                                          >  8:49568/1‑149 (MQ=255)
     aGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGc                                                                                                                                                           <  6:329023/147‑1 (MQ=255)
      gCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGccc                                                                                                                                                         <  8:325877/148‑1 (MQ=255)
        gAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                                       >  7:18783/1‑148 (MQ=255)
          aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCTTCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                                    <  8:69793/149‑1 (MQ=255)
              aGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCACGGGCTTGCCCGCGAAGAgg                                                                                                                                               <  2:153619/149‑1 (MQ=255)
               gCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                              >  5:154584/1‑149 (MQ=255)
                ccAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                                     >  1:172717/1‑142 (MQ=255)
                ccAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                                     <  2:172717/142‑1 (MQ=255)
                 cAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGcc                                                                                                                                             >  8:57494/1‑148 (MQ=255)
                 cAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                                            <  6:368093/149‑1 (MQ=255)
                  aGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                                    >  4:176030/1‑141 (MQ=255)
                  aGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                                    <  3:176030/141‑1 (MQ=255)
                         tGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGc                                                                                                                                    <  3:117721/148‑1 (MQ=255)
                           aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAAc                                                                                                                                 <  2:120274/149‑1 (MQ=255)
                           aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                                    >  7:88808/1‑131 (MQ=255)
                           aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                                    <  8:88808/131‑1 (MQ=255)
                                gtggGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTa                                                                                                                            <  3:259304/149‑1 (MQ=255)
                                    gCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGCCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGaaa                                                                                                                        <  7:85542/149‑1 (MQ=255)
                                     cTACAGCAAGGAGCAGATCATCAACGCCCTCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAg                                                                                                                       <  7:168980/149‑1 (MQ=255)
                                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGCAGAGGCCAGACAGGCAACCGCTAGAAAgg                                                                                                                      >  7:325880/1‑149 (MQ=255)
                                       aCAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                                     >  7:172948/1‑149 (MQ=255)
                                         aGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                                            >  2:346507/1‑125 (MQ=255)
                                         aGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                                            <  1:346507/125‑1 (MQ=255)
                                           cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCAt                                                                                                                 >  2:285784/1‑149 (MQ=255)
                                             aGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                                <  4:123827/117‑1 (MQ=255)
                                             aGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                                >  3:123827/1‑117 (MQ=255)
                                              ggAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAATGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATcc                                                                                                              <  1:242068/149‑1 (MQ=255)
                                                 gCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCaa                                                                                                          <  8:172948/149‑1 (MQ=255)
                                                 gCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCaa                                                                                                          <  5:349967/149‑1 (MQ=255)
                                                    gATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCACATAGCCGGCGTTCCCGGCAATGGGCCAGAAACGCAACCGGGATCAAGGTGCCTACCCAAgg                                                                                                        >  2:70182/1‑148 (MQ=255)
                                                        atcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGt                                                                                                  >  7:296906/1‑149 (MQ=255)
                                                                     tCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCAt                                                                                      >  6:167989/1‑148 (MQ=255)
                                                                     tCGACAAGTAAGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCAGTACTGGCCTCTTAGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCACCCGCAAGAAAGGTGAATCCCCATGGCTGGTCGCCCGGTTCAt                                                                                      >  6:375628/1‑148 (MQ=255)
                                                                       gACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCa                                                                                    <  1:70182/148‑1 (MQ=255)
                                                                          cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTg                                                                                         >  6:358966/1‑140 (MQ=255)
                                                                          cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTg                                                                                         <  5:358966/140‑1 (MQ=255)
                                                                             tACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgccg                                                                             >  2:255881/1‑149 (MQ=255)
                                                                                       acaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTGTTCGCGGGCTTGCCCGCGATGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                                   <  3:270981/149‑1 (MQ=255)
                                                                                         aTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTcggtcggt                                                                  >  4:221039/1‑148 (MQ=255)
                                                                                         aTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTa                                                                 <  8:317297/149‑1 (MQ=255)
                                                                                            cAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAat                                                               >  8:333301/1‑148 (MQ=255)
                                                                                                  cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAc                                                        <  7:333301/149‑1 (MQ=255)
                                                                                                   tGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGAGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAcc                                                       >  6:367038/1‑149 (MQ=255)
                                                                                                     cACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCg                                                      <  1:285784/149‑1 (MQ=255)
                                                                                                     cACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                                      <  4:21693/132‑1 (MQ=255)
                                                                                                     cACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                                      >  3:21693/1‑132 (MQ=255)
                                                                                                        cTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑Caa                                               >  3:128698/1‑148 (MQ=255)
                                                                                                        cTGGGCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatc                                                              <  7:140778/137‑1 (MQ=255)
                                                                                                           gTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAATGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATg                                           >  1:118289/1‑149 (MQ=255)
                                                                                                              cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTg                                        <  6:154584/149‑1 (MQ=255)
                                                                                                               gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑Caaa                                              <  6:305477/142‑1 (MQ=255)
                                                                                                               gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑Caaa                                              >  5:305477/1‑142 (MQ=255)
                                                                                                               gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTgg                                       <  7:57494/149‑1 (MQ=255)
                                                                                                               gCTGACCAAGCTCTCCTGTCCGGGCCTGTTGGCGGGCTTGCCCTCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTgg                                       <  8:18783/149‑1 (MQ=255)
                                                                                                                     cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                                   <  1:45460/120‑1 (MQ=255)
                                                                                                                     cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                                   >  2:45460/1‑119 (MQ=255)
                                                                                                                     cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCAcg                                 >  3:104066/1‑149 (MQ=255)
                                                                                                                      aaGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCAcgc                                <  4:104066/149‑1 (MQ=255)
                                                                                                                              ctttcctgGCCTCTGCCCGGCCTTGCCCGCGAATAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAATGATGGTCGCCCGGTGCATCACCCGCCGGTTTGGTCGTTAATCATCCACCG‑‑‑CAAAATGCTGTGTCACGCGGTTGTcc                        <  8:232646/146‑1 (MQ=255)
                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTg                                        >  2:132827/1‑133 (MQ=255)
                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTg                                        <  1:132827/133‑1 (MQ=255)
                                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTcc                        <  1:255881/149‑1 (MQ=255)
                                                                                                                                gtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTccc                       >  1:256102/1‑149 (MQ=255)
                                                                                                                                gtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTgg                                       >  8:347269/1‑132 (MQ=255)
                                                                                                                                gtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTgg                                       <  7:347269/132‑1 (MQ=255)
                                                                                                                                    tgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAAc                  >  3:326498/1‑149 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGc                >  2:321835/1‑149 (MQ=255)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa              >  1:233100/1‑149 (MQ=255)
                                                                                                                                        ctctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa              >  2:239989/1‑149 (MQ=255)
                                                                                                                                             cgcgGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc          <  3:221039/148‑1 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCa         <  7:49568/148‑1 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt       >  4:240889/1‑149 (MQ=255)
                                                                                                                                               cgGGCGTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt       <  8:296906/149‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGGGAAGAGGCCAGACAGGCAACCGCTTGAAAGTTGCATCCCCAAGCATGGTCGCCCGGTGCATAACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGTTTGTCCCAGAACGCCACGTCATtc      <  5:375628/149‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCATGGGGTCGAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt       >  1:308151/1‑148 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTcc                        >  8:264447/1‑127 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTcc                        <  7:264447/127‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCt  >  7:322334/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                            
ACCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCG‑‑‑CAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCT  >  NC_002947/284438‑284734

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