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

GGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCT  >  NC_000913/630538‑631029
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
ggggCGTATTGGCGAAGTGGCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACAttt                                                                                                                                                                                                                                                     <  2:450506/251‑1 (MQ=255)
 gggCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACAtttt                                                                                                                                                                                                                                                    <  1:106752/251‑1 (MQ=255)
       aTTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGCTAGACATTTTGATTa                                                                                                                                                                                                                                               <  1:374689/250‑1 (MQ=255)
       aTTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGAGTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTAt                                                                                                                                                                                                                                              <  2:32475/251‑1 (MQ=255)
                  gtCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACg                                                                                                                                                                                                                                    <  1:393020/250‑1 (MQ=255)
                       tCATCGGTCTGGTATTACTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGGGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGGCTATCCTCGGTCTGGGGGGAGACAATTTTATTTTTCGCCCGACGCTggcc                                                                                                                                                                                                                              >  1:406483/1‑248 (MQ=255)
                             gTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGGCTCCCACTCTACCCT‑TCTGCAAATTCGGACTATCGTCGGTTTGGCGGGAAGCCTTTTTATTATGCGCCCGCCGCTGGACATGccc                                                                                                                                                                                                                        >  1:515180/1‑250 (MQ=255)
                              tCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGGGGTGCCGGGGTGGTTACTGCCGGGCCCGCGTGTCTACCTCTCTTCCTTTCCT‑AAAATAGGGACTATCGTTGGGCTTGCGGGAAGCATTTTTTTTATGCGCCcgccgccggccctgcccc                                                                                                                                                                                                                       >  2:286074/1‑239 (MQ=255)
                                      tCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTTCTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGAc                                                                                                                                                                                                                <  1:187582/250‑1 (MQ=255)
                                        cTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCTGTCTGGCGGGAGACATTTTGATTATGCGCCCGACGCTGACCCTGCCTTCGCTGACCa                                                                                                                                                                                                              >  1:15358/1‑250 (MQ=255)
                                          gATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCaaa                                                                                                                                                                                                            <  2:143406/250‑1 (MQ=255)
                                          gATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAAt                                                                                                                                                                                                           <  1:240519/251‑1 (MQ=255)
                                                  ccATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGAt                                                                                                                                                                                                   <  1:141203/251‑1 (MQ=255)
                                                      tATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGc                                                                                                                                                                                                <  1:300483/250‑1 (MQ=255)
                                                                         gCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTTACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGCCGCTGACCCTGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGGATGGGCCgg                                                                                                                                                                           >  1:380367/1‑251 (MQ=255)
                                                                                     aCCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCCGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGAACTTTTGATTTTGCGCCCGGAGCTGACCCTGCCTGCGCCGACCAAAATTGTCGATGGCACTGG‑GCCGGTATggggccgtcacccgttccc                                                                                                                                                               >  2:572071/1‑234 (MQ=255)
                                                                                      ccTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCGGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCg                                                                                                                                                             <  1:330996/251‑1 (MQ=255)
                                                                                          ggCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTcc                                                                                                                                                         <  2:99329/251‑1 (MQ=255)
                                                                                           gCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGACCTTTTTATTATTCGCCCGACGCTTACCCAGGCTGCGCTGACCAAATTTTTCGATGGCAATGG‑CCCGGTATGGGACGGGAAACTG‑TTCCCggttcc                                                                                                                                                        >  1:528079/1‑246 (MQ=255)
                                                                                            cACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTg                                                                                                                                                       <  1:484951/251‑1 (MQ=255)
                                                                                                           ttACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGCGCTGCCGGCGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTCCCCT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGAAGCTGACCATGCCCGCGCTGACCAAATTTGTCGATGGCCCTGG‑CCCGGTGTGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCaacccc                                                                                                                                        >  1:246498/1‑246 (MQ=255)
                                                                                                             aCCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCtg                                                                                                                                      <  1:332926/251‑1 (MQ=255)
                                                                                                                   gTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCCGGGTGGTTACTGCTGGCCCCGCGGGACTCCCTCTCTACCCT‑CCTGCAAATCGGGACTATCGTCGGTCTGGCGGGTGGCATTTTTATTATGCGCCCGACGCTCACCATGCCCGCCCCGACCCAATTTGTCCATGGCACTGGGCCCGGT‑TGGGCCGGTAACCTGGTCCCCGTTccgttttttacccacccctgaggcgg                                                                                                                                 >  2:86012/1‑227 (MQ=255)
                                                                                                                     gCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGcgcg                                                                                                                               <  2:333318/250‑1 (MQ=255)
                                                                                                                      cAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTTCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTTTCGATGGCCATGG‑CCCGGGATGGACCGGGAACCTG‑TTCCCGTTCCTGTTTATTACCATCGCCTGTGGCGCggg                                                                                                                             >  1:448592/1‑250 (MQ=255)
                                                                                                                          tGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTc                                                                                                                          <  1:18839/250‑1 (MQ=255)
                                                                                                                                       ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                                                             <  1:489429/250‑1 (MQ=255)
                                                                                                                                          gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                                                                          <  1:143963/250‑1 (MQ=255)
                                                                                                                                          gggTTACGGTTTTGTGGCGGCGGTGCTGCCGGGGTGGTTACTGCTCGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑GCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCCGTGGCGCGGGGTCTGGCTTCCATgcggtgg                                                                                                         >  1:84914/1‑247 (MQ=255)
                                                                                                                                           ggTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGGAGACATTTTGATTATTCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGGAACCTG‑TTTCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                                                                        >  1:9005/1‑251 (MQ=255)
                                                                                                                                                         ggcggcGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCa                                                                                          <  2:435522/251‑1 (MQ=255)
                                                                                                                                                          gcggcgGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTAGCGTCGGTCTGGCGGTAGACAGTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAc                                                                                         <  2:473345/251‑1 (MQ=255)
                                                                                                                                                            ggcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGc                                                                                       >  1:382442/1‑251 (MQ=255)
                                                                                                                                                                     gcCGGTGTGTATACGGCTGGCCCCTCGTGACTTCCTCTCTACCTG‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTTCCCAAATTTGTCAATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCAGGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                                              <  2:124069/251‑1 (MQ=255)
                                                                                                                                                                       cGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                                                            >  1:109537/1‑251 (MQ=255)
                                                                                                                                                                                   tgctgGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATg                                                                                                                                                                                 <  1:152689/139‑1 (MQ=255)
                                                                                                                                                                                   tgctgGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATg                                                                                                                                                                                 >  2:152689/1‑139 (MQ=255)
                                                                                                                                                                                   tgctgGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAggg                                                                >  1:17655/1‑251 (MQ=255)
                                                                                                                                                                                    gctgGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAACATCGGGACTATCGTCGGTCTGGCGGTAGACCTTTTGATTATGCGCCCGACGCTGAACATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGGTCCTGTTTATCACCATCGCCTGTTGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGACGATGCTGGCGAACGAAggg                                                                >  2:71203/1‑250 (MQ=255)
                                                                                                                                                                                    gctgGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAgg                                                                 <  1:410289/249‑1 (MQ=255)
                                                                                                                                                                                       ggCCCCGCGTGACTACCTCTCTACCTT‑CCGGAAAAGCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCa                                                              <  1:252835/249‑1 (MQ=255)
                                                                                                                                                                                         ccccGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCg                                                          <  1:415257/251‑1 (MQ=255)
                                                                                                                                                                                                tGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTCGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTACTGTTTATCAACATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAAGCGTGCtttc                                                   >  2:346986/1‑250 (MQ=255)
                                                                                                                                                                                                  aCTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc                                                 <  2:402986/251‑1 (MQ=255)
                                                                                                                                                                                                     aCCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCgg                                               <  1:286074/250‑1 (MQ=255)
                                                                                                                                                                                                     aCCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCg                                                <  1:443263/249‑1 (MQ=255)
                                                                                                                                                                                                           ctACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCGTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACg                                          <  2:525987/249‑1 (MQ=255)
                                                                                                                                                                                                                 t‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGAt                                   <  1:548097/250‑1 (MQ=255)
                                                                                                                                                                                                                    cTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCCCTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCCTGCGCTGATCTCTTCTTGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTTCTTTATCGGCTACGGTGGGATGtt                                >  2:489403/1‑251 (MQ=255)
                                                                                                                                                                                                                      gAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGGGCTTTATCGGCTACGGGGGGATGTTa                               >  2:517413/1‑250 (MQ=255)
                                                                                                                                                                                                                      gAAAATCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCAGCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGtt                                <  1:47925/249‑1 (MQ=255)
                                                                                                                                                                                                                         aaTCGGGACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCATGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTTGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAAttg                           >  2:230097/1‑249 (MQ=255)
                                                                                                                                                                                                                             gggACTATCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCTGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATc                       <  2:249002/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTCGTCGGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc                   <  2:372331/249‑1 (MQ=255)
                                                                                                                                                                                                                                     cgtcgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGGGGGATGTTAATGGCATCCTTCgggg               >  2:530776/1‑248 (MQ=255)
                                                                                                                                                                                                                                      gtcgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGc              <  1:448029/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tcgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCg             >  1:579307/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cgGTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCg             >  1:418754/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGCGGTAGACATTTTGATTATGCGCCCGATGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAtt          <  2:349526/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAtt          >  1:67159/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tCTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGAtt          >  1:300319/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGtt                                >  1:559318/1‑227 (MQ=255)
                                                                                                                                                                                                                                            cTGGCGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTAGCGGCTACGGTGGGATGtt                                <  2:559318/227‑1 (MQ=255)
                                                                                                                                                                                                                                             gggctgTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGTCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCTATTATg       <  2:373923/246‑1 (MQ=255)
                                                                                                                                                                                                                                                cGGTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGCAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGgcg    >  2:77173/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 ggTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGgcgc   >  1:59372/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gTAGACATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCt  <  2:381286/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
GGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTTTCGCCATTATCTCTGGCGGCTGGGTGGCAGAAAGTCCGACCTGGGCACCGTACTTTGACTTTACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTT‑CCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGATTATGCGCCCGACGCTGACCATGCCTGCGCTGACCAAATTTGTCGATGGCACTGG‑CCCGGTATGGACCGGTAACCTG‑TTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCT  >  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: ‑