Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,133,252 C→T P213S (CCG→TCG)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,133,2520CT96.3% 143.2 / NA 54P213S (CCG→TCG) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/34);  minor base G (1/1);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ATGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGCCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGA  >  NC_000913/4133007‑4133494
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
aTGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATGCCGGGAATTTTGTCGGt                                                                                                                                                                                                                                                    <  1:562484/250‑1 (MQ=255)
     tCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTa                                                                                                                                                                                                                                               <  1:148250/250‑1 (MQ=255)
         aCCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTGAT‑GTGG‑GAATTATTTCGGGGATTTTGTCGGTATCTAACttt                                                                                                                                                                                                                                          >  1:278037/1‑251 (MQ=255)
            gggTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCGCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTTCCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑aaa                                                                                                                                                                                                                                      <  2:79171/250‑1 (MQ=255)
                 acccTGTGAAAAGAAGTGTCAGATTACGATTTTTCCGGGGCGGCGGAACCGGAAGTTCACCCGTAAGCAAAAAGCGCTCAGGCGAAATGGCTTAATCTGAAACTCAAAGGGGATGCCGGAGCCAACCGCGCGATTCCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTGTGTCTTTATCTAACTTT‑AAACatgc                                                                                                                                                                                                                                 <  2:27094/248‑4 (MQ=255)
                 aCGCTGTTAAAAGAAGTGGCAGAGTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACGTT‑AAACAGG‑c                                                                                                                                                                                                                                <  2:258128/251‑1 (MQ=255)
                   gCTGTTAAAAGAAGTGGCAGATTTCGATAACTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTAGTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑Cga                                                                                                                                                                                                                              <  2:189647/251‑1 (MQ=255)
                      ggtAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑Cgaag                                                                                                                                                                                                                            <  1:544108/248‑1 (MQ=255)
                        tAAAAGAAGTGGCAGAGTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCGGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑Cgaagaa                                                                                                                                                                                                                          <  2:559863/250‑1 (MQ=255)
                                   gCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTTGT‑GTGG‑GAATTTTTCCGGGGATTTTGTTGGTATCTAACTTT‑TAACAGA‑CGAAGAAAATTTGccatat                                                                                                                                                                                                              >  1:398894/1‑247 (MQ=255)
                                      gATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAc                                                                                                                                                                                                          <  1:74015/251‑1 (MQ=255)
                                           cGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAAc                                                                                                                                                                                                     <  2:182679/250‑1 (MQ=255)
                                            gATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCGAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGt                                                                                                                                                                                                   <  1:99489/251‑1 (MQ=255)
                                                       ggcggcGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCGCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGGCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCgg                                                                                                                                                                                        <  2:373586/251‑1 (MQ=255)
                                                       ggcggaGTATACGGAAGTGCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTGGCCGATATGAC‑Caatgggggtatttcgg                                                                                                                                                                                        <  2:324195/251‑15 (MQ=255)
                                                                   ggAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATc                                                                                                                                                                                                                                                 <  1:587728/186‑1 (MQ=255)
                                                                      aGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTAAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGTTGTCGAAAGCTACCTGAGATTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑GAATTTTTACGGGAATTTTGTCGGTATATAACTTTCAAAC‑GG‑CGGAGAAA‑TTTTCTGATATGACACAA‑GTGCGTATTCCGGCGTgggaggggcca                                                                                                                                                                          >  1:419457/1‑241 (MQ=255)
                                                                        ttCACCCGGAAGCAAAAAGCGCTCAGGCGGCTTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGTT‑GTGG‑AAATTTTTTCGGGGATTTTTTCGGGATCTCACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGCTATGAA‑CAACGTGCGTATTTCGGCGGGGGTGGGGCaactg                                                                                                                                                                       >  1:543577/1‑248 (MQ=255)
                                                                        ttCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGACACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑CAATTTTTTCGGGGGTTTTGTCTGGATCTAACTTT‑TTACAGAAGGAAAAAA‑TTT‑ACGATATGAA‑CAACGTGCGTTTTCCCTCGAGGATGGCGCCAATg                                                                                                                                                                       >  2:273433/1‑251 (MQ=255)
                                                                              cGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTTTTTCGGGGATTTTTTCGGTATCTATCTTT‑TAACAGG‑CGAAGAAA‑TTTTCCGATATTAC‑CCAAGTGCGTATTCCGGCGTGGCTGGGGCGAATATTCGAc                                                                                                                                                                 >  1:349978/1‑251 (MQ=255)
                                                                                     aaaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGTT‑GTGG‑GAATTTTTTCGGGGATTTTGTCGGTGTCTAAATTT‑TAACAGG‑CGAAGAAA‑ATTGCCGATATGAC‑CAACATGCGTATTCCGCTGTGGCGGGGGCAAATGTTctgcggggtgg                                                                                                                                                          >  2:48081/1‑241 (MQ=255)
                                                                                      aaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCACAGTGGATGCCGGAGCCAACCGCGCTATTACTCAGTTCTTCTTCGAGGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCTGCGTGGATGGCGCAAATGTTCGACGGTCTGga                                                                                                                                                         <  1:520628/251‑1 (MQ=255)
                                                                                        aaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑GAATTTTTTCGGGGATTTTTTCGGTTTCTAACTTT‑TAACAAG‑CGGAGAAA‑ATTTCCGATATGAC‑CCACGTGCGTATTTCGGCGTGGATGGGGCAAATGTTCGGCGGTCTgggtg                                                                                                                                                       >  1:118001/1‑248 (MQ=255)
                                                                                                aGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGcc                                                                                                                                               <  1:67107/251‑1 (MQ=255)
                                                                                                            ttAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGAc                                                                                                                                                                 <  2:173556/221‑1 (MQ=255)
                                                                                                            ttAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGGATCTAACTTT‑AAAAAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTTCGGCGGGGATGGCGCAAATGTTCGAc                                                                                                                                                                 >  1:173556/1‑221 (MQ=255)
                                                                                                             tAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGGTTTGTAGCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATGTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGAGATGTC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCa                                                                                                                                     <  1:568628/248‑1 (MQ=255)
                                                                                                                          aaaGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCGTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATGATTCCGGGAATTTTGTCGGTATCTAACGTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTTTATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAATCTGGTGTGCGCGaa                                                                                                                     <  2:236570/251‑1 (MQ=255)
                                                                                                                              tGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGGATCTAACTTT‑AAACAAG‑CGAAGAAA‑TTTTCCGATATGAC‑CAACGTGCGGATTTCGGCGTGGATGGCGCAAATGTTCGACGGGCTGGGTGATGATGCCGAAACCCGCAAACAGGTTGGCGCGGATAtt                                                                                                                 >  1:18511/1‑251 (MQ=255)
                                                                                                                                      ggAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCAGTGAT‑GTGG‑GAATTATTTCGGGGATTTTGTCGGTATCTAACTTT‑TAACAGG‑CGAAGAAA‑TTTGTCGATATGAC‑AAACCTGCGTATTCCGGCGGGGGTGGCGCCAATGTTCGACGGTCTCGGTGATGGTGCCGAAACACCCAAAATGGTTGGCGCGGATCTTGcacagg                                                                                                          >  1:355376/1‑245 (MQ=255)
                                                                                                                                            aacggcgcgATTACTCATTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGGCGGTCTGGATGATGATGCCGAAACCCGCAAACTGGGTGTCGCGAAGATTGCCATGGATATGGt                                                                                                   <  2:52479/247‑1 (MQ=255)
                                                                                                                                                           cAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGc                                                                                      <  1:244140/249‑1 (MQ=255)
                                                                                                                                                             gTTCTTCGTCGATGTCGAAAGCTACCTGCGTTGTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTAGTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTg                                                                                  <  2:93834/251‑1 (MQ=255)
                                                                                                                                                                  tcttcGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTAGCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGaa                                                                                <  1:356388/248‑1 (MQ=255)
                                                                                                                                                                     tcGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAg                                                                            <  1:169671/249‑1 (MQ=255)
                                                                                                                                                                        aTGTCGAAGGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGGT‑GGGGGAATTTTTTCCGGGA‑CTTTGTTGGGATCTTAATTT‑TAACAAG‑CGGAGGAA‑TTTTACCATATGAA‑CAAAGTGGGTATTTCTGCGGGGGTGGGGCAAAAATTCGCCGCTCTTGATGATGGTGCCGAAAACCCCCAACTGGTTGGCGCGAAGAATGCCATGGCTTTGGTGGAGCATTTAAATCGTGAAAGAGTgga                                                                        >  2:544560/1‑248 (MQ=255)
                                                                                                                                                                        aTGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGGTTGTTTGGGG‑AAAT‑ATTCCGGGAATTTTTTTGGGTTCTCACTTT‑AAACAGG‑CGCAGAGA‑TTTGCCGATATGAT‑AAAGGTGCGTATTTCGGCGGGGCAGGCGCAAAAATTTGGCGGTCTGGGGTATGATACAGAAAACCGCAAACTGGTTGGCGCGAATATTGCCATTGATATGGTGGAAATTTTAAACCggggaagagggac                                                                        >  1:299705/1‑238 (MQ=255)
                                                                                                                                                                          gTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGAGATGGTAAAGATTTTAAGCCTTCAAGGATTGAAAGa                                                                     <  2:40971/251‑1 (MQ=255)
                                                                                                                                                                                 gCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCAc                                                              <  2:114733/251‑1 (MQ=255)
                                                                                                                                                                                    aCCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACt                                                             <  2:278037/249‑1 (MQ=255)
                                                                                                                                                                                          gTTTTCGTGACCGCTGTGTATCGGCGGGCATTTGT‑GTGG‑AAATTATTTCGGGGATTTTTTCGGGTTCTAAATTT‑TAAAAAG‑CGAAGAAA‑TTTTGCGATATGAC‑CCACATGCGTATTTCGGCGGGGGTGGGGCAAATTTTCGACGGGCTGGGTGATGATTCCGAAAACCGCAAACCGGTTGGGGGGGATATTGCCATGGGTATGGGGAAGATTTTTAATCGTGAAGGAGGGAAAGAATTACACTTTTATACg                                                     >  2:322258/1‑251 (MQ=255)
                                                                                                                                                                                                               cggcggGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTGT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAg                                 <  2:536512/250‑1 (MQ=255)
                                                                                                                                                                                                                  cggGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAc                             >  1:172666/1‑251 (MQ=255)
                                                                                                                                                                                                                    ggCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAcgc                           <  2:118001/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAcgc                           <  2:453606/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAcgc                           >  1:41757/1‑251 (MQ=255)
                                                                                                                                                                                                                              gtgG‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGcc                   <  1:445828/250‑1 (MQ=255)
                                                                                                                                                                                                                               tgG‑AAATTATTCCGGGAATTTTGGAGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGGGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAAACGTGCTGAAATGAGTTACGCGATTTGCCAt                 >  2:307642/1‑251 (MQ=255)
                                                                                                                                                                                                                                 g‑AAATTATTCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATAc               <  2:349978/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttattCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTc                                                                                                                                                                    <  2:119954/98‑1 (MQ=255)
                                                                                                                                                                                                                                      ttattCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTc                                                                                                                                                                    >  1:119954/1‑98 (MQ=255)
                                                                                                                                                                                                                                          tCCGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTggggg       >  2:370911/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccGGGAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGt      <  1:20368/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gAATTTTGTCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTGGCCATACGCTGGGGGTTCGa  <  1:203013/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
ATGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAT‑GTGG‑AAATTATTCCGGGAATTTTGCCGGTATCTAACTTT‑AAACAGG‑CGAAGAAA‑TTTGCCGATATGAC‑CAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGA  >  NC_000913/4133007‑4133494

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: ‑