Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,126,949 A→G 5.1% intergenic (+40/+27) 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,9490AG 5.1% 234.6 / 2.0 78intergenic (+40/+27)PP_5374/PP_5375choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (32/42);  new base G (2/2);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

TGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCAC  >  NC_002947/6126801‑6127095
                                                                                                                                                    |                                                                                                                                                  
tGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCa                                                                                                                                                    <  4:61657/149‑1 (MQ=255)
       cTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCg                                                                                                                                             <  3:107660/149‑1 (MQ=255)
           cTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCt                                                                                                                                         >  2:134604/1‑149 (MQ=255)
            tataCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCtt                                                                                                                                        >  7:344766/1‑149 (MQ=255)
              taCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCTtgt                                                                                                                                      <  8:3178/149‑1 (MQ=39)
              taCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTtgt                                                                                                                                      <  3:121159/149‑1 (MQ=255)
                 cTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGtt                                                                                                                                   >  6:18556/1‑149 (MQ=255)
                    ggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTg                                                                                                                                <  1:134604/149‑1 (MQ=255)
                     gATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTgg                                                                                                                               >  5:140746/1‑149 (MQ=255)
                     gATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTgg                                                                                                                               <  5:18556/149‑1 (MQ=255)
                      aTTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTggg                                                                                                                              >  6:30995/1‑149 (MQ=255)
                        tCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTc                                                                                                                                          >  4:57534/1‑135 (MQ=255)
                        tCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTc                                                                                                                                          <  3:57534/135‑1 (MQ=255)
                          tcccGCATTGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGTGGGGCTTGATATGTGGCTGTCAGCGAAGAAAGCGAGCTTTCCTTAAGCCACGCTGGAAGACACAAGACCGAATGCTGTTGGCTGTCGGTCTtgt                                                                                                                                      <  6:357330/136‑1 (MQ=255)
                          aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGCGCAGACAGCGAAGAAAGCCAGCTAACCTAAAGCCAAACTGGAAGACACAAGACAGACAGCTGATGGCTGTCGGACTtgt                                                                                                                                      >  5:357330/1‑137 (MQ=255)
                            ccGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGc                                                                                                                         >  6:294071/1‑148 (MQ=255)
                             cGCAAGGCATCACGCCGGTGGCTGGGCCGTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTa                                                                                                                       >  1:211962/1‑149 (MQ=255)
                                  ggCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCt                                                                                                                  >  3:254466/1‑149 (MQ=255)
                                  ggCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCt                                                                                                                  <  1:359175/149‑1 (MQ=255)
                                      tCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCAc                                                                                                              <  8:251887/149‑1 (MQ=255)
                                           ccGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAg                                                                                                         <  2:211962/149‑1 (MQ=255)
                                            cGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAgg                                                                                                        >  7:225940/1‑149 (MQ=255)
                                               tGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAg                                                                                                     >  6:274345/1‑149 (MQ=255)
                                                 gCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAggg                                                                                                   <  5:236553/149‑1 (MQ=255)
                                                 gCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAggg                                                                                                   <  8:128491/149‑1 (MQ=255)
                                                   tGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGcc                                                                                                 <  3:74596/149‑1 (MQ=255)
                                                   tGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGc                                                                                                  >  7:155642/1‑148 (MQ=255)
                                                        ccTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCaa                                                                                            >  6:363574/1‑149 (MQ=255)
                                                         cTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                           <  5:30995/149‑1 (MQ=255)
                                                            ggCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAcc                                                                                          <  1:165843/147‑1 (MQ=255)
                                                                  cGAGGGGGTTGATATGGGGCTGACATCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTGGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCgg                                                                                  <  5:43587/149‑1 (MQ=255)
                                                                  cGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTccc                                                                                    <  6:281002/147‑1 (MQ=255)
                                                                  cGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTccc                                                                                    <  5:363574/147‑1 (MQ=255)
                                                                         cTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATAAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGcc                                                                           >  1:368548/1‑149 (MQ=255)
                                                                                tGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACCGCCCTCAGCAGTCGGGGGTTGGGTTTGGGGACGCGAATTGAGCACCGAAGGCCGGGGCAGCAACCAGTGCCGGGACACCCCCCCGcg                                                                    >  6:55332/1‑148 (MQ=255)
                                                                                tGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTACCGGAACAGCCGCACGCt                                                                    >  1:258923/1‑149 (MQ=255)
                                                                                 ggggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTggg                                                                                                                              >  8:132814/1‑90 (MQ=255)
                                                                                 ggggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTggg                                                                                                                              <  7:132814/90‑1 (MQ=255)
                                                                                 ggggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGcc                                                                           >  7:355277/1‑141 (MQ=255)
                                                                                 ggggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGcc                                                                           <  8:355277/141‑1 (MQ=255)
                                                                                    gCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCa                                                                <  6:140746/149‑1 (MQ=255)
                                                                                        aCAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAg                                                            <  3:346744/149‑1 (MQ=255)
                                                                                         cAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAgg                                                           <  5:294071/149‑1 (MQ=255)
                                                                                         cAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAgg                                                           <  8:102859/149‑1 (MQ=255)
                                                                                          aGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGa                                                          >  7:233427/1‑149 (MQ=255)
                                                                                          aGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGa                                                          >  3:90595/1‑149 (MQ=255)
                                                                                          aGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGa                                                          >  6:145759/1‑149 (MQ=255)
                                                                                             gaagaaAGCCAGCTAACGTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGCACAGCCGCACGCTTTCACGAGGACCa                                                       >  4:252562/1‑149 (MQ=255)
                                                                                             gaagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCa                                                       <  8:233427/149‑1 (MQ=255)
                                                                                             gaagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCa                                                       >  8:106593/1‑149 (MQ=255)
                                                                                             gaagaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTCCTGCACCGAAGGCAggg                                                                                                   <  6:241546/105‑1 (MQ=255)
                                                                                                     ccAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCTTCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGt                                               <  1:277185/149‑1 (MQ=255)
                                                                                                     ccAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGt                                               <  4:254466/149‑1 (MQ=255)
                                                                                                         cTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                           <  6:32805/149‑1 (MQ=255)
                                                                                                         cTAACCTAAAGACACCCTGGAATAGACAAGACGGACAGCGATCATCTGTCGGTCTTGTGTTTTGGGCAGTCTATTTCTGCACCGAAGGCATGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACTCTCGGTGgac                                           <  4:348487/149‑1 (MQ=255)
                                                                                                          tAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                           >  3:240213/1‑148 (MQ=255)
                                                                                                            aCCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgacga                                         >  7:332606/1‑148 (MQ=255)
                                                                                                            aCCTAAAGCCACCCCGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgacga                                         <  1:205905/148‑1 (MQ=255)
                                                                                                              cTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTc                                      >  1:267051/1‑149 (MQ=255)
                                                                                                                   gCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCcagca                                 <  8:155642/149‑1 (MQ=255)
                                                                                                                       cccTGGAAGACACAAGAGCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcagg                               >  5:29926/1‑147 (MQ=255)
                                                                                                                        ccTGGAAGACACAAGACCGACAGCCATCAGCTGTGGGTCTTGTGTTTTGGGGAGGCTATTTCTGCACCGAAGGCAGGGCCCTCAACCAGTCCCGGAACAGCCGCACGGTTTCACGAGGACCACGCTCGGGGGCCGACTCCAgcagg                               >  3:231276/1‑146 (MQ=255)
                                                                                                                        ccTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAgcagg                               >  6:71519/1‑146 (MQ=255)
                                                                                                                           ggAAGACACAAGACCGACAGCCATCAGCCGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAAc                         >  1:163150/1‑149 (MQ=255)
                                                                                                                            gAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGaa                          <  2:267051/147‑1 (MQ=255)
                                                                                                                            gAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAAc                         <  8:225940/148‑1 (MQ=255)
                                                                                                                              aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                      <  2:258923/149‑1 (MQ=255)
                                                                                                                              aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                      <  7:106593/149‑1 (MQ=255)
                                                                                                                              aGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTg                      <  4:90595/149‑1 (MQ=255)
                                                                                                                                  acaAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGt                  >  8:227831/1‑149 (MQ=39)
                                                                                                                                   caAGCCCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCAGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTc                 <  8:344766/149‑1 (MQ=255)
                                                                                                                                    aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCg                <  3:279734/149‑1 (MQ=255)
                                                                                                                                           aCAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGa         <  6:29926/149‑1 (MQ=255)
                                                                                                                                           aCAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGa         <  5:71519/149‑1 (MQ=255)
                                                                                                                                               ccATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAc     >  6:278880/1‑149 (MQ=255)
                                                                                                                                               ccATCAGCCGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAc     <  2:163150/149‑1 (MQ=255)
                                                                                                                                                cATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTAcc    <  5:274345/149‑1 (MQ=255)
                                                                                                                                                 aTCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACca   >  1:216328/1‑149 (MQ=255)
                                                                                                                                                  tCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACcac  <  4:240213/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCAC  >  NC_002947/6126801‑6127095

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