Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,126,947 T→A 6.8% intergenic (+38/+29) PP_5374 → / ← PP_5375 choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,126,9470TA 6.8% 199.8 / 3.0 68intergenic (+38/+29)PP_5374/PP_5375choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (27/36);  new base A (3/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACC  >  NC_002947/6126805‑6127093
                                                                                                                                              |                                                                                                                                                  
cgcCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACTGCCATCAGcgg                                                                                                                                              >  8:168762/1‑147 (MQ=255)
 gcCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGt                                                                                                                                             <  1:77679/149‑1 (MQ=255)
       cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCt                                                                                                                                       >  3:70622/1‑149 (MQ=255)
             cTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGTAGAAAGCCAGCTAACCTAAAACCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGtt                                                                                                                                 >  1:104652/1‑149 (MQ=255)
             cTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAAACACAAGACCGACAGCTGATGGCTGTCGGTCTTgtgtc                                                                                                                                 <  7:303698/149‑2 (MQ=38)
                ggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGttt                                                                                                                                <  3:178576/147‑1 (MQ=255)
                ggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTg                                                                                                                              <  5:21986/149‑1 (MQ=255)
                ggAGTCACCCTCCAGGCATCCCGCCGGTGGCGGGGCCCTTGGCCGGCGAGGGGCTGGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAACGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGtttt                                                                                                                               <  4:184166/148‑1 (MQ=255)
                 gATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTgg                                                                                                                             <  7:304726/149‑1 (MQ=255)
                  aTTCAGCCGCAAGGCATCACGCCGGTGGCTGGGCGGTTGGCCTGCGTGGGTCTTGATATGGGGCTGACAGGGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCCTCAGGTGTCGGTCTTGTGTTTTggg                                                                                                                            <  4:188489/149‑1 (MQ=255)
                  aTTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTggg                                                                                                                            <  7:67947/149‑1 (MQ=255)
                      aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGTAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAgg                                                                                                                        >  8:169328/1‑149 (MQ=255)
                           cAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTAtt                                                                                                                   >  3:280856/1‑149 (MQ=255)
                             aGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTc                                                                                                                 <  6:38723/149‑1 (MQ=255)
                               gCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTg                                                                                                               >  7:102147/1‑149 (MQ=255)
                               gCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTg                                                                                                               <  4:148880/149‑1 (MQ=255)
                                cATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCCACAGACATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCtt                                                                                                               >  8:59226/1‑147 (MQ=255)
                                 aTCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGTAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCa                                                                                                             <  7:169328/149‑1 (MQ=255)
                                   cACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACAGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCAcc                                                                                                           >  5:153892/1‑149 (MQ=255)
                                      gCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGaa                                                                                                        <  6:154875/149‑1 (MQ=255)
                                        cGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAgg                                                                                                      <  5:80271/149‑1 (MQ=255)
                                        cGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAgg                                                                                                      >  6:64483/1‑149 (MQ=255)
                                          gTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCa                                                                                                    <  4:280856/149‑1 (MQ=255)
                                                 ggCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCGGGCTATTTCTGCACCGAAGGCAGGGCCAt                                                                                             <  1:47380/149‑1 (MQ=255)
                                                 ggCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGcc                                                                                               <  3:247121/147‑1 (MQ=255)
                                                      ttGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                         >  2:201539/1‑148 (MQ=255)
                                                        ggCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCa                                                                                       <  7:56411/148‑1 (MQ=255)
                                                              cGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCg                                                                                 <  6:26464/148‑1 (MQ=255)
                                                                aGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGaa                                                                              <  2:104652/149‑1 (MQ=255)
                                                                aGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCACGGCCATCAACCAGTCCCGGaa                                                                              <  7:72474/149‑1 (MQ=255)
                                                                  gggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACa                                                                            <  1:201539/149‑1 (MQ=255)
                                                                     cTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGcc                                                                         >  7:25021/1‑149 (MQ=255)
                                                                           atGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAATAACGACAGCCCAAAGGGGTCGGTATTGTGTTTTGGGGAGGCTATTTATGCCCCCAACGCAGGGCCCTCAAACACTCCCCGAACAGCCGccagc                                                                   >  6:31107/1‑145 (MQ=255)
                                                                             ggggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCt                                                                  <  8:102147/148‑1 (MQ=255)
                                                                               ggCTGACAGCTAAGAAAGCCATCTAACCTAAATCCACCCTGTAATACACAAGAACGACATCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTc                                                               <  2:82840/149‑1 (MQ=255)
                                                                                  tGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCCCGCTTTCACg                                                            >  7:57826/1‑149 (MQ=255)
                                                                                        cGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGAcc                                                      <  5:64483/149‑1 (MQ=255)
                                                                                         gaagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCa                                                     >  3:293607/1‑149 (MQ=255)
                                                                                            gaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGa                                                           <  4:291646/140‑1 (MQ=255)
                                                                                            gaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGa                                                           <  4:286963/140‑1 (MQ=255)
                                                                                            gaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGa                                                           >  3:291646/1‑140 (MQ=255)
                                                                                            gaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGa                                                           >  3:286963/1‑140 (MQ=255)
                                                                                              aaGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTc                                                >  1:137041/1‑149 (MQ=255)
                                                                                               aGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGGTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCg                                               >  7:123678/1‑149 (MQ=255)
                                                                                                tcctgCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCgg                                              <  1:293909/145‑1 (MQ=255)
                                                                                                 ccAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTTGGCAGGCTATTTCTGCACCGAAGGCAGGGCGATAAACCAGTCCCGGAACCGCCGCACGCTTTCACGAGGACCCAGCTCggg                                             >  4:261322/1‑148 (MQ=255)
                                                                                                          cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTc                                    >  3:17657/1‑149 (MQ=255)
                                                                                                            aaaGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCca                                  <  4:17657/149‑1 (MQ=255)
                                                                                                             aaGCCACCCTGGAAGACACAACACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCcag                                 <  7:223784/149‑1 (MQ=255)
                                                                                                                  aCCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggc                            <  4:70622/149‑1 (MQ=255)
                                                                                                                   cccTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcagg                             >  3:9958/1‑147 (MQ=255)
                                                                                                                       ggAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAAc                       >  4:135398/1‑149 (MQ=255)
                                                                                                                       ggAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAAc                       <  6:50101/149‑1 (MQ=255)
                                                                                                                       ccaaaacacAAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAAc                       >  2:31460/6‑149 (MQ=38)
                                                                                                                          aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                    <  2:268957/149‑1 (MQ=255)
                                                                                                                          aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                    <  3:135398/149‑1 (MQ=255)
                                                                                                                              acaAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                   >  2:243006/1‑130 (MQ=39)
                                                                                                                              acaAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                   <  1:243006/130‑1 (MQ=39)
                                                                                                                                aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCg              >  4:308169/1‑149 (MQ=255)
                                                                                                                                aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCg              >  3:88233/1‑149 (MQ=255)
                                                                                                                                  gACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGc                  <  1:50364/143‑1 (MQ=255)
                                                                                                                                  gACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGc                  >  2:50364/1‑143 (MQ=255)
                                                                                                                                       accgCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTg         >  2:22845/4‑147 (MQ=255)
                                                                                                                                       aCAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGa       <  4:293607/149‑1 (MQ=255)
                                                                                                                                       aCAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGa       <  5:31107/149‑1 (MQ=255)
                                                                                                                                         aGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAg      >  6:121626/1‑148 (MQ=255)
                                                                                                                                            cATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAcc  <  2:137041/149‑1 (MQ=255)
                                                                                                                                            cATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAc   >  8:209701/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACC  >  NC_002947/6126805‑6127093

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