Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,126,947 T→A 6.0% 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.0% 242.9 / 3.4 82intergenic (+38/+29)PP_5374/PP_5375choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base A (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

TGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCA  >  NC_002947/6126801‑6127094
                                                                                                                                                  |                                                                                                                                                   
tGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCa                                                                                                                                                   <  5:357788/149‑1 (MQ=255)
           cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAACCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTc                                                                                                                                         <  2:240680/148‑1 (MQ=255)
            tataCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCtt                                                                                                                                       >  8:98669/1‑149 (MQ=255)
             ataCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTtg                                                                                                                                      >  5:19855/1‑149 (MQ=255)
                 cTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCg                                                                                                                                            >  1:174621/1‑139 (MQ=255)
                 cTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCg                                                                                                                                            <  2:174621/139‑1 (MQ=255)
                    ggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTg                                                                                                                               <  2:353792/149‑1 (MQ=255)
                     gATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTgg                                                                                                                              <  5:181457/149‑1 (MQ=255)
                      aTTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTACCCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTggg                                                                                                                             >  6:346273/1‑149 (MQ=255)
                          aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAgg                                                                                                                         <  5:5905/149‑1 (MQ=255)
                             cGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTa                                                                                                                      >  2:169928/1‑149 (MQ=255)
                                      tCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTc                                                                                                                                         >  8:147829/1‑121 (MQ=39)
                                      tCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTc                                                                                                                                         <  7:147829/121‑1 (MQ=39)
                                      tCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCAc                                                                                                             <  6:81320/149‑1 (MQ=255)
                                            cGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGAAACAAGACCCACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCAATTTCTGCCCCGAAgg                                                                                                       >  4:339844/1‑149 (MQ=255)
                                              gTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCa                                                                                                     <  5:350057/149‑1 (MQ=255)
                                                 gCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAggg                                                                                                  <  7:98669/149‑1 (MQ=255)
                                                 gCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAgg                                                                                                   <  5:109825/148‑1 (MQ=255)
                                                   tGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGc                                                                                                 >  5:221212/1‑148 (MQ=255)
                                                      gCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAggg                                                                                                  <  1:63760/144‑1 (MQ=255)
                                                      gCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAggg                                                                                                  >  2:63760/1‑144 (MQ=255)
                                                               cTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTcc                                                                                    >  1:149116/1‑149 (MQ=255)
                                                                 gCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCg                                                                                  >  4:357614/1‑149 (MQ=255)
                                                                  cGAGGGGCTTGTTATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCgg                                                                                 <  1:136207/149‑1 (MQ=255)
                                                                  agtggggCTTGATTTGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACGCTGGAAGACACAAGCACGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCgg                                                                                 <  2:269215/146‑1 (MQ=255)
                                                                     ggggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGaa                                                                               <  2:349109/148‑1 (MQ=255)
                                                                        gCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGc                                                                           >  2:89614/1‑149 (MQ=255)
                                                                        gCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGc                                                                           >  2:172613/1‑149 (MQ=255)
                                                                         cTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGcc                                                                          >  4:126406/1‑149 (MQ=255)
                                                                         cTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAACACACAAGACCGACAGCCATCAGCTCTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGCGCCATCAACCAGTCCCGGAACAGcc                                                                          >  1:334315/1‑149 (MQ=255)
                                                                            gATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAg                                                                                                                                                  >  2:67612/1‑74 (MQ=255)
                                                                            gATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAg                                                                                                                                                  <  1:67612/74‑1 (MQ=255)
                                                                            gATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTATGCACCGAAGGCAGGGCCAACAACCAGTCCCGGAACAGCCCCa                                                                       >  4:375459/1‑149 (MQ=255)
                                                                             atatGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCAc                                                                      <  8:142324/149‑1 (MQ=255)
                                                                                tGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGc                                                                    <  3:357614/148‑1 (MQ=255)
                                                                                   ggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTc                                                                <  3:103363/149‑1 (MQ=255)
                                                                                       gACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCAc                                                              >  6:108079/1‑147 (MQ=255)
                                                                                        aCAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGa                                                            <  7:200718/148‑1 (MQ=255)
                                                                                          aGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGa                                                         >  4:312279/1‑149 (MQ=255)
                                                                                              aagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCAc                                                     >  3:66165/1‑149 (MQ=255)
                                                                                               agaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACg                                                    >  7:233111/1‑149 (MQ=255)
                                                                                               agaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACg                                                    >  5:80077/1‑149 (MQ=255)
                                                                                                 aaaGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCt                                                  <  5:192220/149‑1 (MQ=255)
                                                                                                       aGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTg                                             >  5:165701/1‑148 (MQ=255)
                                                                                                        gCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCgg                                               >  6:109850/1‑145 (MQ=255)
                                                                                                        gCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCgg                                               <  5:109850/145‑1 (MQ=255)
                                                                                                        gCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGga                                           <  5:108079/149‑1 (MQ=255)
                                                                                                        gCTAACCTAAAGCCACCCTGGAAGACACAACACCCAAAGCAATCAGCTGTGGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGga                                           >  8:284503/1‑149 (MQ=255)
                                                                                                         cTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                          <  3:312279/149‑1 (MQ=255)
                                                                                                             ccTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTAGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACt                                      >  5:268891/1‑149 (MQ=255)
                                                                                                              cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTAAGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCCAGCTTTCACGAGGACCACGCTAGGTGGACGACTc                                     >  2:356665/1‑149 (MQ=255)
                                                                                                              cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCATGTGTTTTGGGCAGGCAAATTCAGCACCGAAGGCAGGGCCATCAACCAGTACCGGAACAGCCGCACGCTTTCACCAGGACCACGCTCGGTGCACGACTc                                     >  2:131232/1‑149 (MQ=255)
                                                                                                               tAAAGCCACCCTGGAAGACACAAGACCCACAGCCATCAGCGGTCGGTCTTGTGTTTTGGGGAGGCTATTTATGCACCGAAGGCAGGGCCCTCAACCAGTCCCGGGACAGCCGCACGCTTTCACGAGGACCACCCTCGGTGgacgccgcc                                    >  2:186067/1‑144 (MQ=255)
                                                                                                                    ccACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGaa                                                                                                         <  5:106894/75‑1 (MQ=255)
                                                                                                                    ccACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGaa                                                                                                         >  6:106894/1‑75 (MQ=255)
                                                                                                                     cACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTTTGTTTTTGGCAGGCTTTTTCTGCACCGAAGGCAGGGCCATCAACCAATCCCGGAACAGCCGCACGCTTTACAGAGGACCACGCACGGTGGAAGACTACAtcagg                              >  8:35415/1‑149 (MQ=255)
                                                                                                                     cACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCGCGGTGGACGACTCCAgcagg                              >  2:228202/1‑149 (MQ=255)
                                                                                                                      aCCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCAcc                                                                                                            <  8:285261/70‑1 (MQ=255)
                                                                                                                      aCCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCAcc                                                                                                            >  7:285261/1‑70 (MQ=255)
                                                                                                                       cccTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggca                            <  5:346273/149‑1 (MQ=255)
                                                                                                                       cccTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggca                            <  6:45280/149‑1 (MQ=255)
                                                                                                                        ccTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcaggcag                           <  2:149116/149‑1 (MQ=255)
                                                                                                                        ccTGGAAGACACAAGACCGACAGCCATCAGCTGCCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAgg                                                          >  6:77165/1‑118 (MQ=255)
                                                                                                                        ccTGGAAGACACAAGACCGACAGCCATCAGCTGCCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAgg                                                          <  5:77165/118‑1 (MQ=255)
                                                                                                                            gAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGCCTCCAGCAGGCAGAACt                       >  3:176332/1‑149 (MQ=255)
                                                                                                                              aGACACAAGACCGACAGCCATCAGCTTTCGGTCTTGTGTTTTGGTCAGTCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGTAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                     <  7:35415/149‑1 (MQ=255)
                                                                                                                              aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                     <  1:169928/149‑1 (MQ=255)
                                                                                                                              aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                     >  3:58752/1‑149 (MQ=255)
                                                                                                                               gACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTgg                    <  1:228202/149‑1 (MQ=255)
                                                                                                                               gACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAACCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTgg                    <  6:221212/149‑1 (MQ=255)
                                                                                                                                 cacaAGACCGCCAGCCAACCGCTGGCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCTCCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCg                  <  4:239194/149‑1 (MQ=255)
                                                                                                                                 cacaAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCg                  >  2:44720/1‑149 (MQ=39)
                                                                                                                                 cacaAGACCGACAGCCATCAGCTGTCGGTGTTGTGTTTTGGGCAGGCTTTTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGGCGCACGCTTTCACGTGGACCACGCTCGGTGGACTACTCCAGCAGGCAGAACTTGGCg                  <  3:339844/149‑1 (MQ=255)
                                                                                                                                 cacaAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCg                  >  1:337664/1‑149 (MQ=255)
                                                                                                                                  acaAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGt                 >  8:135310/1‑149 (MQ=255)
                                                                                                                                    aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCACACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCg               >  5:122785/1‑149 (MQ=255)
                                                                                                                                      gACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGt             >  4:78193/1‑149 (MQ=255)
                                                                                                                                       aCCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCAGGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTc            <  7:135310/149‑1 (MQ=255)
                                                                                                                                          gACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTgg         <  1:172613/149‑1 (MQ=255)
                                                                                                                                            cAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAg       <  3:78193/149‑1 (MQ=255)
                                                                                                                                             aGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAg       >  7:47679/1‑148 (MQ=255)
                                                                                                                                             aGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGt      <  3:375459/149‑1 (MQ=255)
                                                                                                                                                  tCAGCTGGCGGTCTTTTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCACGTAAAAGCCGCACTCTTTCAAGAGTAACACGCTCGGTTGAAGACTCCCGAAGGCAGCACTTGGCGTCTGTCTGGAGTAcaa  >  1:103090/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCA  >  NC_002947/6126801‑6127094

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