Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,038,952 G→A 5.1% G236D (GGC→GAC)  PP_3559 → glycine betaine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,038,9520GA 5.1% 422.1 / 8.0 137G236D (GGC→GAC) PP_3559glycine betaine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (58/72);  new base A (3/4);  total (61/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 138 total aligned reads displayed.

GTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCAGTCGCTGAGCCG  >  NC_002947/4038803‑4039098
                                                                                                                                                     |                                                                                                                                                  
gTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTCTCGATGGTGGTCGTGGCGTCGAGGATCgg                                                                                                                                                    >  1:947157/1‑150 (MQ=255)
     tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                               >  2:498159/1‑150 (MQ=255)
     tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                               >  1:1762982/1‑150 (MQ=255)
     tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                               >  1:1715590/1‑150 (MQ=255)
     tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGCCGATCATCGACGcca                                                                                                                                               <  1:415242/150‑2 (MQ=255)
     ggcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                               <  2:1535518/149‑1 (MQ=255)
        aGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTg                                                                                                                                            <  2:1762982/150‑1 (MQ=255)
          gCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGa                                                                                                                                          >  1:762857/1‑150 (MQ=255)
            cGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGAtt                                                                                                                                        <  1:1724291/150‑1 (MQ=255)
            cGCACACGGCCTTGGCGCAAGCCGCTGGCACCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGAtt                                                                                                                                        <  1:1209796/150‑1 (MQ=255)
              cacaCGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTgg                                                                                                                                      <  1:378127/150‑1 (MQ=255)
                 acGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTg                                                                                                                                   >  1:1895979/1‑150 (MQ=255)
                  cGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGa                                                                                                                                  >  2:2064642/1‑150 (MQ=255)
                   ggCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGaa                                                                                                                                 >  1:2030766/1‑150 (MQ=255)
                      cTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGAt                                                                                                                              <  2:156233/150‑1 (MQ=255)
                          gcgcAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGc                                                                                                                          <  2:2180074/150‑1 (MQ=255)
                            gcAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTg                                                                                                                        >  1:1409683/1‑150 (MQ=255)
                            gcAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGCCGATCATCGACGcca                                                                                                                                               <  1:1542686/127‑2 (MQ=255)
                             cAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTgg                                                                                                                       <  2:1577992/150‑1 (MQ=255)
                                gCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGCAGATGTGCTGGCgg                                                                                                                    >  2:224667/1‑150 (MQ=255)
                                gCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCgg                                                                                                                    <  1:866328/150‑1 (MQ=255)
                                gCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCgg                                                                                                                    <  1:2154961/150‑1 (MQ=255)
                                    cTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAAt                                                                                                                <  1:1045989/150‑1 (MQ=255)
                                     tGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTCGCTTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATc                                                                                                               <  1:2242756/150‑1 (MQ=255)
                                        cAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGAAGAGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAg                                                                                                            >  2:942636/1‑150 (MQ=255)
                                               tGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGa                                                                                                     <  2:1409683/150‑1 (MQ=255)
                                                 cggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAAt                                                                                                   <  1:1798379/150‑1 (MQ=255)
                                                  ggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCt                                                                                                                         <  1:119058/127‑1 (MQ=255)
                                                   gcgcATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGt                                                                                                 >  2:2378237/1‑150 (MQ=255)
                                                     gcATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTg                                                                                               <  2:1405592/150‑1 (MQ=255)
                                                     gcATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTg                                                                                               <  2:936107/150‑1 (MQ=255)
                                                       aTCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGc                                                                                             >  1:566916/1‑150 (MQ=255)
                                                          gCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCaa                                                                                          >  1:774552/1‑150 (MQ=255)
                                                          gCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTAGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGCCCTTGAAGGTTGGCCaa                                                                                          <  2:897207/150‑1 (MQ=255)
                                                           cccTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATcc                                                                                                              <  1:241065/129‑1 (MQ=255)
                                                           cccTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAg                                                                                         <  2:1639869/150‑1 (MQ=255)
                                                            ccTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCCTGTCGAGGGTGGTCGTGGTGTCGATGAGCGGCGCCCGTGGATTTCGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAgg                                                                                        <  1:1855536/150‑1 (MQ=255)
                                                             cTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATTTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAggg                                                                                       >  1:2016861/1‑150 (MQ=255)
                                                             cTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCg                                                                                                                     <  2:2001996/120‑1 (MQ=255)
                                                             cTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCg                                                                                                                     >  1:2001996/1‑120 (MQ=255)
                                                                ccccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATg                                                                                    >  1:959953/1‑150 (MQ=255)
                                                                 cccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATgg                                                                                   <  2:1037777/150‑1 (MQ=255)
                                                                  ccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTg                                                                                                                                       >  2:617289/1‑97 (MQ=255)
                                                                  ccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTg                                                                                                                                       <  1:617289/97‑1 (MQ=255)
                                                                  ccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGa                                                                                  <  2:959953/150‑1 (MQ=255)
                                                                        gcTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCgg                                                                            <  2:2381762/150‑1 (MQ=255)
                                                                        gcTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGAAAGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCgg                                                                            >  1:637179/1‑150 (MQ=255)
                                                                          tACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGt                                                                                                                                    >  2:475272/1‑92 (MQ=255)
                                                                          tACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGt                                                                                                                                    <  2:1905893/92‑1 (MQ=255)
                                                                          tACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGt                                                                                                                                    >  1:1905893/1‑92 (MQ=255)
                                                                          tACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGt                                                                                                                                    <  1:475272/92‑1 (MQ=255)
                                                                             cAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCg                                                                             >  2:1320574/1‑144 (MQ=255)
                                                                             cAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCg                                                                             <  1:1320574/144‑1 (MQ=255)
                                                                              aaGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGc                                                                                                                      <  2:1003585/102‑1 (MQ=255)
                                                                              aaGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGc                                                                                                                      >  1:1003585/1‑102 (MQ=255)
                                                                              aaGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGAGGGAAGCGGGGCTgtc                                                                      >  2:625492/1‑148 (MQ=255)
                                                                                gCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCaa                                                                    >  2:1918551/1‑150 (MQ=255)
                                                                                  atcatGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAAtc                                                                  <  1:224667/150‑2 (MQ=255)
                                                                                     atGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATg                                                                                                  <  2:1085816/115‑1 (MQ=255)
                                                                                     atGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTc                                                               <  2:1492714/150‑1 (MQ=255)
                                                                                          cGGGCTCAACCAGTCGGGCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCt                                                          <  2:762857/150‑1 (MQ=255)
                                                                                           gggCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTg                                                         <  1:2378237/150‑1 (MQ=255)
                                                                                             gCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGc                                                       >  1:553583/1‑150 (MQ=255)
                                                                                               tCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCa                                                     >  2:989103/1‑150 (MQ=255)
                                                                                               tCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGACGTGCTGGCGGGAATCCAGACCTTGAATGTTAGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCa                                                     >  1:1364717/1‑150 (MQ=255)
                                                                                                  aCCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTgg                                                                       <  1:2041120/129‑1 (MQ=255)
                                                                                                  aCCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTgg                                                                       >  2:2041120/1‑129 (MQ=255)
                                                                                                   ccAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGt                                                 >  1:577642/1‑150 (MQ=255)
                                                                                                    cAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTgg                                                                       <  2:146899/127‑1 (MQ=255)
                                                                                                    cAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTgg                                                                       >  1:146899/1‑127 (MQ=255)
                                                                                                    cAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTg                                                >  1:1176903/1‑150 (MQ=255)
                                                                                                       gcGGGCATGACGCCCTGGGCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATc                                             <  2:1715590/149‑1 (MQ=255)
                                                                                                         ggTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTgg                                                        >  1:996332/1‑137 (MQ=255)
                                                                                                          gTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGAc                                          >  1:811470/1‑150 (MQ=255)
                                                                                                           tCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCTAAGGGAGTCATGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGAcc                                         >  1:1142477/1‑150 (MQ=255)
                                                                                                           tCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGAcc                                         <  1:2064642/150‑1 (MQ=255)
                                                                                                             atgatgGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGc                                                                      <  2:2277443/119‑1 (MQ=255)
                                                                                                             atgatgGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGc                                                                      >  1:2277443/1‑119 (MQ=255)
                                                                                                                atgGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATc                                    <  2:774552/150‑1 (MQ=255)
                                                                                                                 tgGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCa                                   <  2:2296152/150‑1 (MQ=255)
                                                                                                                 tgGCCTTGTCGAGGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCa                                   <  2:1247081/150‑1 (MQ=255)
                                                                                                                   gCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGc                                 <  1:942636/150‑1 (MQ=255)
                                                                                                                           gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCt                                                          <  2:207929/117‑1 (MQ=255)
                                                                                                                           gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAg                              <  2:1198762/145‑1 (MQ=255)
                                                                                                                             tggtggTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATg                       <  1:2156154/150‑1 (MQ=255)
                                                                                                                              ggtggtCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGc                                                       <  2:1884456/117‑1 (MQ=255)
                                                                                                                                  gtCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGt                  <  1:1918551/150‑1 (MQ=255)
                                                                                                                                    cGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAg                <  1:1071993/150‑1 (MQ=255)
                                                                                                                                    acgggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGa                                               <  2:258185/116‑1 (MQ=255)
                                                                                                                                    acgggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGa                                               >  1:258185/4‑119 (MQ=255)
                                                                                                                                      tGGCGTCGATGATCGGCGCCCGTGGATTGTGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCTAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCa              >  1:442637/1‑150 (MQ=255)
                                                                                                                                      tGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCa              <  2:1176903/150‑1 (MQ=255)
                                                                                                                                      gggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCg                                                              <  1:1308058/101‑1 (MQ=255)
                                                                                                                                      gggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCg                                                              >  2:1308058/2‑102 (MQ=255)
                                                                                                                                             gatgatCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAg                                                                                         <  2:612201/68‑1 (MQ=255)
                                                                                                                                             gatgatCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAg                                                                                         >  1:612201/1‑68 (MQ=255)
                                                                                                                                             gatgatCGGCGCCCGTGGATTGGGTGAAGATGTGCTAGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATCGCGGCCGAGGTGGCGATTGTCTCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCAGTCGCTg       >  1:2349599/1‑150 (MQ=255)
                                                                                                                                                  tCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATgg                      <  1:1215285/130‑1 (MQ=255)
                                                                                                                                                  tCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATgg                      >  2:1215285/1‑130 (MQ=255)
                                                                                                                                                  tCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCAGTCGCTGAGCCg  >  2:1342840/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCAGTCGCTGAGCCG  >  NC_002947/4038803‑4039098

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