Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,079,215 T→A 7.1% E9V (GAA→GTA) ‡ sdaC ← serine:H+ symport permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,079,2150TA 7.1% 258.1 / 9.1 85E9V (GAA→GTA) ‡sdaCserine:H+ symport permease
Reads supporting (aligned to +/- strand):  ref base T (34/45);  new base A (4/2);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

CAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCT  >  NC_002947/4079074‑4079363
                                                                                                                                             |                                                                                                                                                    
cAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATTGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCtt                                                                                                                                              <  1:1132158/150‑1 (MQ=255)
 aGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTg                                                                                                                                             <  1:129930/150‑1 (MQ=255)
  gAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGAACAACGCtttg                                                                                                                                            <  2:632362/150‑3 (MQ=255)
   aaGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTggg                                                                                                                                           >  1:578729/1‑150 (MQ=255)
      ccgccAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGcc                                                                                                                                        >  1:609598/1‑150 (MQ=255)
           aaGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCtc                                                                                                                                   >  2:1417183/1‑150 (MQ=255)
             ggCCGGCGTTGAACGACATGAATGGTGTCCCGGCGCGGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCtcat                                                                                                                                 <  1:533472/150‑1 (MQ=255)
                     ttGATCGGCAGGAACAGTGTCCCGGCACCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCCGCTCATTCATTCCg                                                                                                                         <  1:1134151/150‑1 (MQ=255)
                      tGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCgg                                                                                                                        <  2:1245223/150‑1 (MQ=255)
                      tGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCgg                                                                                                                        <  1:1417183/150‑1 (MQ=255)
                      tGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCgg                                                                                                                        <  1:935719/150‑1 (MQ=255)
                        aTCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGt                                                                                                                      <  2:1380277/150‑1 (MQ=255)
                        aTCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGt                                                                                                                      >  2:2121955/1‑150 (MQ=255)
                            gCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAg                                                                                                                  <  1:1535361/150‑1 (MQ=255)
                             cAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGc                                                                                                                 <  1:2121955/150‑1 (MQ=255)
                              aGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCt                                                                                                                <  2:1510151/150‑1 (MQ=255)
                              aGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCt                                                                                                                >  1:1447944/1‑150 (MQ=255)
                               ggAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGGCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTc                                                                                                               >  1:662895/1‑150 (MQ=255)
                                gAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGTTGGCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTcc                                                                                                              >  1:1750784/1‑150 (MQ=255)
                                gAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTcc                                                                                                              >  2:1014381/1‑150 (MQ=255)
                                gAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTcc                                                                                                              >  2:132142/1‑150 (MQ=255)
                                 aaCAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCa                                                                                                             <  1:1404422/150‑1 (MQ=255)
                                    aGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTCGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTc                                                                                                          <  1:1161753/150‑1 (MQ=255)
                                      tgtCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCtt                                                                                                                                              <  1:1806127/112‑1 (MQ=255)
                                      tgtCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCtt                                                                                                                                              >  2:1806127/1‑112 (MQ=255)
                                         cccGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAg                                                                                                     >  1:1091524/1‑150 (MQ=255)
                                           cGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGAc                                                                                                   <  2:1627772/150‑1 (MQ=255)
                                            ggCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACt                                                                                                  <  1:2263763/150‑1 (MQ=255)
                                                  gATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACa                                                                                            <  2:609598/150‑1 (MQ=255)
                                                   aTGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAg                                                                                           >  1:346276/1‑150 (MQ=255)
                                                    tGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAgg                                                                                          <  2:346276/150‑1 (MQ=255)
                                                      gCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCCCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCt                                                                                        <  1:430003/150‑1 (MQ=255)
                                                        tgtgCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTgg                                                                                                                                            <  2:564100/96‑1 (MQ=255)
                                                        tgtgCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTgg                                                                                                                                            >  1:564100/1‑96 (MQ=255)
                                                            ccGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAAt                                                                                  >  1:889469/1‑150 (MQ=255)
                                                             cGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAgg                                                                                          <  2:590318/141‑1 (MQ=255)
                                                             cGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAgg                                                                                          >  1:590318/1‑141 (MQ=255)
                                                              gAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATgg                                                                                >  1:1290857/1‑150 (MQ=255)
                                                               aaCAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACTGGCTGAGAATGGc                                                                               <  2:600861/150‑1 (MQ=255)
                                                               aaCAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGGTGAGAATTGc                                                                               <  2:599898/150‑1 (MQ=255)
                                                               aaCAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGc                                                                               <  1:558756/150‑1 (MQ=255)
                                                               aaCAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGACGGGGCCAAGTGCCACCCGGCATGGTTCAACGCTTTGGTCCTGCTCCTTCCTTACGGGTTTAGTGTCGCTCTCAAGACTGCAACAGGCTGAGAATGGc                                                                               <  1:551428/150‑1 (MQ=255)
                                                                     cccAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATa                                                                         >  2:764434/1‑150 (MQ=255)
                                                                     cccAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATa                                                                         <  1:519009/150‑1 (MQ=255)
                                                                      ccAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAg                                                                        >  2:2044543/1‑150 (MQ=255)
                                                                       cAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGa                                                                       >  1:464318/1‑150 (MQ=255)
                                                                           aTCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGAc                                                                      <  1:1925419/147‑1 (MQ=255)
                                                                           aTCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGAc                                                                      >  2:1925419/1‑147 (MQ=255)
                                                                            tCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTg                                                                                                 <  2:315486/119‑1 (MQ=255)
                                                                            tCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTg                                                                                                 >  1:315486/1‑119 (MQ=255)
                                                                             ccAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTc                                                                 >  1:745591/1‑150 (MQ=255)
                                                                               aGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGc                                                               >  1:2402191/1‑150 (MQ=255)
                                                                                 gtggtgTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGccc                                                             <  2:1447944/150‑1 (MQ=255)
                                                                                   ggtgTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCaa                                                           <  2:578729/150‑1 (MQ=255)
                                                                                         gTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCaa                                                                                                                                                   <  1:117598/56‑1 (MQ=255)
                                                                                         gTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCaa                                                                                                                                                   >  2:117598/1‑56 (MQ=255)
                                                                                         gTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAAc                                                     >  1:1120637/1‑150 (MQ=255)
                                                                                         gTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAAc                                                     >  2:2200296/1‑150 (MQ=255)
                                                                                          tGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAAcc                                                    <  1:132142/150‑1 (MQ=255)
                                                                                            aCGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGc                                                                                                                                     >  1:1028087/1‑67 (MQ=255)
                                                                                            aCGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGc                                                                                                                                     <  2:1028087/67‑1 (MQ=255)
                                                                                                gcCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGaaaa                                              >  2:1630710/1‑150 (MQ=255)
                                                                                                gcCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGaaaa                                              >  1:1386262/1‑150 (MQ=255)
                                                                                                  ccAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGGTAAGCTCCACTCGCAAGACTGCAACTGACTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGaaagag                                            >  1:162310/1‑147 (MQ=255)
                                                                                                      gTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAgggg                                        <  1:2200296/150‑1 (MQ=255)
                                                                                                                   cggGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCAc                           <  2:1750784/150‑1 (MQ=255)
                                                                                                                    gggggTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGc                                                                               <  1:799821/97‑1 (MQ=255)
                                                                                                                    gggggTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGc                                                                               >  2:799821/1‑97 (MQ=255)
                                                                                                                        gTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGa                                   <  2:1760280/137‑1 (MQ=255)
                                                                                                                        gTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGa                                   >  1:1760280/1‑137 (MQ=255)
                                                                                                                          cGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAgg                    <  1:458606/150‑1 (MQ=255)
                                                                                                                           gATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGaaaa                                              >  1:1175373/1‑123 (MQ=255)
                                                                                                                           gATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGaaaa                                              <  2:1175373/123‑1 (MQ=255)
                                                                                                                           gATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGt                   <  1:1630710/150‑1 (MQ=255)
                                                                                                                           gATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGt                   <  1:1079964/150‑1 (MQ=255)
                                                                                                                            aTTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGtt                  <  2:1763786/150‑1 (MQ=255)
                                                                                                                            aTTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGtt                  <  2:1290857/150‑1 (MQ=255)
                                                                                                                             ttCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTg                 <  2:889469/150‑1 (MQ=255)
                                                                                                                               cTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCa               <  1:742720/150‑1 (MQ=255)
                                                                                                                                 cccaaGCGTTGAACAACGCTTTGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAAATCACCCGGAAAAGAGGGGCTCTATTGTGCACGGGAAGGTTGCACt             >  2:1599323/4‑150 (MQ=255)
                                                                                                                                  ccaaGCGTTGAACAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGATCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACtt            >  1:1042443/3‑150 (MQ=255)
                                                                                                                                   caaGCGTTGAACAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGcc                                                              >  2:1442098/2‑99 (MQ=255)
                                                                                                                                   caaGCGTTGAACAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGcc                                                              <  1:1442098/98‑1 (MQ=255)
                                                                                                                                   caaGCGTTGAACAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAGCTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTg           >  1:1145195/2‑150 (MQ=255)
                                                                                                                                     aGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTACGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGACAAGAGGGGAGTGATTGGGCACGGAAATGTTGCACTTGCg         >  2:1162311/1‑150 (MQ=255)
                                                                                                                                        gttgttCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAg      >  2:2397645/1‑150 (MQ=255)
                                                                                                                                            ttCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCt  >  1:411534/1‑150 (MQ=255)
                                                                                                                                            ttCAACGCTTGGGGCCTGCTCATTCATTCCCGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGCAATAGACGTTTCGCCCAACTCAACCGGACAAGAGGGGCGCGATTGTGCCAGGAAAGGTTGCCCTTGCGAGCCCCt  >  2:107460/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCT  >  NC_002947/4079074‑4079363

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: