Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 630,783 G→A G297D (GGC→GAC)  cstA → pyruvate transporter CstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913630,7830GA95.5% 219.7 / ‑5.7 67G297D (GGC→GAC) cstApyruvate transporter CstA
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (26/38);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 7.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCT  >  NC_000913/630538‑631029
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
ggggCGTATTGGCGAAGTGGCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACAttt                                                                                                                                                                                                                                                     <  2:450506/251‑1 (MQ=255)
 gggCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACAtttt                                                                                                                                                                                                                                                    <  1:106752/251‑1 (MQ=255)
       aTTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGCTAGACATTTTGATTa                                                                                                                                                                                                                                               <  1:374689/250‑1 (MQ=255)
       aTTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGAGTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTAt                                                                                                                                                                                                                                              <  2:32475/251‑1 (MQ=255)
                  gtCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACg                                                                                                                                                                                                                                    <  1:393020/250‑1 (MQ=255)
                       tCATCGGTCTGGTATTACTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGGGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGGCTATCCTCGGTCTGGGGGGAGACAATTTTATTTTTCGCCCGACGCTggcc                                                                                                                                                                                                                              >  1:406483/1‑248 (MQ=255)
                             gTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGGCTCCCACTCTACCCTTCTGCAAATTCGGACTATCGTCGGTTTGGCGGGAAGCCTTTTTATTATGCGCCCGCCGCTGGACATGccc                                                                                                                                                                                                                        >  1:515180/1‑250 (MQ=255)
                              tCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGGTGCCGGGGTGGTTACTGCCGGGCCCGCGTGTCTACCTCTCTTCCTTTCCTAAAATAGGGACTATCGTTGGGCTTGCGGGAAGCATTTTTTTTATGCGCCcgccgccggccctgcccc                                                                                                                                                                                                                       >  2:286074/1‑239 (MQ=255)
                                      tCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTTCTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGAc                                                                                                                                                                                                                <  1:187582/250‑1 (MQ=255)
                                        cTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCTGTCTGGCGGGAGACATTTTGATTATGCGCCCGACGCTGACCCTGCCTTCGCTGACCa                                                                                                                                                                                                              >  1:15358/1‑250 (MQ=255)
                                          gATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCaaa                                                                                                                                                                                                            <  2:143406/250‑1 (MQ=255)
                                          gATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAAt                                                                                                                                                                                                           <  1:240519/251‑1 (MQ=255)
                                                  ccATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGAt                                                                                                                                                                                                   <  1:141203/251‑1 (MQ=255)
                                                      tATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGc                                                                                                                                                                                                <  1:300483/250‑1 (MQ=255)
                                                                         gCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTTACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGCCGCTGACCCTGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGGATGGGCCgg                                                                                                                                                                           >  1:380367/1‑251 (MQ=255)
                                                                                     aCCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCCGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGAACTTTTGATTTTGCGCCCGGAGCTGACCCTGCCTGCGCCGACCAAAATTGTCGATGGCACTGGGCCGGTATggggccgtcacccgttccc                                                                                                                                                               >  2:572071/1‑234 (MQ=255)
                                                                                      ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCGGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCg                                                                                                                                                             <  1:330996/251‑1 (MQ=255)
                                                                                          ggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTcc                                                                                                                                                         <  2:99329/251‑1 (MQ=255)
                                                                                           gCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGACCTTTTTATTATTCGCCCGACGCTTACCCAGGCTGCGCTGACCAAATTTTTCGATGGCAATGGCCCGGTATGGGACGGGAAACTGTTCCCggttcc                                                                                                                                                        >  1:528079/1‑246 (MQ=255)
                                                                                            cACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTg                                                                                                                                                       <  1:484951/251‑1 (MQ=255)
                                                                                                           ttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGCGCTGCCGGCGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTCCCCTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGAAGCTGACCATGCCCGCGCTGACCAAATTTGTCGATGGCCCTGGCCCGGTGTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCaacccc                                                                                                                                        >  1:246498/1‑246 (MQ=255)
                                                                                                             aCCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCtg                                                                                                                                      <  1:332926/251‑1 (MQ=255)
                                                                                                                   gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCCGGGTGGTTACTGCTGGCCCCGCGGGACTCCCTCTCTACCCTCCTGCAAATCGGGACTATCGTCGGTCTGGCGGGTGGCATTTTTATTATGCGCCCGACGCTCACCATGCCCGCCCCGACCCAATTTGTCCATGGCACTGGGCCCGGTTGGGCCGGTAACCTGGTCCCCGTTccgttttttacccacccctgaggcgg                                                                                                                                 >  2:86012/1‑227 (MQ=255)
                                                                                                                     gCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGcgcg                                                                                                                               <  2:333318/250‑1 (MQ=255)
                                                                                                                      cAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTTCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTTTCGATGGCCATGGCCCGGGATGGACCGGGAACCTGTTCCCGTTCCTGTTTATTACCATCGCCTGTGGCGCggg                                                                                                                             >  1:448592/1‑250 (MQ=255)
                                                                                                                          tGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTc                                                                                                                          <  1:18839/250‑1 (MQ=255)
                                                                                                                                       ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                                                             <  1:489429/250‑1 (MQ=255)
                                                                                                                                          gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                                                                          <  1:143963/250‑1 (MQ=255)
                                                                                                                                          gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTGCTCGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGGCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCCGTGGCGCGGGGTCTGGCTTCCATgcggtgg                                                                                                         >  1:84914/1‑247 (MQ=255)
                                                                                                                                           ggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGACATTTTGATTATTCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGGAACCTGTTTCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                                                                        >  1:9005/1‑251 (MQ=255)
                                                                                                                                                         ggcggcGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCa                                                                                          <  2:435522/251‑1 (MQ=255)
                                                                                                                                                          gcggcgGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTAGCGTCGGTCTGGCGGTAGACAGTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAc                                                                                         <  2:473345/251‑1 (MQ=255)
                                                                                                                                                            ggcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGc                                                                                       >  1:382442/1‑251 (MQ=255)
                                                                                                                                                                     gcCGGTGTGTATACGGCTGGCCCCTCGTGACTTCCTCTCTACCTGCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTTCCCAAATTTGTCAATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCAGGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                                              <  2:124069/251‑1 (MQ=255)
                                                                                                                                                                       cGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                                                            >  1:109537/1‑251 (MQ=255)
                                                                                                                                                                                   tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATg                                                                                                                                                                                 <  1:152689/139‑1 (MQ=255)
                                                                                                                                                                                   tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATg                                                                                                                                                                                 >  2:152689/1‑139 (MQ=255)
                                                                                                                                                                                   tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAggg                                                                >  1:17655/1‑251 (MQ=255)
                                                                                                                                                                                    gctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAACATCGGGACTATCGTCGGTCTGGCGGTAGACCTTTTGATTATGCGCCCGACGCTGAACATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGGTCCTGTTTATCACCATCGCCTGTTGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGACGATGCTGGCGAACGAAggg                                                                >  2:71203/1‑250 (MQ=255)
                                                                                                                                                                                    gctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAgg                                                                 <  1:410289/249‑1 (MQ=255)
                                                                                                                                                                                       ggCCCCGCGTGACTACCTCTCTACCTTCCGGAAAAGCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCa                                                              <  1:252835/249‑1 (MQ=255)
                                                                                                                                                                                         ccccGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCg                                                          <  1:415257/251‑1 (MQ=255)
                                                                                                                                                                                                tGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTCGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTACTGTTTATCAACATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAAGCGTGCtttc                                                   >  2:346986/1‑250 (MQ=255)
                                                                                                                                                                                                  aCTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc                                                 <  2:402986/251‑1 (MQ=255)
                                                                                                                                                                                                     aCCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCgg                                               <  1:286074/250‑1 (MQ=255)
                                                                                                                                                                                                     aCCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCg                                                <  1:443263/249‑1 (MQ=255)
                                                                                                                                                                                                           ctACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCGTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACg                                          <  2:525987/249‑1 (MQ=255)
                                                                                                                                                                                                                 tCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGAt                                   <  1:548097/250‑1 (MQ=255)
                                                                                                                                                                                                                    cTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCCCTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCCTGCGCTGATCTCTTCTTGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTTCTTTATCGGCTACGGTGGGATGtt                                >  2:489403/1‑251 (MQ=255)
                                                                                                                                                                                                                      gAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGGGCTTTATCGGCTACGGGGGGATGTTa                               >  2:517413/1‑250 (MQ=255)
                                                                                                                                                                                                                      gAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCAGCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGtt                                <  1:47925/249‑1 (MQ=255)
                                                                                                                                                                                                                         aaTCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCATGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTTGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAAttg                           >  2:230097/1‑249 (MQ=255)
                                                                                                                                                                                                                             gggACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCTGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATc                       <  2:249002/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc                   <  2:372331/249‑1 (MQ=255)
                                                                                                                                                                                                                                     cgtcgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGGGGGATGTTAATGGCATCCTTCgggg               >  2:530776/1‑248 (MQ=255)
                                                                                                                                                                                                                                      gtcgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGc              <  1:448029/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tcgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCg             >  1:579307/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCg             >  1:418754/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGCGGTAGACATTTTGATTATGCGCCCGATGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAtt          <  2:349526/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAtt          >  1:67159/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAtt          >  1:300319/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGtt                                >  1:559318/1‑227 (MQ=255)
                                                                                                                                                                                                                                            cTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTAGCGGCTACGGTGGGATGtt                                <  2:559318/227‑1 (MQ=255)
                                                                                                                                                                                                                                             gggctgTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGTCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCTATTATg       <  2:373923/246‑1 (MQ=255)
                                                                                                                                                                                                                                                cGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGCAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGgcg    >  2:77173/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ggTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGgcgc   >  1:59372/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCt  <  2:381286/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
GGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCT  >  NC_000913/630538‑631029

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