Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,772,157 A→C 6.4% 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,1570AC 6.4% 191.4 / 4.4 63I269S (ATC→AGC) betT‑Icholine transporter
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base C (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTT  >  NC_002947/5772011‑5772304
                                                                                                                                                  |                                                                                                                                                   
cGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgatg                                                                                                                                                   >  8:142450/1‑149 (MQ=255)
  gCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGtt                                                                                                                                                 >  2:165906/1‑149 (MQ=255)
       gCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAgg                                                                                                                                             <  2:26537/148‑1 (MQ=255)
        ccTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTc                                                                                                                                                <  4:226614/144‑1 (MQ=255)
        ccTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTc                                                                                                                                                >  3:226614/1‑144 (MQ=255)
            cATCCACATCGAGGTCGACGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCTTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGa                                                                                                                                        <  6:130130/148‑1 (MQ=255)
              tACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGAGGTTCAGGTTGGACAg                                                                                                                                     >  4:174887/1‑149 (MQ=255)
               acacATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAgg                                                                                                                                    >  8:11773/1‑149 (MQ=255)
                         gTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGc                                                                                                                          >  7:131240/1‑149 (MQ=255)
                            gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                       >  3:286701/1‑149 (MQ=255)
                            gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                       >  2:72696/1‑149 (MQ=255)
                            gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGACGt                                                                                                                       >  6:67095/1‑149 (MQ=255)
                            gAAGGTTTTCAGCACGATGCCGTTCAAGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGt                                                                                                                       >  1:9676/1‑149 (MQ=255)
                                gTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGa                                                                                                                                        <  6:46170/128‑1 (MQ=255)
                                gTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGa                                                                                                                                        >  5:46170/1‑128 (MQ=255)
                                   ttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCAc                                                                                                                >  5:247633/1‑149 (MQ=255)
                                   ttCAGCACGATGCCGTTCAGGTAGTCGCCAACCTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCa                                                                                                                 <  1:122650/148‑1 (MQ=255)
                                    tCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCAc                                                                                                                <  5:128494/148‑1 (MQ=255)
                                      aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGcc                                                                                                             >  3:161002/1‑149 (MQ=255)
                                      aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGcc                                                                                                             <  3:292795/149‑1 (MQ=255)
                                       gCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGccc                                                                                                            <  7:11773/149‑1 (MQ=255)
                                        cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                            <  2:41945/116‑3 (MQ=255)
                                        cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                            <  2:37214/116‑3 (MQ=255)
                                        cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                            >  1:41945/1‑114 (MQ=255)
                                        cACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                            >  1:37214/1‑114 (MQ=255)
                                           gATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATa                                                                                                        >  7:73840/1‑149 (MQ=255)
                                                     aGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCa                                                                                               <  7:142450/148‑1 (MQ=255)
                                                       gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                            >  5:144041/1‑149 (MQ=255)
                                                       gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                            <  4:286701/149‑1 (MQ=255)
                                                        tAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                            >  7:187408/1‑148 (MQ=255)
                                                        tAGTCGCCAAAGTTCTGCACGATGCCGTTGAGCAGGTGCAGGGTTTCGACGCCCTGCAGCACGTAGATAAACAGGCCACTGATCAGCATGATGTTCAGGTTGGACAGGCGTCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGt                                                                                           <  2:9676/149‑1 (MQ=255)
                                                           tCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGc                                                                                        >  5:188245/1‑149 (MQ=255)
                                                             gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCGGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCa                                                                                      >  2:325776/1‑149 (MQ=255)
                                                             gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCa                                                                                      >  8:210642/1‑149 (MQ=255)
                                                                   gTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                          <  2:24221/139‑1 (MQ=255)
                                                                   gTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCgg                                                                                >  3:235447/1‑149 (MQ=255)
                                                                   gTTCTGCACGAAGCCGTTGAGCAGGTGCACGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                          >  1:24221/1‑139 (MQ=255)
                                                                            gAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCa                                                                       <  3:174887/149‑1 (MQ=255)
                                                                                   ttGAGCAGGTGCAGGGTTTCGCCGCCCAGCATCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcac                                                                <  7:43227/149‑1 (MQ=255)
                                                                                   ttGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcac                                                                <  4:235447/149‑1 (MQ=255)
                                                                                   ttGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCAACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcac                                                                <  6:155386/149‑1 (MQ=255)
                                                                                       gCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcaccac                                                             >  4:58714/1‑148 (MQ=255)
                                                                                         aGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACcag                                                          >  8:204682/1‑149 (MQ=255)
                                                                                            tGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGACATGACCAGGATCACCACcagcag                                                       >  7:292768/1‑149 (MQ=255)
                                                                                               aGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGtt                                                    <  6:188245/149‑1 (MQ=255)
                                                                                                 ggTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTgg                                                  <  3:23244/149‑1 (MQ=255)
                                                                                                   tttCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTc                                                >  8:185466/1‑149 (MQ=255)
                                                                                                   tttCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTc                                                >  2:168698/1‑149 (MQ=255)
                                                                                                       gcggccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTg                                             <  5:67095/145‑1 (MQ=255)
                                                                                                       gccgccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTg                                             <  1:72696/148‑1 (MQ=255)
                                                                                                       gcagccCAGCAGGAGGAAGATCAACAGGCCACTGAACAGGATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATTCCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTTGTCTTg                                             <  8:292768/145‑1 (MQ=255)
                                                                                                        ccgccCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCt                                           <  6:247633/149‑1 (MQ=255)
                                                                                                              agcagcACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTc                                     >  1:163057/1‑149 (MQ=255)
                                                                                                              agcagcACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTc                                     <  3:58714/149‑1 (MQ=255)
                                                                                                                    aCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGcc                               <  8:73840/149‑1 (MQ=255)
                                                                                                                       aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGc                                <  6:304888/145‑1 (MQ=255)
                                                                                                                       aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGc                                >  5:304888/1‑145 (MQ=255)
                                                                                                                       aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGaa                            <  4:230637/149‑1 (MQ=255)
                                                                                                                       aaGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGaa                            <  1:168698/149‑1 (MQ=255)
                                                                                                                         gATCAACAGGCCACTGAACAGCATTATGTTCAGGTTGGTCAGGCGTCGTATGCGGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACa                          <  2:163057/149‑1 (MQ=255)
                                                                                                                          aTCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAg                         >  2:165608/1‑149 (MQ=255)
                                                                                                                             aaCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTa                      <  5:124772/149‑1 (MQ=255)
                                                                                                                                           cAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGAATTCCAGGCCGGACg        <  8:131240/149‑1 (MQ=255)
                                                                                                                                              cATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATa     >  5:138254/1‑149 (MQ=255)
                                                                                                                                                 gatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTg                                                   <  6:250052/100‑1 (MQ=255)
                                                                                                                                                 gatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTg                                                   >  5:250052/1‑100 (MQ=255)
                                                                                                                                                 gatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt  >  5:249223/1‑149 (MQ=255)
                                                                                                                                                 gatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt  >  4:18551/1‑149 (MQ=255)
                                                                                                                                                  atgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt  <  7:210642/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTT  >  NC_002947/5772011‑5772304

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