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

ATGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGA  >  NC_000913/4133007‑4133494
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
aTGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATGCCGGGAATTTTGTCGGt                                                                                                                                                                                                                                                    <  1:562484/250‑1 (MQ=255)
     tCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTa                                                                                                                                                                                                                                               <  1:148250/250‑1 (MQ=255)
         aCCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTGATGTGGGAATTATTTCGGGGATTTTGTCGGTATCTAACttt                                                                                                                                                                                                                                          >  1:278037/1‑251 (MQ=255)
            gggTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCGCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTTCCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTaaa                                                                                                                                                                                                                                      <  2:79171/250‑1 (MQ=255)
                 acccTGTGAAAAGAAGTGTCAGATTACGATTTTTCCGGGGCGGCGGAACCGGAAGTTCACCCGTAAGCAAAAAGCGCTCAGGCGAAATGGCTTAATCTGAAACTCAAAGGGGATGCCGGAGCCAACCGCGCGATTCCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTGTGTCTTTATCTAACTTTAAACatgc                                                                                                                                                                                                                                 <  2:27094/248‑4 (MQ=255)
                 aCGCTGTTAAAAGAAGTGGCAGAGTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACGTTAAACAGGc                                                                                                                                                                                                                                <  2:258128/251‑1 (MQ=255)
                   gCTGTTAAAAGAAGTGGCAGATTTCGATAACTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTAGTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCga                                                                                                                                                                                                                              <  2:189647/251‑1 (MQ=255)
                      ggtAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCgaag                                                                                                                                                                                                                            <  1:544108/248‑1 (MQ=255)
                        tAAAAGAAGTGGCAGAGTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCGGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCgaagaa                                                                                                                                                                                                                          <  2:559863/250‑1 (MQ=255)
                                   gCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTTGTGTGGGAATTTTTCCGGGGATTTTGTTGGTATCTAACTTTTAACAGACGAAGAAAATTTGccatat                                                                                                                                                                                                              >  1:398894/1‑247 (MQ=255)
                                      gATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGAc                                                                                                                                                                                                          <  1:74015/251‑1 (MQ=255)
                                           cGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAAc                                                                                                                                                                                                     <  2:182679/250‑1 (MQ=255)
                                            gATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCGAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGt                                                                                                                                                                                                   <  1:99489/251‑1 (MQ=255)
                                                       ggcggcGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCGCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGGCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCgg                                                                                                                                                                                        <  2:373586/251‑1 (MQ=255)
                                                       ggcggaGTATACGGAAGTGCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTGGCCGATATGACCaatgggggtatttcgg                                                                                                                                                                                        <  2:324195/251‑15 (MQ=255)
                                                                   ggAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATc                                                                                                                                                                                                                                                 <  1:587728/186‑1 (MQ=255)
                                                                      aGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTAAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGTTGTCGAAAGCTACCTGAGATTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGGAATTTTTACGGGAATTTTGTCGGTATATAACTTTCAAACGGCGGAGAAATTTTCTGATATGACACAAGTGCGTATTCCGGCGTgggaggggcca                                                                                                                                                                          >  1:419457/1‑241 (MQ=255)
                                                                        ttCACCCGGAAGCAAAAAGCGCTCAGGCGGCTTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGTTGTGGAAATTTTTTCGGGGATTTTTTCGGGATCTCACTTTAAACAGGCGAAGAAATTTGCCGCTATGAACAACGTGCGTATTTCGGCGGGGGTGGGGCaactg                                                                                                                                                                       >  1:543577/1‑248 (MQ=255)
                                                                        ttCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGACACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGCAATTTTTTCGGGGGTTTTGTCTGGATCTAACTTTTTACAGAAGGAAAAAATTTACGATATGAACAACGTGCGTTTTCCCTCGAGGATGGCGCCAATg                                                                                                                                                                       >  2:273433/1‑251 (MQ=255)
                                                                              cGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTTTTTCGGGGATTTTTTCGGTATCTATCTTTTAACAGGCGAAGAAATTTTCCGATATTACCCAAGTGCGTATTCCGGCGTGGCTGGGGCGAATATTCGAc                                                                                                                                                                 >  1:349978/1‑251 (MQ=255)
                                                                                     aaaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGTTGTGGGAATTTTTTCGGGGATTTTGTCGGTGTCTAAATTTTAACAGGCGAAGAAAATTGCCGATATGACCAACATGCGTATTCCGCTGTGGCGGGGGCAAATGTTctgcggggtgg                                                                                                                                                          >  2:48081/1‑241 (MQ=255)
                                                                                      aaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCACAGTGGATGCCGGAGCCAACCGCGCTATTACTCAGTTCTTCTTCGAGGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCTGCGTGGATGGCGCAAATGTTCGACGGTCTGga                                                                                                                                                         <  1:520628/251‑1 (MQ=255)
                                                                                        aaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGGAATTTTTTCGGGGATTTTTTCGGTTTCTAACTTTTAACAAGCGGAGAAAATTTCCGATATGACCCACGTGCGTATTTCGGCGTGGATGGGGCAAATGTTCGGCGGTCTgggtg                                                                                                                                                       >  1:118001/1‑248 (MQ=255)
                                                                                                aGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGcc                                                                                                                                               <  1:67107/251‑1 (MQ=255)
                                                                                                            ttAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGAc                                                                                                                                                                 <  2:173556/221‑1 (MQ=255)
                                                                                                            ttAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGGATCTAACTTTAAAAAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTTCGGCGGGGATGGCGCAAATGTTCGAc                                                                                                                                                                 >  1:173556/1‑221 (MQ=255)
                                                                                                             tAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGGTTTGTAGCGGCGGGCATTGATGTGGAAATTATTCCGGGAATGTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGAGATGTCCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCa                                                                                                                                     <  1:568628/248‑1 (MQ=255)
                                                                                                                          aaaGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCGTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATGATTCCGGGAATTTTGTCGGTATCTAACGTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTTTATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAATCTGGTGTGCGCGaa                                                                                                                     <  2:236570/251‑1 (MQ=255)
                                                                                                                              tGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGGATCTAACTTTAAACAAGCGAAGAAATTTTCCGATATGACCAACGTGCGGATTTCGGCGTGGATGGCGCAAATGTTCGACGGGCTGGGTGATGATGCCGAAACCCGCAAACAGGTTGGCGCGGATAtt                                                                                                                 >  1:18511/1‑251 (MQ=255)
                                                                                                                                      ggAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCAGTGATGTGGGAATTATTTCGGGGATTTTGTCGGTATCTAACTTTTAACAGGCGAAGAAATTTGTCGATATGACAAACCTGCGTATTCCGGCGGGGGTGGCGCCAATGTTCGACGGTCTCGGTGATGGTGCCGAAACACCCAAAATGGTTGGCGCGGATCTTGcacagg                                                                                                          >  1:355376/1‑245 (MQ=255)
                                                                                                                                            aacggcgcgATTACTCATTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGGCGGTCTGGATGATGATGCCGAAACCCGCAAACTGGGTGTCGCGAAGATTGCCATGGATATGGt                                                                                                   <  2:52479/247‑1 (MQ=255)
                                                                                                                                                           cAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGc                                                                                      <  1:244140/249‑1 (MQ=255)
                                                                                                                                                             gTTCTTCGTCGATGTCGAAAGCTACCTGCGTTGTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTAGTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTg                                                                                  <  2:93834/251‑1 (MQ=255)
                                                                                                                                                                  tcttcGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTAGCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGaa                                                                                <  1:356388/248‑1 (MQ=255)
                                                                                                                                                                     tcGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAg                                                                            <  1:169671/249‑1 (MQ=255)
                                                                                                                                                                        aTGTCGAAGGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGGTGGGGGAATTTTTTCCGGGACTTTGTTGGGATCTTAATTTTAACAAGCGGAGGAATTTTACCATATGAACAAAGTGGGTATTTCTGCGGGGGTGGGGCAAAAATTCGCCGCTCTTGATGATGGTGCCGAAAACCCCCAACTGGTTGGCGCGAAGAATGCCATGGCTTTGGTGGAGCATTTAAATCGTGAAAGAGTgga                                                                        >  2:544560/1‑248 (MQ=255)
                                                                                                                                                                        aTGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGGTTGTTTGGGGAAATATTCCGGGAATTTTTTTGGGTTCTCACTTTAAACAGGCGCAGAGATTTGCCGATATGATAAAGGTGCGTATTTCGGCGGGGCAGGCGCAAAAATTTGGCGGTCTGGGGTATGATACAGAAAACCGCAAACTGGTTGGCGCGAATATTGCCATTGATATGGTGGAAATTTTAAACCggggaagagggac                                                                        >  1:299705/1‑238 (MQ=255)
                                                                                                                                                                          gTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGAGATGGTAAAGATTTTAAGCCTTCAAGGATTGAAAGa                                                                     <  2:40971/251‑1 (MQ=255)
                                                                                                                                                                                 gCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCAc                                                              <  2:114733/251‑1 (MQ=255)
                                                                                                                                                                                    aCCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACt                                                             <  2:278037/249‑1 (MQ=255)
                                                                                                                                                                                          gTTTTCGTGACCGCTGTGTATCGGCGGGCATTTGTGTGGAAATTATTTCGGGGATTTTTTCGGGTTCTAAATTTTAAAAAGCGAAGAAATTTTGCGATATGACCCACATGCGTATTTCGGCGGGGGTGGGGCAAATTTTCGACGGGCTGGGTGATGATTCCGAAAACCGCAAACCGGTTGGGGGGGATATTGCCATGGGTATGGGGAAGATTTTTAATCGTGAAGGAGGGAAAGAATTACACTTTTATACg                                                     >  2:322258/1‑251 (MQ=255)
                                                                                                                                                                                                               cggcggGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTGTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAg                                 <  2:536512/250‑1 (MQ=255)
                                                                                                                                                                                                                  cggGCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAc                             >  1:172666/1‑251 (MQ=255)
                                                                                                                                                                                                                    ggCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAcgc                           <  2:118001/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAcgc                           <  2:453606/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggCATTGATGTGGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTAcgc                           >  1:41757/1‑251 (MQ=255)
                                                                                                                                                                                                                              gtgGAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGcc                   <  1:445828/250‑1 (MQ=255)
                                                                                                                                                                                                                               tgGAAATTATTCCGGGAATTTTGGAGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGGGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAAACGTGCTGAAATGAGTTACGCGATTTGCCAt                 >  2:307642/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gAAATTATTCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATAc               <  2:349978/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttattCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTc                                                                                                                                                                    <  2:119954/98‑1 (MQ=255)
                                                                                                                                                                                                                                      ttattCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTc                                                                                                                                                                    >  1:119954/1‑98 (MQ=255)
                                                                                                                                                                                                                                          tCCGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTggggg       >  2:370911/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccGGGAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGt      <  1:20368/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gAATTTTGTCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTGGCCATACGCTGGGGGTTCGa  <  1:203013/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
ATGCTTCTGACCTGGTGACGCTGTTAAAAGAAGTGGCAGATTTCGATATCTCCGTGGCGGCGTATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGCCGGTATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAATATTGCCATGGATATGGTGAAGATTTTAAGCCGTGAAGGAGTGAAAGATTTCCACTTCTATACGCTTAACCGTGCTGAAATGAGTTACGCGATTTGCCATACGCTGGGGGTTCGA  >  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: