Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,131,952 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,131,9520TA 6.7% 802.9 / 26.8 252intergenic (+38/+29)PP_5374/PP_5375choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (114/120);  new base A (7/11);  total (121/131)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

GGGTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCAC  >  NC_002947/6131803‑6132100
                                                                                                                                                     |                                                                                                                                                    
 ggTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATc                                                                                                                                                     <  1:5638919/150‑1 (MQ=255)
 ggTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATc                                                                                                                                                     <  1:4167283/150‑1 (MQ=255)
   tGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAg                                                                                                                                                   <  2:6830477/150‑1 (MQ=255)
       cgcCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGt                                                                                                                                               >  2:4351588/1‑150 (MQ=255)
        gcCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTc                                                                                                                                              <  2:3058118/150‑1 (MQ=255)
         ccTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCg                                                                                                                                             <  2:6193622/150‑1 (MQ=255)
            tCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTc                                                                                                                                          <  1:6377388/150‑1 (MQ=255)
            tCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTc                                                                                                                                          <  2:1557963/150‑1 (MQ=255)
             ccTATACCTGGGATTCACCCGCAAGTCATCACGCCGGTGGCCGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTTACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCt                                                                                                                                         <  1:771439/150‑1 (MQ=255)
               tataCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTtg                                                                                                                                       >  2:3417748/1‑150 (MQ=255)
                atcccagggCTTCACCCGCAAGGCAACCCGCCGGTGGCTGGGCCCGTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTtgt                                                                                                                                      <  1:3650390/144‑1 (MQ=255)
                   ccTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGtt                                                                                                                                   >  1:2510303/1‑150 (MQ=255)
                    cTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATAAGCTGTCGGTGTTGTGttt                                                                                                                                  >  1:592940/1‑150 (MQ=255)
                      gggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTg                                                                                                                                <  1:3937798/150‑1 (MQ=255)
                      gggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTg                                                                                                                                <  2:4259184/150‑1 (MQ=255)
                      gggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTg                                                                                                                                <  2:4619406/150‑1 (MQ=255)
                           tCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCGTGTGTTTTGGGCAg                                                                                                                           >  2:5884061/1‑150 (MQ=255)
                             aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGc                                                                                                                         >  2:6576777/1‑150 (MQ=255)
                             aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGc                                                                                                                         >  2:258392/1‑150 (MQ=255)
                              cccGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCt                                                                                                                        >  2:6665119/1‑150 (MQ=255)
                               ccGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTa                                                                                                                       >  1:3441915/1‑150 (MQ=255)
                               ccGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTGTTGGGCAGGCTa                                                                                                                       >  2:3461608/1‑150 (MQ=255)
                                cGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTAt                                                                                                                      >  1:2071657/1‑150 (MQ=255)
                                 gCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCTAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTAtt                                                                                                                     <  1:4409764/150‑1 (MQ=255)
                                 gCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTAtt                                                                                                                     >  1:1847950/1‑150 (MQ=255)
                                   aaGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTc                                                                                                                   >  1:535111/1‑150 (MQ=255)
                                      gCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGc                                                                                                                <  2:1982536/150‑1 (MQ=255)
                                           aCGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGa                                                                                                           >  2:5146563/1‑150 (MQ=255)
                                              ccGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAgg                                                                                                        <  1:1065356/150‑1 (MQ=255)
                                                ggTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCa                                                                                                      <  1:268263/150‑1 (MQ=255)
                                                ggTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCa                                                                                                      <  2:6308951/150‑1 (MQ=255)
                                                    gCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGc                                                                                                  <  1:1789226/150‑1 (MQ=255)
                                                    gCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGc                                                                                                  <  1:1687043/150‑1 (MQ=255)
                                                       gggCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCAt                                                                                               <  2:5179748/150‑1 (MQ=255)
                                                          cccTTGGCCTGCGAGGGGTTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCaa                                                                                            >  2:3874834/1‑150 (MQ=255)
                                                           ccTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTgg                                                                                                                               >  2:4936271/1‑114 (MQ=255)
                                                               ggCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGt                                                                                       <  2:3127873/150‑1 (MQ=255)
                                                                gCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGtn                                                                                      <  1:1196473/150‑2 (MQ=255)
                                                                gCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTc                                                                                      <  2:4525026/150‑1 (MQ=255)
                                                                  cTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTGTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTccc                                                                                    >  1:4356189/1‑150 (MQ=255)
                                                                        ggggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACa                                                                              >  1:6799420/1‑150 (MQ=255)
                                                                        ggggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACACCCAACAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCCCCGAAGGCAGGGCCAACAACCACTCCCCGAACa                                                                              >  1:5761125/1‑150 (MQ=255)
                                                                           gCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGcc                                                                           <  2:6652264/150‑1 (MQ=255)
                                                                             ttGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGc                                                                         >  2:3602856/1‑150 (MQ=255)
                                                                              tGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCGTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCa                                                                        >  1:2677521/1‑150 (MQ=255)
                                                                                 tatGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGc                                                                     >  1:6149292/1‑150 (MQ=255)
                                                                                     gggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTc                                                                 <  2:2804759/150‑1 (MQ=255)
                                                                                     gggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTc                                                                 >  2:6811560/1‑150 (MQ=255)
                                                                                     gggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCAACAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTATGCACCGAAGGCAGGGCCAACAACCAGTCCCGGAACAGCCCCACCCTTTc                                                                 >  1:3481708/1‑150 (MQ=255)
                                                                                      ggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCa                                                                <  1:6665119/150‑1 (MQ=255)
                                                                                      ggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCa                                                                <  1:3864050/150‑1 (MQ=255)
                                                                                         tGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGa                                                             <  1:890003/150‑1 (MQ=255)
                                                                                         tGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCGGTCGGTCGTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCGTTCACGa                                                             >  2:164468/1‑150 (MQ=255)
                                                                                          gACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAg                                                            <  2:5772717/150‑1 (MQ=255)
                                                                                           aCAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAgg                                                           <  1:3579140/150‑1 (MQ=255)
                                                                                             aGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGAc                                                         >  2:523259/1‑150 (MQ=255)
                                                                                              gCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGAcc                                                        >  2:3903129/1‑150 (MQ=255)
                                                                                                 aagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCACCTGTCGGTCTTGTGTTTTTGGCACGCTATTTCTGCACCGAAGGCAGGGCCATCAAACAGTCCCCGAAAACCCCCAACCTTTCACGAGGACCACg                                                     >  1:4178320/1‑150 (MQ=255)
                                                                                                    aaaGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTc                                                  >  2:4139623/1‑150 (MQ=255)
                                                                                                      aGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCgg                                                <  1:6433617/150‑1 (MQ=255)
                                                                                                       gCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGt                                               <  2:6521156/150‑1 (MQ=255)
                                                                                                       gCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGt                                               <  2:326559/150‑1 (MQ=255)
                                                                                                         cAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTgg                                             >  2:861305/1‑150 (MQ=255)
                                                                                                                 cTAAAGCCACCCTGGAAGACCCAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCGTTCACGAGGCCCACGCTCGGTGGACGACTcc                                     >  1:4958197/1‑150 (MQ=255)
                                                                                                                 cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                     >  2:3682942/1‑150 (MQ=255)
                                                                                                                 cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                     >  2:379714/1‑150 (MQ=255)
                                                                                                                 cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                     >  1:2151575/1‑150 (MQ=255)
                                                                                                                 cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCGTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCCTCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTcc                                     >  2:2609036/1‑150 (MQ=255)
                                                                                                                  tAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCca                                    <  1:3901319/150‑1 (MQ=255)
                                                                                                                  tAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCca                                    <  1:480912/150‑1 (MQ=255)
                                                                                                                    aaGCCACCCTGGAAGACACAAGACCAACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCcagc                                  >  2:29629/1‑150 (MQ=255)
                                                                                                                         aCCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggca                             >  2:4462984/1‑150 (MQ=255)
                                                                                                                                aaGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                      <  1:1983264/150‑1 (MQ=255)
                                                                                                                                aaGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                      <  1:2366209/150‑1 (MQ=255)
                                                                                                                                  gACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGc                    <  1:499283/150‑1 (MQ=255)
                                                                                                                                  gACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGc                    <  1:6283838/150‑1 (MQ=255)
                                                                                                                                  gACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGc                    <  2:6035933/150‑1 (MQ=255)
                                                                                                                                   acacAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTggtg                   >  2:1470536/1‑148 (MQ=255)
                                                                                                                                     acaAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTc                 >  2:1807301/1‑150 (MQ=255)
                                                                                                                                      caAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCg                <  2:1055828/150‑1 (MQ=255)
                                                                                                                                       aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGga                                            <  2:5626159/121‑1 (MQ=255)
                                                                                                                                       aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCgg               >  2:2487740/1‑150 (MQ=255)
                                                                                                                                       aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCgg               >  1:3376604/1‑150 (MQ=255)
                                                                                                                                         gACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCg                >  1:4465952/1‑147 (MQ=255)
                                                                                                                                         gACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTc             >  1:3057438/1‑150 (MQ=255)
                                                                                                                                           ccGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTg           <  1:6576777/150‑1 (MQ=255)
                                                                                                                                             gACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGa         <  1:4909167/150‑1 (MQ=255)
                                                                                                                                             gACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGa         <  1:4520656/150‑1 (MQ=255)
                                                                                                                                              aCAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAg        >  1:325463/1‑150 (MQ=255)
                                                                                                                                              aCAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAg        <  2:1359810/150‑1 (MQ=255)
                                                                                                                                                aGCCATCAGCTGTCGGTCGTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTa      >  1:16204/1‑150 (MQ=255)
                                                                                                                                                 gCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAc     <  2:1748864/150‑1 (MQ=255)
                                                                                                                                                  ccATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGATCTTGGCGTCGGTCTGGAGTAcc    <  2:4820997/150‑1 (MQ=255)
                                                                                                                                                  ccATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAcc    >  1:4699337/1‑150 (MQ=255)
                                                                                                                                                  ccATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAcc    >  1:133786/1‑150 (MQ=255)
                                                                                                                                                  ccATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAcc    >  1:5171890/1‑150 (MQ=255)
                                                                                                                                                   cATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACca   >  1:6874626/1‑150 (MQ=255)
                                                                                                                                                    aTCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACcac  >  2:4289071/1‑150 (MQ=255)
                                                                                                                                                    aTCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACcac  <  1:591381/150‑1 (MQ=255)
                                                                                                                                                    aTCAGCTGTCGGTCTTGTGTTTTGGGAAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGACGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACcac  >  1:3841836/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCAC  >  NC_002947/6131803‑6132100

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: