Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,128,925 T→A 6.7% 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,128,9250TA 6.7% 438.5 / 13.0 140intergenic (+38/+29)PP_5374/PP_5375choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (52/77);  new base A (3/7);  total (56/84)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

GGGTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCAC  >  NC_002947/6128776‑6129073
                                                                                                                                                     |                                                                                                                                                    
gggTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCAt                                                                                                                                                      <  1:1868092/150‑1 (MQ=255)
 ggTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATc                                                                                                                                                     <  2:3264195/150‑1 (MQ=255)
  gTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCa                                                                                                                                                    <  1:4905409/150‑1 (MQ=255)
   tGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAg                                                                                                                                                   <  1:907801/150‑1 (MQ=255)
        gcCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTc                                                                                                                                              <  1:2095639/150‑1 (MQ=255)
         ccTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCg                                                                                                                                             <  2:1428615/150‑1 (MQ=255)
           ttCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGt                                                                                                                                           <  2:907809/150‑1 (MQ=255)
           ttCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGt                                                                                                                                           <  1:5317170/150‑1 (MQ=255)
            tCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTc                                                                                                                                          <  2:1748997/150‑1 (MQ=255)
             ccgacaCCTGGGAGTCACCCGCAAGGCCTCACGCCGGTGGCGGGGCCCGCGGCCTGCGAGGGGCTTGATATCGGGCTGACAGCGAAGAAAGCCAGCTAACCAAAAGCCACCCTGGACGACACAAGACCGACAGCCATCAGCTGTCGGTCt                                                                                                                                         <  2:4778335/145‑1 (MQ=255)
              cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATCTGGGGCTGACAGCGAAGAAAGCCAGCTAACCTCAAGCCACCCTGGAAGACACAAGACCGACAGCCCTCAGCTGTCGGTcgg                                                                                                                                        >  2:363820/1‑148 (MQ=255)
              cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAACGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGACGGTCtt                                                                                                                                        <  2:5162904/150‑1 (MQ=255)
              cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACCGCCATCAGCTGTCGGGCtt                                                                                                                                        >  2:4149989/1‑150 (MQ=255)
              cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCtt                                                                                                                                        <  2:3511468/150‑1 (MQ=255)
              cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCtt                                                                                                                                        >  2:3231382/1‑150 (MQ=255)
              cTATACCTGGGATTCACCCGAAAGGCATCACGCCGGTGGCTGGGCCCTTGTCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTTAAGCCACCCTGGAAGACAAAAGACCGTCAGCCATCAGCTGTCGGTCtt                                                                                                                                        <  2:3016684/150‑1 (MQ=255)
              cTATACCTGGGAGTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCtt                                                                                                                                        <  2:3772074/150‑1 (MQ=255)
                ataCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCTtgt                                                                                                                                      <  2:2070488/150‑1 (MQ=255)
                  aCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCTTgtgt                                                                                                                                    <  1:297114/150‑1 (MQ=39)
                     tGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAAACACAAGACCGACAGCCATCAGCTGTCGGTc                                                                                                                                          >  1:3558946/1‑141 (MQ=255)
                             aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGGGGTGTGTTTTTGGGAGGc                                                                                                                         >  1:1808626/1‑150 (MQ=255)
                                cGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTAt                                                                                                                      >  2:1275046/1‑150 (MQ=255)
                                cGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACACTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTAt                                                                                                                      >  1:3859591/1‑150 (MQ=255)
                                agCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTAt                                                                                                                      >  1:3586852/2‑150 (MQ=255)
                                    aGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCt                                                                                                                  <  2:4973102/150‑1 (MQ=255)
                                     ggCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTg                                                                                                                 <  2:5133377/150‑1 (MQ=255)
                                       cATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTACAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCa                                                                                                               <  2:1449398/150‑1 (MQ=255)
                                         tCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAATACCGACAGCCATCCGCTGTCGGGCTTGTGTTTTGGGCAGGCTCTTTCTGCAcc                                                                                                             >  1:4811342/1‑150 (MQ=255)
                                           aCGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGa                                                                                                           >  1:4269963/1‑150 (MQ=255)
                                            cGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGaa                                                                                                          >  2:3397462/1‑150 (MQ=255)
                                              ccGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAgg                                                                                                        <  1:4807154/150‑1 (MQ=255)
                                                      tGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCGAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGGTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCa                                                                                                >  1:4872657/1‑150 (MQ=255)
                                                       gggCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCAt                                                                                               <  1:330295/150‑1 (MQ=255)
                                                        ggCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATc                                                                                              >  2:1213198/1‑150 (MQ=255)
                                                        ggCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATc                                                                                              <  2:4493999/150‑1 (MQ=255)
                                                          cccTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCaa                                                                                            <  1:1762905/150‑1 (MQ=255)
                                                             ttGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCa                                                                                         <  1:4028688/150‑1 (MQ=255)
                                                             ttGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCa                                                                                         <  1:654955/150‑1 (MQ=255)
                                                               ggCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGt                                                                                       <  1:2999068/150‑1 (MQ=255)
                                                                gCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTc                                                                                      <  2:534965/150‑1 (MQ=255)
                                                                   tGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCg                                                                                   <  2:1035903/150‑1 (MQ=255)
                                                                    gCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCgg                                                                                  <  2:4684668/150‑1 (MQ=255)
                                                                    gCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCgg                                                                                  <  2:4684798/150‑1 (MQ=255)
                                                                      gAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGaa                                                                                <  2:4045474/150‑1 (MQ=255)
                                                                        ggggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCtt                                                                                                                                        >  2:934012/1‑92 (MQ=255)
                                                                        ggggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCtt                                                                                                                                        <  1:934012/92‑1 (MQ=255)
                                                                        ggggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACa                                                                              <  1:1636432/150‑1 (MQ=255)
                                                                         gggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAg                                                                             <  2:1611989/150‑1 (MQ=255)
                                                                          ggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGc                                                                            <  1:5142550/150‑1 (MQ=255)
                                                                           gCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGcc                                                                           <  2:2310063/150‑1 (MQ=255)
                                                                                  atGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCt                                                                    <  1:1275551/150‑1 (MQ=255)
                                                                                           aCAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTATGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAgg                                                           >  1:3773982/1‑150 (MQ=255)
                                                                                              gCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGAcc                                                        <  2:2772822/150‑1 (MQ=255)
                                                                                               cGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGc                                                                                                       >  2:5321877/1‑102 (MQ=255)
                                                                                               cGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCa                                                       >  1:2690831/1‑150 (MQ=255)
                                                                                                gaagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCAc                                                      >  2:407816/1‑150 (MQ=255)
                                                                                                gaagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCGTGGGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCGTTCACGAGGACCAc                                                      >  2:233579/1‑150 (MQ=255)
                                                                                                 aagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGaa                                                                                <  1:4969180/123‑1 (MQ=255)
                                                                                                 aagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACg                                                     >  1:1294690/1‑150 (MQ=255)
                                                                                                   gaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCa                                                                <  2:4495104/137‑1 (MQ=255)
                                                                                                   gaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCa                                                                >  1:4495104/1‑137 (MQ=255)
                                                                                                      aGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATAAGCTGTCGGTCGTGTGTTTTGGGCAGGCTATTTCTTCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCAAGAGGACCACGCACgg                                                >  1:1427312/1‑150 (MQ=255)
                                                                                                       gCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGTCGCACGCTTTCACGAGGACCACGCTCGGt                                               <  1:2482014/150‑1 (MQ=255)
                                                                                                        ccAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTg                                              <  1:3397290/150‑1 (MQ=255)
                                                                                                         cAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTgg                                             <  1:2167668/150‑1 (MQ=255)
                                                                                                          aGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGga                                            <  2:1017522/150‑1 (MQ=255)
                                                                                                           gCTACCCTACAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                           <  2:3036458/150‑1 (MQ=255)
                                                                                                           gCTAACCTACCGCCACCCTGGAAGACACAAGACCGACCGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                           <  2:4854334/150‑1 (MQ=255)
                                                                                                           gCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                           <  2:4576169/150‑1 (MQ=255)
                                                                                                              aaCCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCAGCAGCTGTCGGTCGTGTGTTTTTGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGGAAAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgacgac                                        >  1:1096427/1‑150 (MQ=255)
                                                                                                                ccTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTc                                      <  2:2797806/150‑1 (MQ=255)
                                                                                                                ccTAAAGCCACCCTGGAAGACACAAGACCGACAGCCAACAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCAACAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTc                                      >  1:732424/1‑150 (MQ=255)
                                                                                                                 cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                     >  2:1772373/1‑150 (MQ=255)
                                                                                                                 cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                     <  2:3773982/150‑1 (MQ=255)
                                                                                                                  tAAAGCCACCCTGGAAGACACATGACCGACCGCCATCAGCTGTCGGTCTTGTGGTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAGCCGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTccg                                    >  2:2140605/1‑149 (MQ=255)
                                                                                                                     aGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCcagca                                 <  1:1766624/150‑1 (MQ=255)
                                                                                                                       ccACCCTGGAAGACACAAGACCGACAGCCCAAAGCGGTCGGTCTTGTGTTTTTGGCAGGCGATTTCTTCACCCAAAGCAGGGCCATCAAACACGCCCCGAACAGCCCCACCCTTTAACGAGGCCCACGCACCGTGGACGACGCCAgcagg                               >  1:616843/1‑150 (MQ=255)
                                                                                                                        cACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggc                              <  1:1216792/150‑1 (MQ=255)
                                                                                                                         aCCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggca                             <  2:3563755/150‑1 (MQ=255)
                                                                                                                         aCCCTGGAAGACACAAGACCGACAGCCATCAGCGGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggca                             <  2:3202444/150‑1 (MQ=255)
                                                                                                                          cccTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggcag                            >  2:3588742/1‑150 (MQ=255)
                                                                                                                            cTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGaa                          >  2:2276387/1‑150 (MQ=255)
                                                                                                                                aaaacacAAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                      >  1:3732668/4‑150 (MQ=39)
                                                                                                                                aaGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                           >  2:1587723/1‑81 (MQ=255)
                                                                                                                                aaGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                           <  1:1587723/81‑1 (MQ=255)
                                                                                                                                 aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTgg                     <  1:3162914/150‑1 (MQ=255)
                                                                                                                                 aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTgg                     <  2:2804589/150‑1 (MQ=255)
                                                                                                                                     acaAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCaa                                                                                            >  2:4542473/1‑75 (MQ=255)
                                                                                                                                     acaAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCaa                                                                                            <  1:4542473/75‑1 (MQ=255)
                                                                                                                                     acaAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGAAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTc                 >  1:2058925/1‑150 (MQ=255)
                                                                                                                                      caAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCg                >  2:4647970/1‑150 (MQ=255)
                                                                                                                                        aGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGt              >  1:1743782/1‑150 (MQ=255)
                                                                                                                                          aCCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCt            <  1:4585335/150‑1 (MQ=255)
                                                                                                                                           ccGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTg           >  2:1278200/1‑150 (MQ=255)
                                                                                                                                            cGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTgg          >  2:1427243/1‑150 (MQ=255)
                                                                                                                                            cGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTgg          >  1:1097823/1‑150 (MQ=255)
                                                                                                                                             gACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGa         <  1:4499971/150‑1 (MQ=255)
                                                                                                                                                 gCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAc     >  1:1055279/1‑150 (MQ=255)
                                                                                                                                                  ccATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAAc                                                                               <  1:4589803/75‑1 (MQ=255)
                                                                                                                                                  ccATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAcc    >  1:1480153/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCAC  >  NC_002947/6128776‑6129073

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