Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,038,947 G→C 5.3% M234I (ATG→ATC 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,9470GC 5.3% 426.9 / 11.8 133M234I (ATG→ATCPP_3559glycine betaine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (56/70);  new base C (3/4);  total (59/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

Only 100 of 140 total aligned reads displayed.

CTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCAGTCGCTGAGC  >  NC_002947/4038800‑4039096
                                                                                                                                                   |                                                                                                                                                     
cTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCgatgat                                                                                                                                                     >  2:1154444/1‑150 (MQ=255)
cTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCgatgat                                                                                                                                                     >  1:1577992/1‑150 (MQ=255)
 tCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGGCCGGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATc                                                                                                                                                    <  1:1233533/150‑1 (MQ=255)
   gTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTCTCGATGGTGGTCGTGGCGTCGAGGATCgg                                                                                                                                                  >  1:947157/1‑150 (MQ=255)
        tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                             >  2:498159/1‑150 (MQ=255)
        tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                             >  1:1762982/1‑150 (MQ=255)
           aGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTg                                                                                                                                          >  2:794977/1‑150 (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)
                  acacGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTggg                                                                                                                                   >  2:2181265/1‑150 (MQ=255)
                     cGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGa                                                                                                                                >  2:2064642/1‑150 (MQ=255)
                      ggCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTgga                                                                                                                               >  1:2335211/1‑148 (MQ=255)
                        ccTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGa                                                                                                                             >  1:381746/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                                                                                                                                             >  2:1542686/1‑126 (MQ=255)
                                   gCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGCAGATGTGCTGGCgg                                                                                                                  >  2:224667/1‑150 (MQ=255)
                                   gCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCgg                                                                                                                  <  1:866328/150‑1 (MQ=255)
                                       cTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAAt                                                                                                              <  1:1045989/150‑1 (MQ=255)
                                           cAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGAAGAGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAg                                                                                                          >  2:942636/1‑150 (MQ=255)
                                                    cggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGGAGATGTGCTGGCGGGAATCCAGACCTTGAAt                                                                                                 <  2:2335211/150‑1 (MQ=255)
                                                    cggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAAt                                                                                                 <  1:1365834/150‑1 (MQ=255)
                                                     ggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCt                                                                                                                       >  2:119058/1‑127 (MQ=255)
                                                        gcATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGttc                                                                                             <  2:1653576/150‑2 (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                                                                                                            >  2:241065/1‑129 (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)
                                                                   ccccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATg                                                                                  >  1:959953/1‑150 (MQ=255)
                                                                    cccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATgg                                                                                 <  1:1728965/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)
                                                                        ggCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGc                                                                             <  1:2027786/150‑1 (MQ=255)
                                                                         gcgcTACCAAGCATCATGGCGGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGcc                                                                            <  2:658828/150‑1 (MQ=255)
                                                                           gcTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCgg                                                                          >  1:157317/1‑150 (MQ=255)
                                                                           gcTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCgg                                                                          <  2:2381762/150‑1 (MQ=255)
                                                                           gcTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGAAAGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCgg                                                                          >  1:637179/1‑150 (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                                                                                                                    >  1:1003585/1‑102 (MQ=255)
                                                                                 aaGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGc                                                                                                                    <  2:1003585/102‑1 (MQ=255)
                                                                                 aaGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGAGGGAAGCGGGGCTgtc                                                                    >  2:625492/1‑148 (MQ=255)
                                                                                   gCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCaa                                                                  >  2:1918551/1‑150 (MQ=255)
                                                                                     atcatGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAAtt                                                                <  1:2181265/150‑1 (MQ=255)
                                                                                     atcatGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAAtc                                                                <  1:224667/150‑2 (MQ=255)
                                                                                        atGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATg                                                                                                <  2:1085816/115‑1 (MQ=255)
                                                                                        atGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTc                                                             <  2:1492714/150‑1 (MQ=255)
                                                                                              gggCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTg                                                       <  1:2378237/150‑1 (MQ=255)
                                                                                               ggCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTgg                                                      <  2:1455307/150‑1 (MQ=255)
                                                                                                gCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGc                                                     >  1:553583/1‑150 (MQ=255)
                                                                                                gCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGc                                                     >  2:1369970/1‑150 (MQ=255)
                                                                                                  tCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCa                                                   >  2:989103/1‑150 (MQ=255)
                                                                                                     aCCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTgg                                                                     <  1:2041120/129‑1 (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)
                                                                                                        aGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGa                                             <  2:1861662/150‑1 (MQ=255)
                                                                                                          gcGGGCATGACGCCCTGGGCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATc                                           <  2:1715590/149‑1 (MQ=255)
                                                                                                            ggTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTgg                                                      >  1:996332/1‑137 (MQ=255)
                                                                                                            ggTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGa                                         >  2:1298455/1‑150 (MQ=255)
                                                                                                             gTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGAc                                        >  1:811470/1‑150 (MQ=255)
                                                                                                             gTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGATGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGAc                                        <  2:574340/150‑1 (MQ=255)
                                                                                                              tCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCTAAGGGAGTCATGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGAcc                                       >  1:1142477/1‑150 (MQ=255)
                                                                                                              tCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGAcc                                       <  1:2064642/150‑1 (MQ=255)
                                                                                                                atgatgGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGc                                                                    >  1:2277443/1‑119 (MQ=255)
                                                                                                                   atgGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATc                                  <  2:774552/150‑1 (MQ=255)
                                                                                                                    tgGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCa                                 <  2:2296152/150‑1 (MQ=255)
                                                                                                                      gCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGc                               <  1:942636/150‑1 (MQ=255)
                                                                                                                              gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCt                                                        >  1:207929/1‑117 (MQ=255)
                                                                                                                              gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATg                                                 >  2:1998186/1‑124 (MQ=255)
                                                                                                                              gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATg                                                 <  1:1998186/124‑1 (MQ=255)
                                                                                                                              gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAg                            <  2:1198762/145‑1 (MQ=255)
                                                                                                                              gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAg                            >  1:1198762/1‑145 (MQ=255)
                                                                                                                                tggtggTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATg                     <  1:2156154/150‑1 (MQ=255)
                                                                                                                                 ggtggtCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGc                                                     <  2:1884456/117‑1 (MQ=255)
                                                                                                                                 ggtggtCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGc                                                     >  1:1884456/1‑117 (MQ=255)
                                                                                                                                     gtCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGt                <  1:1918551/150‑1 (MQ=255)
                                                                                                                                      tCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTa               >  2:1071993/1‑150 (MQ=255)
                                                                                                                                       cGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAg              <  1:1071993/150‑1 (MQ=255)
                                                                                                                                       acgggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGa                                             <  2:258185/116‑1 (MQ=255)
                                                                                                                                       acgggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGa                                             >  1:258185/4‑119 (MQ=255)
                                                                                                                                        gTGGCGTCGATGATCGGCGCCCGTGGATTGGGCGAAGATGTGCTGGCGGGAATCAAGACCTTCAATGTTGGACAAGGGATGGAGGCAGGGCTGGCAATTGTCGCCCTGGCCATGGTGAACGACCGAATCAGCCAGGGCTATGGGAGTaac             >  2:20597/1‑148 (MQ=255)
                                                                                                                                         tGGCGTCGATGATCGGCGCCCGTGGATTGTGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCTAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCa            >  1:442637/1‑150 (MQ=255)
                                                                                                                                         gggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCg                                                            <  1:1308058/101‑1 (MQ=255)
                                                                                                                                         gggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCg                                                            >  2:1308058/2‑102 (MQ=255)
                                                                                                                                                gatgatCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAg                                                                                       >  1:612201/1‑68 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCAGTCGCTGAGC  >  NC_002947/4038800‑4039096

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: