Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,772,213 A→C 5.8% I269S (ATC→AGC)  betT‑I ← choline transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,772,2130AC 5.8% 251.6 / 5.8 86I269S (ATC→AGC) betT‑Icholine transporter
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base C (2/3);  total (42/44)
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.77e-01

GCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTTGC  >  NC_002947/5772065‑5772362
                                                                                                                                                     |                                                                                                                                                     
gCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgat                                                                                                                                                      >  1:2124108/1‑150 (MQ=255)
gCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgat                                                                                                                                                      <  2:213728/150‑1 (MQ=255)
gCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgat                                                                                                                                                      >  1:1634067/1‑150 (MQ=255)
   ggCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGtt                                                                                                                                                   >  2:1293161/1‑150 (MQ=255)
     cTTCGCCTTCATACACATAGGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAg                                                                                                                                                <  2:117471/150‑1 (MQ=255)
      ttCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAg                                                                                                                                                <  1:1293161/150‑1 (MQ=255)
       tCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgt                                                                                                                                               <  1:2172892/150‑2 (MQ=255)
         gacTTCAGACACTTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGtt                                                                                                                                             >  2:1744939/3‑150 (MQ=255)
                   acaTAGAGGTCGAAGGTTTTCAGCACGATGCAGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGAGGCAGGCCTTTATGCACACTAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcgg                                                                                                                                   >  1:756191/1‑150 (MQ=255)
                     aTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGgcggcg                                                                                                                                 >  1:1704210/1‑150 (MQ=255)
                      tAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcgttgg                                                                                                                                <  2:1896314/150‑5 (MQ=255)
                          ggTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGc                                                                                                                            <  2:1339150/150‑1 (MQ=255)
                            tCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAg                                                                                                                                       >  1:547755/1‑137 (MQ=255)
                            tCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAg                                                                                                                                       <  2:547755/137‑1 (MQ=255)
                            tCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATg                                                                                                                             >  1:943010/1‑147 (MQ=255)
                            tCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATg                                                                                                                             <  2:943010/147‑1 (MQ=255)
                              gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGgcggcg                                                                                                                                 <  1:467090/141‑1 (MQ=255)
                              gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGgcggcg                                                                                                                                 >  2:467090/1‑141 (MQ=255)
                                 ggTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTtctc                                                                                                                     <  1:925315/150‑1 (MQ=255)
                                   ttttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGACGGCGGATGCCGTTCTCCa                                                                                                                   >  1:239348/1‑150 (MQ=255)
                                   ttttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGCTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAg                                                                                                                                                <  2:2123988/121‑1 (MQ=255)
                                   ttttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGCTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAg                                                                                                                                                >  1:2123988/1‑121 (MQ=255)
                                        aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgt                                                                                                                                               <  2:574017/117‑2 (MQ=255)
                                        aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgt                                                                                                                                               >  1:574017/1‑116 (MQ=255)
                                         gCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCg                                                                                                             <  1:1152633/150‑1 (MQ=255)
                                          catcaTGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgatg                                                                                                                                                     <  1:136464/105‑1 (MQ=255)
                                          catcaTGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgatg                                                                                                                                                     >  2:136464/5‑109 (MQ=255)
                                          cACGATGCTGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGa                                                                                                            >  2:1810175/1‑150 (MQ=255)
                                            cGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATa                                                                                                          <  2:1704210/150‑1 (MQ=255)
                                               tGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCg                                                                                                       >  2:1534482/1‑150 (MQ=255)
                                                  cGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTg                                                                                                    <  1:1744939/150‑1 (MQ=255)
                                                   gTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGc                                                                                                                >  2:2220123/1‑137 (MQ=255)
                                                   gTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGc                                                                                                                <  1:2220123/137‑1 (MQ=255)
                                                   gTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGc                                                                                                   >  1:868235/1‑150 (MQ=255)
                                                   gTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGc                                                                                                   >  1:550089/1‑150 (MQ=255)
                                                        ggTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGACCAGCATGATGTTCAGGTCGGACCGGCGGCGGATGCCCTCCTCCACGCCCGATACCGCTGCCACgg                                                                                              <  2:1893989/150‑1 (MQ=255)
                                                        ggTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                              <  1:1129911/150‑1 (MQ=255)
                                                        ggTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                              >  2:1611455/1‑150 (MQ=255)
                                                         gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGt                                                                                             <  2:586217/150‑1 (MQ=255)
                                                          tAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                            <  1:1200517/150‑1 (MQ=255)
                                                           aGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAg                                                                                           <  2:179774/150‑1 (MQ=255)
                                                               gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCAc                                                                                       >  2:905535/1‑150 (MQ=255)
                                                               gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCAc                                                                                       <  2:868235/150‑1 (MQ=255)
                                                                ccAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACg                                                                                      >  1:1298126/1‑150 (MQ=255)
                                                                 cAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGc                                                                                     <  2:2124108/150‑1 (MQ=255)
                                                                      ttCTGCACGAAGCCGTTGATCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGAc                                                                                                                                         <  1:155323/93‑1 (MQ=255)
                                                                      ttCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGAc                                                                                                                                         >  2:155323/1‑93 (MQ=255)
                                                                            aCGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAcc                                                                          >  2:1600343/1‑150 (MQ=255)
                                                                             cGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCATGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATCCCGCCGCCACGGTAGCCACCCCCGCCATGACCa                                                                         <  2:1993433/150‑1 (MQ=255)
                                                                              gAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCAACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAg                                                                        >  1:1628326/1‑150 (MQ=255)
                                                                               aaGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCACGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAgg                                                                       >  1:725382/1‑150 (MQ=255)
                                                                                  ccGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATc                                                                    <  1:1534482/150‑1 (MQ=255)
                                                                                   cGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAgg                                                                                                                                               >  2:1811941/1‑74 (MQ=255)
                                                                                   cGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAgg                                                                                                                                               <  1:1811941/74‑1 (MQ=255)
                                                                                     ttGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcacc                                                                 <  1:1417477/150‑1 (MQ=255)
                                                                                        aGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCaccacc                                                              >  2:561667/1‑150 (MQ=255)
                                                                                         gCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCAccacca                                                             <  1:1810175/150‑1 (MQ=255)
                                                                                                gCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGt                                                       >  2:1754464/1‑150 (MQ=255)
                                                                                                   gggTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATc                                                                    >  2:1865315/1‑134 (MQ=255)
                                                                                                   gggTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATc                                                                    <  1:1865315/134‑1 (MQ=255)
                                                                                                    ggTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGt                                                   >  2:1020597/1‑150 (MQ=255)
                                                                                                          gccgccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCt                                             >  1:954091/1‑150 (MQ=255)
                                                                                                              cccAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTgg                                         <  2:725382/150‑1 (MQ=255)
                                                                                                               ccAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGt                                        <  2:954091/150‑1 (MQ=255)
                                                                                                                cagcagCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTc                                       <  2:1634067/150‑1 (MQ=255)
                                                                                                                 agcagcACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTcc                                      <  1:1263121/150‑1 (MQ=255)
                                                                                                                      caCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGcc                                 >  1:1962679/1‑150 (MQ=255)
                                                                                                                       aCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                           >  2:1005218/1‑107 (MQ=255)
                                                                                                                       aCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                           <  1:1005218/107‑1 (MQ=255)
                                                                                                                        cGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTtct                                                                                                                      <  2:755812/63‑1 (MQ=255)
                                                                                                                        cGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTtct                                                                                                                      >  1:755812/1‑63 (MQ=255)
                                                                                                                         gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTgg                                         >  2:282128/1‑139 (MQ=255)
                                                                                                                         gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTgg                                         <  1:282128/139‑1 (MQ=255)
                                                                                                                         gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGaa                              <  2:756191/150‑1 (MQ=255)
                                                                                                                         gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGCTCCATGCCGaa                              >  2:1117510/1‑150 (MQ=255)
                                                                                                                         gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGATGTCGTGCGTGCCGac                              >  2:2248509/1‑149 (MQ=255)
                                                                                                                          aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTCTGGTCCATGCCGAAc                             >  2:1631299/1‑150 (MQ=255)
                                                                                                                              tCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATa                                                                                                          >  1:1311822/1‑69 (MQ=255)
                                                                                                                              tCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATa                                                                                                          <  2:1311822/69‑1 (MQ=255)
                                                                                                                                aaCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTAt                       >  1:205806/1‑150 (MQ=255)
                                                                                                                                 aCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTAtt                      <  1:1754464/150‑1 (MQ=255)
                                                                                                                                   aGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTcc                    >  2:1915650/1‑150 (MQ=255)
                                                                                                                                    ggCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCa                   <  1:80016/150‑1 (MQ=255)
                                                                                                                                     gCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGt                                                                                             <  2:1209914/75‑1 (MQ=255)
                                                                                                                                     gCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGt                                                                                             >  1:1209914/1‑75 (MQ=255)
                                                                                                                                     gCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAg                  <  2:205806/150‑1 (MQ=255)
                                                                                                                                        aCTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTcgagatc               >  1:1662890/1‑144 (MQ=255)
                                                                                                                                          tGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCgg             >  1:238862/1‑150 (MQ=255)
                                                                                                                                          tGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCgg             >  1:1350050/1‑150 (MQ=255)
                                                                                                                                           gAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGa            <  1:1611455/150‑1 (MQ=255)
                                                                                                                                            aaCAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGAc           <  1:905535/150‑1 (MQ=255)
                                                                                                                                                   tgatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt    <  1:309169/150‑1 (MQ=255)
                                                                                                                                                     atgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTTGc  >  2:797127/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTTGC  >  NC_002947/5772065‑5772362

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