Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,929,270 C→A L348M (CTG→ATG)  sdaC → L‑serine:H(+) symporter SdaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,929,2700CA88.6% 181.8 / 8.7 70L348M (CTG→ATG) sdaCL‑serine:H(+) symporter SdaC
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base A (23/39);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

GTCTTCAGCTGTGTACTGAGCCTGACTCCGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTACTCTCTGTTC  >  NC_000913/2929024‑2929512
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
gTCTTCAGCTGTGTACTGAGCCTGACTCCGGCAGACCTGGCTGCGGCTATAGCGCAGAACATCTCAATTCTGTATATCCTGGCTAACCACTTTAACGCACCGGTTATAGAGTGGACGGCTACGATTATCGCGATTATCGCTATCACTAAATCCTTCCTCGGTCACTCCCTGGGCGCACGTGAAGGCCTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGa                                                                                                                                                                                                                                                 >  1:362693/1‑250 (MQ=255)
    tCAGCTGTGTACTGAGCCTGACTCCGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAAc                                                                                                                                                                                                                                               <  1:348759/248‑1 (MQ=255)
         tgtgTACTGAGCCTGACTCCGTCAGACCTGGCTGCGGCTAAAGAGCATAACATCTCGATTCTGTCTTACCTGGCTTATCATTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCTCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCa                                                                                                                                                                                                                                        <  1:337842/250‑1 (MQ=255)
          gtgtACTGAGCCTGACTCCGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACt                                                                                                                                                                                                                                      >  1:501163/1‑251 (MQ=255)
              aCTGAGCCTGACTCCGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTgcgc                                                                                                                                                                                                                                  >  2:603302/1‑251 (MQ=255)
                tGAGCCTGACTCCGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTg                                                                                                                                                                                                                                >  1:467916/1‑251 (MQ=255)
                   gCCTGACTCCGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTGCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCGACAAGAGGAACCGTAGCACTGCGCTGTTc                                                                                                                                                                                                                             <  1:198633/251‑1 (MQ=255)
                    ccTGACTCCGGCAGTCCTGTCTGCTGCTAAAGAGCAGAACATCTCGATTCTGTCTTTCCTGGCTAACCACTTTTACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGTTATCACCAAATCCTGCCTCGGTCACTACCTGGGCGCACGTGAAGTCTTCAACGGTATGGTGATTAAATCTCTGCTAGGTAAAGGTAAGTCTATCGAAATCAACAAGATGATCCTTATCACTGCGCTGTTCa                                                                                                                                                                                                                            <  2:357771/251‑1 (MQ=255)
                           ccGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGt                                                                                                                                                                                                                     <  2:340040/251‑1 (MQ=255)
                               cAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATGCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATGATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCTCTGCGCGTTTCATGCTGGTAACg                                                                                                                                                                                                                 <  2:629342/251‑1 (MQ=255)
                                aGACCTGGCTGCGGCTAAAGAGCAGAACCTCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACg                                                                                                                                                                                                                 >  2:123878/1‑250 (MQ=255)
                                    cgggCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCTGTTATCGCGTGGATGGCGCCGATTATCGCGATTATCGCTAGCACCAAATCCTTCCTCGGTCACGACCTGGGCGCACGTGAAGTCTTCAACGGTATGGTGATTAAATCTCTGCTTGGAAAAGGTAAGTCTATCGAAACCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTg                                                                                                                                                                                                            <  2:153153/249‑1 (MQ=255)
                                     gggCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTgg                                                                                                                                                                                                           <  2:179788/250‑1 (MQ=255)
                                      ggCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCATTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGa                                                                                                                                                                                                          <  1:238591/251‑1 (MQ=255)
                                                   gatcAGAACATCTCGATTCTGTCTTACCTGGCTAACCCCTTTAACGCACCGGTTAACTCGTGGATTGCTCCGATTATCGCGATTATCGCTATCACCAAAGCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGGAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCt                                                                                                                                                                                             <  1:610392/248‑1 (MQ=255)
                                                   gagCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCt                                                                                                                                                                                             <  1:490594/251‑1 (MQ=255)
                                                   gagCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCt                                                                                                                                                                                             <  1:524782/251‑1 (MQ=255)
                                                                          ttACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGAGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCCTCCTGGGTATg                                                                                                                                                                       >  2:302759/1‑250 (MQ=255)
                                                                           tACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTAt                                                                                                                                                                        <  1:159960/248‑1 (MQ=255)
                                                                              cTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTg                                                                                                                                                                   >  1:618255/1‑250 (MQ=255)
                                                                                  cTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAAcc                                                                                                                                                              >  1:321837/1‑251 (MQ=255)
                                                                                        aCTTTAACGCACCGGTTATCACGTGGATGGCTCCGATTATCGCGCTTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGc                                                                                                                                                        >  1:398553/1‑251 (MQ=255)
                                                                                             aaCGCACCGTTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTTCCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGTTATGATTGAAACCCTGGGCGGTcc                                                                                                                                                   <  2:167217/251‑1 (MQ=255)
                                                                                                 cACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGTCCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAatc                                                                                                                                               <  2:219282/251‑1 (MQ=255)
                                                                                                   ccGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAatcat                                                                                                                                             <  1:333225/251‑1 (MQ=255)
                                                                                                    cGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAatcatc                                                                                                                                            <  2:270866/251‑1 (MQ=255)
                                                                                                      gTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATcgc                                                                                                                                          <  1:514153/251‑1 (MQ=255)
                                                                                                        tATCGCGTGGATGGCTCCGCTTATCGCGATTATCGCTATCACCACATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGAATCACTGCGCTGTTCATGCTGGTAACGAGCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCAAATCATcgcg                                                                                                                                         >  1:581309/1‑250 (MQ=255)
                                                                                                           cgcgTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCgatg                                                                                                                                      <  1:347693/250‑1 (MQ=255)
                                                                                                                              atcgcgattatcgcTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGAt                                                                                                                  >  2:210924/1‑251 (MQ=255)
                                                                                                                              atcgcgattatcgcTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGAt                                                                                                                  >  2:210902/1‑251 (MQ=255)
                                                                                                                                  cgattatcgcTATCACCAAATCCTTCCTCGGTCACTTCCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGAATAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAGCCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTAc                                                                                                              <  2:539429/251‑1 (MQ=255)
                                                                                                                                       atcgcTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAAt                                                                                                         <  2:398553/251‑1 (MQ=255)
                                                                                                                                         cgcTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGAAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTc                                                                                                       >  2:350/1‑251 (MQ=255)
                                                                                                                                          gcTATCAATAAATCCTTCCTCGGTCAATCCCTGGGCGCACGCGAAGGCTTCAACGGTCTCGTGATTCAATCTCTGCGTGGTAACGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTATTCATGCTGGTACCGACCTGGATTGTTGCCACCCTGACCCCGAGCATCCTGGGCATGTTTGAAACCCTGGGCAGTCCATTCATCGCGATGATCCTGTTCCTGATGACGATGTCCGCAATTc                                                                                                       >  2:64567/1‑250 (MQ=255)
                                                                                                                                               cACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGaa                                                                                                   >  1:459647/1‑249 (MQ=255)
                                                                                                                                                 ccAAATCCTTCCGCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTa                                                                                               <  1:387809/251‑1 (MQ=255)
                                                                                                                                                    aaTCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGAGTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCg                                                                                            <  1:92321/251‑1 (MQ=255)
                                                                                                                                                     aTCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCgg                                                                                           >  1:233285/1‑251 (MQ=255)
                                                                                                                                                       ccttcctCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGc                                                                                          >  2:369277/1‑250 (MQ=255)
                                                                                                                                                       ccttcctCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGc                                                                                          <  1:221057/250‑1 (MQ=255)
                                                                                                                                                         ttcctCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGa                                                                                                                                                                                                                                                 <  2:506205/97‑1 (MQ=255)
                                                                                                                                                         ttcctCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGa                                                                                                                                                                                                                                                 >  1:506205/1‑97 (MQ=255)
                                                                                                                                                                      aCCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCTGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACa                                                                            <  1:85094/249‑1 (MQ=255)
                                                                                                                                                                          gggcgaaagTGAAGGCTGCAACTGGATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTAGGATTGAAACCCTGTGCTGTCCAATCATCGGGATGAGTCTGTGCCTGATGCCGATGTATGCAATGCATAAAGTACCGGCAATGCGTAAGTACAGCgg                                                                        <  2:59016/241‑1 (MQ=255)
                                                                                                                                                                          gggCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGTACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTc                                                                      >  1:441412/1‑251 (MQ=255)
                                                                                                                                                                          gggCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGACCCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCAGCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGTCAATGCGTAAGTACAGCGGt                                                                       <  2:64162/250‑1 (MQ=255)
                                                                                                                                                                                aCGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCa                                                                <  1:292508/251‑1 (MQ=255)
                                                                                                                                                                                  gTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGc                                                              <  1:623035/251‑1 (MQ=255)
                                                                                                                                                                                    gAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACGGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAATCCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGTGGTCACATCAGCaa                                                            <  1:234076/251‑1 (MQ=255)
                                                                                                                                                                                          ttCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTAtt                                                      >  1:198156/1‑251 (MQ=255)
                                                                                                                                                                                           tCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTAtt                                                      <  2:501163/250‑1 (MQ=255)
                                                                                                                                                                                            cAACGGTATGGTGATTAAATCGCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCAGGCTTGTAACGACCTGGATTGTTGCCACCCGGAACCCGAGCATCCTGGGTATGATTTAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATTTACGCAATTCATAAAGTACCGGCAATGCGGAAGTACAGCGGTCACATCAGCAACGTATTc                                                     <  2:581975/250‑1 (MQ=255)
                                                                                                                                                                                                ggTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGt                                                <  2:321837/251‑1 (MQ=255)
                                                                                                                                                                                                  tATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCg                                              >  2:547083/1‑251 (MQ=255)
                                                                                                                                                                                                     ggTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAATCCCTGGGCGGTCCAATCATCGCGATGAGCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCTTTGTCGTGa                                           <  2:313564/251‑1 (MQ=255)
                                                                                                                                                                                                      gttcttAAATCTCCGCGTGTTAAAGGTAAGTCTATCGAAATCATCAAGATGACCCGTACCACCGCGCTGTTCATGCTGGTAACGACCGGGATTGTTGCCACCCTGAACCCTAGCATCCTGGGTATGATTGAAACCCTGGGCTGTCCAATCATCGCGAGTATCCTGTACCTGATGCCGATGTACGCAATTCAGAAAATACCGGTAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGCCGTGAt                                          <  2:47061/247‑1 (MQ=255)
                                                                                                                                                                                                                  ctGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCGGTTCATGCTGGTAACGACCTGGAGTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGAt                                 <  2:340972/248‑1 (MQ=255)
                                                                                                                                                                                                                    gCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGc                              <  1:302759/249‑1 (MQ=255)
                                                                                                                                                                                                                           aaaGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCg                      <  1:603302/250‑1 (MQ=255)
                                                                                                                                                                                                                           aaaGGTAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGc                     <  1:123878/251‑1 (MQ=255)
                                                                                                                                                                                                                                tAAGTCTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAAtct                >  1:410094/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aGTCTATCGAAATCAACAAGCTGAACCGTATCGCTGCGCGGTTCATGCTGGGAACGACCGGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGGCGTGATGGGTCTGATTGCAATCTCCGCAAtcttt              <  2:270420/251‑2 (MQ=255)
                                                                                                                                                                                                                                     cTATCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTAc           >  1:297554/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tCGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTActct        >  2:194551/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cGAAATCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGCCCTGGATTGTTGCCACCCTGAACCCGCGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGAATTCGTTGTCGAGATGGGTCTGATTGCAATCTCCGCAATCTTCTActctc       >  2:516850/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aTCAACAAGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCATGGGTATGATTGAAACCATGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTACTCTCTGtt   >  1:497882/1‑251 (MQ=255)
                                                                                                                                                                                                                                               caacaaGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTACTCTCTGTTc  <  2:432368/250‑1 (MQ=255)
                                                                                                                                                                                                                                               caacaaGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTACTCTCTGTTc  >  1:454473/1‑250 (MQ=255)
                                                                                                                                                                                                                                               caacaaGATGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTACTCTCTGTTc  >  1:456252/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  
GTCTTCAGCTGTGTACTGAGCCTGACTCCGGCAGACCTGGCTGCGGCTAAAGAGCAGAACATCTCGATTCTGTCTTACCTGGCTAACCACTTTAACGCACCGGTTATCGCGTGGATGGCTCCGATTATCGCGATTATCGCTATCACCAAATCCTTCCTCGGTCACTACCTGGGCGCACGTGAAGGCTTCAACGGTATGGTGATTAAATCTCTGCGTGGTAAAGGTAAGTCTATCGAAATCAACAAGCTGAACCGTATCACTGCGCTGTTCATGCTGGTAACGACCTGGATTGTTGCCACCCTGAACCCGAGCATCCTGGGTATGATTGAAACCCTGGGCGGTCCAATCATCGCGATGATCCTGTTCCTGATGCCGATGTACGCAATTCAGAAAGTACCGGCAATGCGTAAGTACAGCGGTCACATCAGCAACGTATTCGTTGTCGTGATGGGTCTGATTGCAATCTCCGCAATCTTCTACTCTCTGTTC  >  NC_000913/2929024‑2929512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: