Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,468,949 G→C 6.6% intergenic (+3/+25) PP_3956 → / ← betT‑III hypothetical protein/choline/carnitine/betaine transporter BetTC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,468,9490GC 6.6% 1265.1 / 56.5 377intergenic (+3/+25)PP_3956/betT‑IIIhypothetical protein/choline/carnitine/betaine transporter BetTC
Reads supporting (aligned to +/- strand):  ref base G (187/165);  new base C (4/21);  total (191/186)
Fisher's exact test for biased strand distribution p-value = 3.02e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 378 total aligned reads displayed.

GAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCAC  >  NC_002947/4468800‑4469093
                                                                                                                                                     |                                                                                                                                                
   gCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGa                                                                                                                                               <  1:6076228/150‑1 (MQ=255)
   gCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGa                                                                                                                                               <  2:8493158/150‑1 (MQ=255)
     tGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACa                                                                                                                                             >  1:8752687/1‑150 (MQ=255)
        tGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGTCCTGCCCTGTCGCCGTCAggg                                                                                                                                          <  1:2251312/150‑1 (MQ=255)
          aCTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGTCCTGCCCTGTCGCCGTCAGGGCa                                                                                                                                        <  1:8329582/150‑1 (MQ=255)
          aCTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCa                                                                                                                                        <  2:7700882/150‑1 (MQ=255)
          aCTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCa                                                                                                                                        <  1:3101985/150‑1 (MQ=255)
           cTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAg                                                                                                                                       <  2:7401229/150‑1 (MQ=255)
            tACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTCTTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAgg                                                                                                                                      >  2:712130/1‑150 (MQ=255)
            tACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAgg                                                                                                                                      <  2:4994110/150‑1 (MQ=255)
            tACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAgg                                                                                                                                      >  1:5451821/1‑150 (MQ=255)
             aCCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGCAggg                                                                                                                                     <  1:1241858/150‑2 (MQ=255)
             aCCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGa                                                                                                                                     <  1:7573226/150‑1 (MQ=255)
              ccAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGAc                                                                                                                                    >  2:3888135/1‑150 (MQ=255)
                aGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCg                                                                                                                                  >  2:8633596/1‑150 (MQ=255)
                 ggTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGa                                                                                                                                 <  2:4863657/150‑1 (MQ=255)
                  gTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGaa                                                                                                                                <  1:6763525/150‑1 (MQ=255)
                          cgaaGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCCCTTc                                                                                                                        >  2:4460304/4‑150 (MQ=255)
                                   gTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCccacca                                                                                                               <  2:3760184/150‑1 (MQ=255)
                                     gccagccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGTAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACg                                                                                                             <  1:729854/150‑1 (MQ=255)
                                        agccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGc                                                                                                          >  1:1763063/1‑150 (MQ=255)
                                         gccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCa                                                                                                         <  2:18992/150‑1 (MQ=255)
                                          ccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCaa                                                                                                        <  1:403052/150‑1 (MQ=255)
                                              tGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTaa                                                                                                    <  2:6373/150‑1 (MQ=255)
                                              tGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTaa                                                                                                    >  1:4979052/1‑150 (MQ=255)
                                                tCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATg                                                                                                  >  1:1454686/1‑150 (MQ=255)
                                                   tATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAg                                                                                               >  2:2467999/1‑150 (MQ=255)
                                                   tATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAg                                                                                               >  1:7755920/1‑150 (MQ=255)
                                                   tATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAg                                                                                               >  2:5714366/1‑150 (MQ=255)
                                                   tATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAg                                                                                               >  1:325413/1‑150 (MQ=255)
                                                   tATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAg                                                                                               >  2:1489787/1‑150 (MQ=255)
                                                        gTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCACGTAATGCAGGTGGc                                                                                          >  2:6693355/1‑150 (MQ=255)
                                                          tAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGt                                                                                        <  2:2235643/150‑1 (MQ=255)
                                                          tAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGt                                                                                        <  2:6828794/150‑1 (MQ=255)
                                                            ggACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGACCCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTc                                                                                      <  2:7581218/150‑1 (MQ=255)
                                                            ggACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTc                                                                                      <  1:1892865/150‑1 (MQ=255)
                                                            ggACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTc                                                                                      >  1:6991016/1‑150 (MQ=255)
                                                               ccgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTa                                                                                   <  2:5087822/150‑1 (MQ=255)
                                                               ccgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTa                                                                                   <  1:7840690/150‑1 (MQ=255)
                                                                cgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                                                  <  1:6531314/150‑1 (MQ=255)
                                                                 gccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACt                                                                                 >  2:4392782/1‑150 (MQ=255)
                                                                 gccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACt                                                                                 >  2:6967517/1‑150 (MQ=255)
                                                                 gccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGGTCGTACt                                                                                 >  2:789238/1‑150 (MQ=255)
                                                                   cgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGa                                                                               <  1:329110/150‑1 (MQ=255)
                                                                   cgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGa                                                                               <  2:4053436/150‑1 (MQ=255)
                                                                    gAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGAt                                                                              >  2:4076353/1‑150 (MQ=255)
                                                                       cGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCa                                                                           <  1:2401684/150‑1 (MQ=255)
                                                                        gAACGCATTGCTGCATTAGAATTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAg                                                                          <  1:998456/150‑1 (MQ=255)
                                                                         aaCGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAgg                                                                         >  1:7896758/1‑150 (MQ=255)
                                                                         aaCGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAgg                                                                         >  1:5768117/1‑150 (MQ=255)
                                                                          aCGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGa                                                                        <  2:7790202/150‑1 (MQ=255)
                                                                          aCGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGa                                                                        <  2:8375828/150‑1 (MQ=255)
                                                                             cATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGt                                                                     <  2:4753349/150‑1 (MQ=255)
                                                                              aTTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTc                                                                    >  1:5641092/1‑150 (MQ=255)
                                                                              aTTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTc                                                                    >  1:1196522/1‑150 (MQ=255)
                                                                                tgctgcAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGc                                                                  >  1:6275038/1‑150 (MQ=255)
                                                                                 gctgcAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGcc                                                                 >  2:5224452/1‑150 (MQ=255)
                                                                                 gctgcAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGcc                                                                 >  2:5459694/1‑150 (MQ=255)
                                                                                        tAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATcaca                                                          >  2:2288354/1‑150 (MQ=255)
                                                                                        tAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATcaca                                                          >  2:3557663/1‑150 (MQ=255)
                                                                                        tAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATcaca                                                          >  2:8114776/1‑150 (MQ=255)
                                                                                          taCTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGTATGTCGCCAATCACATc                                                        <  2:3359082/149‑1 (MQ=255)
                                                                                           aCTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCCAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCg                                                       >  2:7867161/1‑150 (MQ=255)
                                                                                             ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTc                                                     >  2:3216934/1‑150 (MQ=255)
                                                                                             ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTc                                                     >  1:2658293/1‑150 (MQ=255)
                                                                                             ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTc                                                     <  1:5096781/150‑1 (MQ=255)
                                                                                               cTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCAc                                                   >  1:4789041/1‑150 (MQ=255)
                                                                                               cTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCAc                                                   <  1:1563531/150‑1 (MQ=255)
                                                                                                    ccTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTc                                              <  1:2947517/150‑1 (MQ=255)
                                                                                                    ccTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTc                                              <  2:7151872/150‑1 (MQ=255)
                                                                                                      tgtgAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCa                                            <  1:6838715/150‑1 (MQ=255)
                                                                                                       gtgAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATc                                                                                                                              <  2:2907451/67‑1 (MQ=255)
                                                                                                       gtgAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCAc                                           >  1:3636234/1‑150 (MQ=255)
                                                                                                        tgAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCAcc                                          >  1:5967356/1‑150 (MQ=255)
                                                                                                         gAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCAccc                                         >  2:2356628/1‑150 (MQ=255)
                                                                                                          aaaaTCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGTCGCCGTCAGGGc                                                                                                                                         >  1:8166080/1‑53 (MQ=25)
                                                                                                              tCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCatga                                    >  2:6206794/1‑150 (MQ=255)
                                                                                                              tCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCatga                                    >  2:1233136/1‑150 (MQ=255)
                                                                                                               ccGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCatgat                                   <  2:3126594/150‑1 (MQ=255)
                                                                                                                 gCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGt                                 >  1:3899677/1‑150 (MQ=255)
                                                                                                                 gCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGt                                 >  1:408051/1‑150 (MQ=255)
                                                                                                                   cTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCg                               <  1:3549454/150‑1 (MQ=255)
                                                                                                                   cTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCg                               <  1:1690123/150‑1 (MQ=255)
                                                                                                                     ttGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTa                             <  2:1645289/150‑1 (MQ=255)
                                                                                                                      tGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAg                            <  1:2805711/150‑1 (MQ=255)
                                                                                                                        gACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTc                          >  2:8691948/1‑150 (MQ=255)
                                                                                                                           agcagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTg                       <  2:2818995/150‑1 (MQ=255)
                                                                                                                           agcagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTg                       <  1:3770960/150‑1 (MQ=255)
                                                                                                                            gcagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCTCGGCTCCACCCCATGATGTCGTAGTCCTgg                      >  2:6495967/1‑150 (MQ=255)
                                                                                                                            gcagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTgg                      <  1:5181188/150‑1 (MQ=255)
                                                                                                                             cagGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTCGTCCTGGc                     >  2:8287326/1‑150 (MQ=255)
                                                                                                                              agGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGcc                    >  2:3037668/1‑150 (MQ=255)
                                                                                                                              agGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGcc                    >  2:2509170/1‑150 (MQ=255)
                                                                                                                                gCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGcccc                  >  1:7098378/1‑150 (MQ=255)
                                                                                                                                      cGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTtctc            >  2:604889/1‑150 (MQ=255)
                                                                                                                                      cGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTtctc            <  2:1112196/150‑1 (MQ=255)
                                                                                                                                         cctgccctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCa         <  1:101671/150‑1 (MQ=255)
                                                                                                                                              cctgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGc    >  1:4388747/1‑150 (MQ=255)
                                                                                                                                                tgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCAc  >  1:2129175/1‑150 (MQ=255)
                                                                                                                                                tgACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCAc  <  1:1119834/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGCCCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCAC  >  NC_002947/4468800‑4469093

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