Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,102,105 C→T A139V (GCT→GTT)  kdgT → 2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,102,1050CT97.9% 142.9 / NA 49A139V (GCT→GTT) kdgT2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/18);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGCTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCAATACCATTGATTTGACTGTGATT  >  NC_000913/4101873‑4102346
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
aTAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGAt                                                                                                                                                                                                                                 <  2:39033/251‑1 (MQ=255)
 tAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTAAGCTGGTGGTAACTAAAATCGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGAGCTGGTGGCGGCGATGGATATGAGCAATGGCGGACTTTACGCTTCCATCATGAAGCAGTACGGCACCAAAGAAGAAGTTGGGGCATATGTGTTGagg                                                                                                                                                                                                                                >  2:2700/1‑249 (MQ=255)
     aTTAAGCTCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGAGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGAAGTCAGGCACAAACGAAGACGCTGGGGCATTTGTGTTGATGTCGt                                                                                                                                                                                                                            >  2:145322/1‑251 (MQ=255)
        aaGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTgg                                                                                                                                                                                                                         >  2:191565/1‑251 (MQ=255)
        aaGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTgg                                                                                                                                                                                                                         >  1:332034/1‑251 (MQ=255)
           cgcgACGAGCACGGTACTGCGTCCATACGGTACGCTGGTGGTAACTACAATTGCCGTAGCGTGGGTGGTTGCGGCAATTGCCCCGCGTATTATTCCGGAACCTGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCCTGCAGCAGTACGGCACACAAGACGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGt                                                                                                                                                                                                                      >  2:269971/1‑251 (MQ=255)
                 gggAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAACGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGc                                                                                                                                                                                                                >  2:513324/1‑251 (MQ=255)
                                          aCGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCCATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGACGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTggg                                                                                                                                                                                       >  2:118933/1‑251 (MQ=255)
                                          aCGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTgg                                                                                                                                                                                        <  1:70924/250‑1 (MQ=255)
                                               ggtggtAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACt                                                                                                                                                                                   <  1:14566/250‑1 (MQ=255)
                                                    tAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGattatt                                                                                                                                                                                            >  1:623185/1‑236 (MQ=255)
                                                    tAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGattatt                                                                                                                                                                                            <  2:623185/236‑1 (MQ=255)
                                                     aaCTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCggg                                                                                                                                                                             >  1:315380/1‑250 (MQ=255)
                                                                ccGTCGCGTGGGTTGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGttt                                                                                                                                                                 <  1:275519/251‑1 (MQ=255)
                                                                    ggcgTGGGTGGTGGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAAc                                                                                                                                                             <  2:29003/250‑1 (MQ=255)
                                                                             ggtTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGttt                                                                                                                                                    >  2:366252/1‑251 (MQ=255)
                                                                             ggtTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATGCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTGGAACCGCATGttt                                                                                                                                                    <  1:82178/251‑1 (MQ=255)
                                                                             ggtTGCGGCAATGTCCTCGAGTATGATTCCGGAACATGGTGTTGAATGTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGTCGGCGTTGGATATGACCAATGGCGGACTGTAGGCGGCCATCATGCGGCAGGACGGCACAAAAAAATAAGTGGGGTCAATTGTGTTGTGGTCGTTGGAGTGCGTGACGCTCATGACGATGATTATTCTGGGCACTGCCGTGATTGCCTCGTGTGAACCGCATGtt                                                                                                                                                     <  2:282612/250‑1 (MQ=255)
                                                                                  cGGCACTTGCCTCGCGTATTATTCCGGAACATGGTGTTGCAGCTGTATTCTTTGCCGGACTTTCAACGATGGCGCTGGTGGCGGCGCTGGCTATGACCCATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGCGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCAGGGATCGACTCGTTTGAACCGCATGTTTtngtc                                                                                                                                               >  2:586921/1‑251 (MQ=255)
                                                                                  cGGCAATTGCCTCGCGTATGATGCCGGAACATGGTGTTGAAGTTGGATTCTTTGCTTGACTTTCAACTCTGGAGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGCTTGAGTCCGGGCCGCTCATGACGATGATTATGCTGGCCCCTGCCGGGAGTGCCTCGTTTGAGCCGCATGTTTtcgtc                                                                                                                                               <  2:360624/251‑1 (MQ=255)
                                                                                         tGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAATTTGGGGCATTTGTGTTGATGTCGTGGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCa                                                                                                                                         <  2:404556/250‑1 (MQ=255)
                                                                                         tGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGgcgc                                                                                                                                          >  1:553876/1‑249 (MQ=255)
                                                                                                     tattCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTcc                                                                                                                            >  1:59841/1‑251 (MQ=255)
                                                                                                        tCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGGCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCtg                                                                                                                          <  1:366252/250‑1 (MQ=255)
                                                                                                                  ggTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTg                                                                                                                <  1:213117/250‑1 (MQ=255)
                                                                                                                              ggATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCt                                                                                                   >  1:446547/1‑251 (MQ=255)
                                                                                                                                aTTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTg                                                                                                 >  2:414857/1‑251 (MQ=255)
                                                                                                                                      tGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTg                                                                                           >  2:553680/1‑251 (MQ=255)
                                                                                                                                         cGGACTTTCAACGCTGGCGCTGGTGGCGGAGATGGATATGACCGATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGGTGGAGTCCGAGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGTGAACCTTGACCCTGAAt                                                                                        <  2:159184/251‑1 (MQ=255)
                                                                                                                                             cTTTCCACGCTGGCGCTGGTGGCGGCGTTGGATATGACCCATTGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACCAACGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCCTGACTATGATTCTTCTTGGCACTGCCTGGCTTGCCTCGTTTTACCCGCATGTTTTCGTCCGCGCAGTCTTACCGTTCCTTCTGTGCTTTGCCCTGGTGACCCTTGACCCTGAAttccc                                                                                    >  2:50385/1‑248 (MQ=255)
                                                                                                                                                    cgctggcgctggTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAAttt                                                                             <  2:420759/251‑1 (MQ=255)
                                                                                                                                                     gctggcgctggTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAAtttt                                                                            >  1:436378/1‑251 (MQ=255)
                                                                                                                                                               gtggCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTTAGCAAAGc                                                                  >  2:388293/1‑251 (MQ=255)
                                                                                                                                                                             aatGACCAATGGCGGTCTTTACGCTTCCATCATGCATCAGTACGGCACAAAAGAAGAAGTTCGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACTATGATTATTCTGGGCACTGCCGGTATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGTTGGGCTTTGCCCTGGGGAACCTTGACCCTTAATTGCGAGAATTTTTCATCAAAGCGGTGCAAACGCtta                                                    <  2:612409/250‑3 (MQ=255)
                                                                                                                                                                                gACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTc                                                 <  1:89815/251‑1 (MQ=255)
                                                                                                                                                                                         cGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGAt                                                                                                                                                                           <  1:25663/120‑1 (MQ=255)
                                                                                                                                                                                         cGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGAt                                                                                                                                                                           >  2:25663/1‑120 (MQ=255)
                                                                                                                                                                                          ggACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGc                                       >  1:481719/1‑251 (MQ=255)
                                                                                                                                                                                             cTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCtt                                    >  2:56025/1‑251 (MQ=255)
                                                                                                                                                                                               ttACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGCTGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGTCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCtt                                    <  2:166286/249‑1 (MQ=255)
                                                                                                                                                                                                         atcatGCAGCAGTACGGCACAAAAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCa                         >  2:215631/1‑250 (MQ=255)
                                                                                                                                                                                                             tGCAGCAGTACGGCACAACAGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCAATAcc                    >  2:262806/1‑251 (MQ=255)
                                                                                                                                                                                                                 gcagTACGGCACAATAGAGGAAGTTGGGGCATTTGTGTTGATAACGATGGAGTCAGGGCCGCTCATGACGATGATTATTCTGGGCCCTGCCGGGATTGACTCGTTTTATCCGCATGTTTTCGTCGGCGCAGAATTACCGTTCCTGGTGGGCTTTGCCCTGGGGACACTGGACCCTGACTTGCGAGAATTTTTCAGCACAGCGGTGCACACGCTGATTACATTATTTGCCTTCGCGCTGGGCAATACCATTg                >  2:622323/1‑251 (MQ=255)
                                                                                                                                                                                                                      aCGGCACAAAAGAAGAAGGTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCAATACCATTGATTTg           <  2:332034/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggCACAAAAGAAGAGGTTGGGGCATTTGTGTTGATGTCGTTGGAGTTCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCAATACCATTGATTTGAc         >  2:464556/1‑251 (MQ=255)
                                                                                                                                                                                                                              aaaGAAGAAGTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCAATACCATTGATTTGACtgngct   >  2:608592/1‑249 (MQ=255)
                                                                                                                                                                                                                               aagaagaagTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCAATACCATTGATTTGACTGTGAtt  >  1:591397/1‑251 (MQ=255)
                                                                                                                                                                                                                                agaagaagTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGttt                                                                                                                                                    >  1:597632/1‑104 (MQ=255)
                                                                                                                                                                                                                                agaagaagTTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGttt                                                                                                                                                    <  2:597632/104‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
ATAAAATTAAGCGCGACGGGAACGGTACTGCGTAAATCCGGTACGCTGGTGGTAACTAAAATTGCCGTCGCGTGGGTGGTTGCGGCAATTGCCTCGCGTATTATTCCGGAACATGGTGTTGAAGTTGGATTCTTTGCCGGACTTTCAACGCTGGCGCTGGTGGCGGCGATGGATATGACCAATGGCGGACTTTACGCTTCCATCATGCAGCAGTACGGCACAAAAGAAGAAGCTGGGGCATTTGTGTTGATGTCGTTGGAGTCCGGGCCGCTCATGACGATGATTATTCTGGGCACTGCCGGGATTGCCTCGTTTGAACCGCATGTTTTCGTCGGCGCAGTATTACCGTTCCTGGTGGGCTTTGCCCTGGGGAACCTTGACCCTGAATTGCGAGAATTTTTCAGCAAAGCGGTGCAAACGCTGATTCCATTCTTTGCCTTCGCGCTGGGCAATACCATTGATTTGACTGTGATT  >  NC_000913/4101873‑4102346

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