Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,038,942 T→C 5.2% S233P (TCG→CCG)  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,9420TC 5.2% 411.9 / 7.9 137S233P (TCG→CCG) PP_3559glycine betaine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (60/69);  new base C (3/4);  total (63/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 = 2.29e-01

Only 100 of 141 total aligned reads displayed.

TCGACCCCTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCA  >  NC_002947/4038793‑4039086
                                                                                                                                                     |                                                                                                                                                
tCGACCCCTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCTAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGt                                                                                                                                                  <  1:665732/150‑1 (MQ=255)
tCGACCCCTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGt                                                                                                                                                  <  1:1368978/150‑1 (MQ=255)
       cTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCgatgat                                                                                                                                           >  2:875979/1‑150 (MQ=255)
        tCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGGCCGGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATc                                                                                                                                          <  1:1233533/150‑1 (MQ=255)
        tCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATc                                                                                                                                          >  2:1365834/1‑150 (MQ=255)
          gTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTCTCGATGGTGGTCGTGGCGTCGAGGATCgg                                                                                                                                        >  1:947157/1‑150 (MQ=255)
               tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                   >  1:1715590/1‑150 (MQ=255)
               tgcAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGccc                                                                                                                                   >  2:498159/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)
                  aGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTg                                                                                                                                >  2:794977/1‑150 (MQ=255)
                    gCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGa                                                                                                                              >  1:762857/1‑150 (MQ=255)
                      cgtgtgcGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGt                                                                                                                                                  >  1:2227571/7‑128 (MQ=255)
                      cgtgtgcGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGt                                                                                                                                                  <  2:2227571/122‑1 (MQ=255)
                      cGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGAtt                                                                                                                            <  1:1724291/150‑1 (MQ=255)
                      cGCACACGGCCTTGGCGCAAGCCGCTGGCACCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGAtt                                                                                                                            <  1:1209796/150‑1 (MQ=255)
                           acGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTg                                                                                                                       >  1:1895979/1‑150 (MQ=255)
                             ggCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTgga                                                                                                                     >  1:2335211/1‑148 (MQ=255)
                             ggCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGaa                                                                                                                     >  1:2030766/1‑150 (MQ=255)
                               ccTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGa                                                                                                                   >  1:381746/1‑150 (MQ=255)
                                 ttGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGAtg                                                                                                                 <  1:875979/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)
                                      gcAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGCCGATCATCGACGcca                                                                                                                                   <  1:1542686/127‑2 (MQ=255)
                                       cAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTgg                                                                                                           <  2:1577992/150‑1 (MQ=255)
                                          gCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGCAGATGTGCTGGCgg                                                                                                        >  2:224667/1‑150 (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)
                                                           cggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGGAGATGTGCTGGCGGGAATCCAGACCTTGAAt                                                                                       <  2:2335211/150‑1 (MQ=255)
                                                           cggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAAt                                                                                       <  1:1365834/150‑1 (MQ=255)
                                                           cggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAAt                                                                                       <  1:1798379/150‑1 (MQ=255)
                                                            ggcgCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCt                                                                                                             >  2:119058/1‑127 (MQ=255)
                                                             gcgcATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGt                                                                                     >  2:2378237/1‑150 (MQ=255)
                                                               gcATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTg                                                                                   <  2:936107/150‑1 (MQ=255)
                                                                 aTCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGc                                                                                 >  1:566916/1‑150 (MQ=255)
                                                                     cccTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATcc                                                                                                  >  2:241065/1‑129 (MQ=255)
                                                                     cccTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAg                                                                             <  2:1639869/150‑1 (MQ=255)
                                                                      ccTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCCTGTCGAGGGTGGTCGTGGTGTCGATGAGCGGCGCCCGTGGATTTCGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAgg                                                                            <  1:1855536/150‑1 (MQ=255)
                                                                       cTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATTTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAggg                                                                           >  1:2016861/1‑150 (MQ=255)
                                                                       cTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCg                                                                                                         >  1:2001996/1‑120 (MQ=255)
                                                                       cTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCg                                                                                                         <  2:2001996/120‑1 (MQ=255)
                                                                           cccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATgg                                                                       <  1:1728965/150‑1 (MQ=255)
                                                                           cccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATgg                                                                       <  2:1037777/150‑1 (MQ=255)
                                                                            ccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTg                                                                                                                           <  1:617289/97‑1 (MQ=255)
                                                                            ccTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGa                                                                      <  2:959953/150‑1 (MQ=255)
                                                                                gcgcTACCAAGCATCATGGCGGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGcc                                                                  <  2:658828/150‑1 (MQ=255)
                                                                                gcgcTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTAGAGACTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGcc                                                                  <  1:24953/150‑1 (MQ=255)
                                                                                  gcTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCgg                                                                >  1:157317/1‑150 (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)
                                                                                       cAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCg                                                                 <  1:1320574/144‑1 (MQ=255)
                                                                                       cAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCg                                                                 >  2:1320574/1‑144 (MQ=255)
                                                                                        aaGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGc                                                                                                          >  1:1003585/1‑102 (MQ=255)
                                                                                        aaGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGAGGGAAGCGGGGCTgtc                                                          >  2:625492/1‑148 (MQ=255)
                                                                                          gCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCaa                                                        >  2:1918551/1‑150 (MQ=255)
                                                                                            atcatGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAAtt                                                      <  1:794977/150‑1 (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)
                                                                                                           aaCCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATg                                       <  1:989103/150‑1 (MQ=255)
                                                                                                            aCCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTgg                                                           <  1:2041120/129‑1 (MQ=255)
                                                                                                            aCCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTgg                                                           >  2:2041120/1‑129 (MQ=255)
                                                                                                             ccAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGt                                     >  1:585529/1‑150 (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)
                                                                                                                       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)
                                                                                                                                     gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCt                                              >  1:207929/1‑117 (MQ=255)
                                                                                                                                     gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATg                                       <  1:1998186/124‑1 (MQ=255)
                                                                                                                                     gATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAg                  >  1:1198762/1‑145 (MQ=255)
                                                                                                                                       tggtggTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATg           <  1:2156154/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)
                                                                                                                                                tGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCa  <  2:1176903/150‑1 (MQ=255)
                                                                                                                                                gggCGCCGATCATCGACGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCg                                                  <  1:1308058/101‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TCGACCCCTCGTTGCTGCAGGCCGCACACGGCCTTGGCGCAAGCCGCTGGCAGCGGTTGCGGCGCATCGCCCTGCCCCTGGCGCTACCAAGCATCATGGCCGGGCTCAACCAGTCGGTCATGATGGCCTTGTCGATGGTGGTCGTGGCGTCGATGATCGGCGCCCGTGGATTGGGTGAAGATGTGCTGGCGGGAATCCAGACCTTGAATGTTGGCCAAGGGATGGAGGCCGGGCTGGCAATTGTCGCCCTGGCCATGGTGATCGACCGGATCAGCCAGGCCTATGGGCGTAGCA  >  NC_002947/4038793‑4039086

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: