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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,079,2140TA 7.1% 260.5 / 9.4 85E9D (GAA→GAT) ‡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 = 2.99e-01

AGCAGCGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCT  >  NC_002947/4079065‑4079363
                                                                                                                                                     |                                                                                                                                                     
agcagcGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCgttgt                                                                                                                                                       >  2:1079964/1‑150 (MQ=255)
         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:935719/150‑1 (MQ=255)
                               tGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCgg                                                                                                                        <  1:1417183/150‑1 (MQ=255)
                                 aTCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGt                                                                                                                      >  2:2121955/1‑150 (MQ=255)
                                 aTCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGt                                                                                                                      <  2:1380277/150‑1 (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:132142/1‑150 (MQ=255)
                                         gAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTcc                                                                                                              >  2:1014381/1‑150 (MQ=255)
                                          aaCAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCa                                                                                                             <  1:1404422/150‑1 (MQ=255)
                                             aGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTCGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTc                                                                                                          <  1:1161753/150‑1 (MQ=255)
                                               tgtCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCtt                                                                                                                                              >  2:1806127/1‑112 (MQ=255)
                                               tgtCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCtt                                                                                                                                              <  1:1806127/112‑1 (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                                                                                          >  1:590318/1‑141 (MQ=255)
                                                                      cGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAgg                                                                                          <  2:590318/141‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
AGCAGCGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCT  >  NC_002947/4079065‑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: