Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,345 A→T 5.7% intergenic (+30/+51) PP_3628 → / ← PP_3629 choline/carnitine/betaine transporter family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,125,3450AT 5.7% 276.7 / 6.5 87intergenic (+30/+51)PP_3628/PP_3629choline/carnitine/betaine transporter family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (42/40);  new base T (2/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTGGC  >  NC_002947/4125198‑4125493
                                                                                                                                                   |                                                                                                                                                    
gATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGACGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAag                                                                                                                                                    >  1:639734/1‑150 (MQ=255)
   cGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcag                                                                                                                                                 <  1:1803544/150‑1 (MQ=255)
    gagaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                                >  1:721035/1‑150 (MQ=255)
    gagaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGGCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                                <  1:488132/150‑1 (MQ=255)
        aGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGt                                                                                                                                            <  2:1904108/150‑1 (MQ=255)
         ggCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGTCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGtg                                                                                                                                           <  2:1662011/150‑2 (MQ=255)
         ggCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGAAGGAGAAGAAGGCATCGCCCAAGACAGCGGCAGCGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGcccctg                                                                                                                                           <  2:678226/150‑4 (MQ=255)
            atctcTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTcgc                                                                                                                                        <  1:530765/146‑1 (MQ=255)
             gTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTcgcg                                                                                                                                       >  1:1429026/1‑150 (MQ=255)
              tACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGa                                                                                                                                      >  2:2092713/1‑150 (MQ=255)
                 tGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGc                                                                                                                                   >  2:2238845/1‑150 (MQ=255)
                 tGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGc                                                                                                                                   <  2:1501482/150‑1 (MQ=255)
                 tGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGc                                                                                                                                   <  2:1632253/150‑1 (MQ=255)
                   tATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGTCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCga                                                                                                                                 <  1:1078053/150‑1 (MQ=255)
                         gaggTGAAGCATGGGGAGGAGGACTACATTGTGGTGCAGGAGAAGAAGGCATCGAACAAGACAACAGCAGAGGCCGAAGTAACGAGCAACTGACGGCGTCAACGAGAAAACCACGGGCTGGTAAGCAAACCGTTCGCTGTGCgagtcagc                                                                                                                           >  2:978461/3‑145 (MQ=255)
                           ggTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCgcgc                                                                                                                         >  2:2289690/1‑150 (MQ=255)
                             tGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGtg                                                                                                                                           >  1:1912583/1‑129 (MQ=255)
                             tGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCGtg                                                                                                                                           <  2:1912583/130‑2 (MQ=255)
                               aaGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCctgc                                                                                                                     <  1:2392699/150‑1 (MQ=255)
                                 gCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCctgctg                                                                                                                   >  2:1501812/1‑150 (MQ=255)
                                 gCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGACAGATGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCctgctg                                                                                                                   >  2:329962/1‑150 (MQ=255)
                                     ggggAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAAc                                                                                                               >  1:1547084/1‑150 (MQ=255)
                                      gggAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACa                                                                                                              >  2:1917018/1‑150 (MQ=255)
                                       ggaggaGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAGGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACAt                                                                                                             <  2:1741220/150‑1 (MQ=255)
                                       ggaggaGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGaaa                                                                                                                >  2:1283453/1‑147 (MQ=255)
                                       ggaggaGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGaaa                                                                                                                <  1:1283453/147‑1 (MQ=255)
                                         aggagTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATgg                                                                                                           <  1:981699/150‑1 (MQ=255)
                                          ggagTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGa                                                                                                          >  1:700226/1‑150 (MQ=255)
                                           gagTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGAt                                                                                                         <  2:573778/150‑1 (MQ=255)
                                                cTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCt                                                                                                    >  2:367599/1‑150 (MQ=255)
                                                 tACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCtt                                                                                                   >  1:1130259/1‑150 (MQ=255)
                                                     ttGTGGTGCAGGAGAAGAAGGCATCGGCGAATGCAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGACTTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAAc                                                                                               <  1:761668/150‑1 (MQ=255)
                                                             cAGGAGCCGAAGGCGTCGACCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                       >  2:1353642/1‑150 (MQ=255)
                                                             cAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGTATCTGCTTCAACCTTCCAAc                                                                                       <  2:1698913/150‑1 (MQ=255)
                                                             cAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                       <  2:487411/150‑1 (MQ=255)
                                                             cAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                       <  1:885024/150‑1 (MQ=255)
                                                              aGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                                <  2:1909804/92‑1 (MQ=255)
                                                              aGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                                >  1:1909804/1‑92 (MQ=255)
                                                                gagaAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgcc                                                                                    >  2:887643/1‑150 (MQ=255)
                                                                gagaAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgcc                                                                                    >  2:1968394/1‑150 (MQ=255)
                                                                   aagaagGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccgcc                                                                                 >  1:2287300/1‑150 (MQ=255)
                                                                   aagaagGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccgcc                                                                                 >  1:826757/1‑150 (MQ=255)
                                                                     gaagGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccgccgc                                                                               >  2:415432/1‑150 (MQ=255)
                                                                       agGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTa                                                                             <  1:1968394/150‑1 (MQ=255)
                                                                           aTCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCaa                                                                                                >  1:1415861/1‑127 (MQ=255)
                                                                           aTCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCaa                                                                                                <  2:1415861/127‑1 (MQ=255)
                                                                           aTCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCtt                                                                         >  1:1314649/1‑150 (MQ=255)
                                                                            tCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCTCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        <  1:2238845/150‑1 (MQ=255)
                                                                            tCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        <  1:1036854/150‑1 (MQ=255)
                                                                            tCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        <  2:2312098/150‑1 (MQ=255)
                                                                            gcGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        <  2:1146863/149‑1 (MQ=255)
                                                                              ggCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGc                                                                      <  1:1940145/150‑1 (MQ=255)
                                                                               gCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGcc                                                                     >  2:2242703/1‑150 (MQ=255)
                                                                                       aGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGAt                                                             <  1:329962/150‑1 (MQ=255)
                                                                                         cGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCt                                                           >  2:1292869/1‑150 (MQ=255)
                                                                                          ggCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCtt                                                          <  1:1292869/150‑1 (MQ=255)
                                                                                          ggCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCtt                                                          <  2:1314649/150‑1 (MQ=255)
                                                                                            cAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGc                                                        <  2:700226/150‑1 (MQ=255)
                                                                                             agagGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGcc                                                       >  2:805359/1‑150 (MQ=255)
                                                                                                   cGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCGGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTg                                                 >  1:1081254/1‑150 (MQ=255)
                                                                                                     aaGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCtt                                                          >  1:1873292/1‑139 (MQ=255)
                                                                                                     aaGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCtt                                                          <  2:1873292/139‑1 (MQ=255)
                                                                                                     aaGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCGATCTGca                                               >  1:1037358/1‑150 (MQ=255)
                                                                                                       gTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGc                                                                                                     <  1:1633393/94‑1 (MQ=255)
                                                                                                       gTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGc                                                                                                     >  2:1633393/1‑94 (MQ=255)
                                                                                                       gTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGcatc                                             <  2:721035/150‑1 (MQ=255)
                                                                                                          aCCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATg                                          >  1:8862/1‑150 (MQ=255)
                                                                                                           cccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATgg                                         <  1:805359/150‑1 (MQ=255)
                                                                                                             cGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGcc                                       <  1:367599/150‑1 (MQ=255)
                                                                                                                  cTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGaa                                  >  1:2367328/1‑150 (MQ=255)
                                                                                                                   tGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAac                                 <  1:1266514/150‑1 (MQ=255)
                                                                                                                       ggCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAacgcac                             >  2:1607183/1‑150 (MQ=255)
                                                                                                                         cGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAacgcacgc                           <  1:2092713/150‑1 (MQ=255)
                                                                                                                          gTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        <  2:2214293/104‑1 (MQ=255)
                                                                                                                          gTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        >  1:2214293/1‑104 (MQ=255)
                                                                                                                                 cGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGaa                                  >  2:512984/1‑135 (MQ=255)
                                                                                                                                 cGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGaa                                  <  1:512984/135‑1 (MQ=255)
                                                                                                                                  gAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgccgcc                  >  2:2213366/1‑150 (MQ=255)
                                                                                                                                  gAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGgccgcc                  >  2:828793/1‑150 (MQ=255)
                                                                                                                                        cACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACgg            >  1:48620/1‑150 (MQ=255)
                                                                                                                                        aaCGGGCTGCTTAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACgg            >  1:1609921/2‑150 (MQ=255)
                                                                                                                                            ggCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCgg        <  1:887643/150‑1 (MQ=255)
                                                                                                                                             gctgctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCggg       >  1:919232/1‑150 (MQ=255)
                                                                                                                                             gctgctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCggg       <  2:1429026/150‑1 (MQ=255)
                                                                                                                                                gctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCtt                                                                         >  2:579358/1‑81 (MQ=255)
                                                                                                                                                gctAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCtt                                                                         <  1:579358/81‑1 (MQ=255)
                                                                                                                                                  tAAGAAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTGGc  >  2:1429979/1‑150 (MQ=255)
                                                                                                                                                   aaGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        >  2:671033/1‑79 (MQ=255)
                                                                                                                                                   aaGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTg                                                                        <  1:671033/79‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTGGC  >  NC_002947/4125198‑4125493

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