Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,001,389 G→A 5.4% intergenic (‑333/‑37) yigE ← / → corA DUF2233 family protein/magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,3890GA 5.4% 231.6 / 3.5 77intergenic (‑333/‑37)yigE/corADUF2233 family protein/magnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base A (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

CGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTG  >  NC_000913/4001153‑4001624
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           
cGTTTTGACGGAAAATTTTTTCGGGCTTAACTTGCGGGACAGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGaa                                                                                                                                                                                                                                           <  2:572288/239‑1 (MQ=255)
 ctgtTGACGGCAAAATTTTGCTGGCGTAACATGCCCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGaa                                                                                                                                                                                                                                           <  1:479896/237‑1 (MQ=255)
  tgtTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTACGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGCAAGCTACCTCAATTTTTCCCTTCCCGCTGCGTCAGGCTGTTTAAGTGTCTGAAACCCAATTTTGTTCGTTGTGCCAACCGa                                                                                                                                                                                                                                            >  2:308349/1‑236 (MQ=255)
  tgtTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGGGTCAGGCTGTTTAATGGGCTGAAACCCAATTTGTTTCTTTGTGCCACCCGa                                                                                                                                                                                                                                            >  1:299272/1‑236 (MQ=255)
        ccgCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTGGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCg                                                                                                                                                                                                                                    <  1:184970/236‑1 (MQ=255)
        cGGCAAAATTTTGCTGGCGTAACCTGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCg                                                                                                                                                                                                                                    <  1:772256/238‑1 (MQ=255)
        cGGCAAAAGTTTGGTTGCGTAACATGCGCGACCGATCCCTCTAAGAGGACATTCGCCTGGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCg                                                                                                                                                                                                                                    <  1:829545/238‑1 (MQ=255)
              aaTTTTGCTGGCGTAACCTGCGCGCACGATCCCTTTAAGAGGCCATTCGCCTTGGACACCCCCAGTAGATACTGGCTCTCTATCCTGTCATCCAGGATCAACTCCGAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATAtttt                                                                                                                                                                                                                             <  2:151364/239‑1 (MQ=255)
              aaTTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCTCTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATAtttt                                                                                                                                                                                                                             <  1:132958/239‑1 (MQ=255)
               aTTTTGCTGGCGTAACAGGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTa                                                                                                                                                                                                                            <  1:308349/239‑1 (MQ=255)
                       ggccgaaCAGGCCCGCACGATCCCTCTAAGAGGACATTCGCCGTGGACACACCCAGTGGATCCTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTAGCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGGGACCTCAATTTTTCCCGTCGCGCTGAGTCAGTCTGTTTAATGTTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTgg                                                                                                                                                                                                                    <  1:573811/234‑1 (MQ=255)
                                                  taggggACTTTCGGCTTTGACCCCCCCAGTACATTCTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTc                                                                                                                                                                                         <  2:812507/235‑1 (MQ=255)
                                                  aaGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTc                                                                                                                                                                                         <  1:371279/239‑1 (MQ=255)
                                                   agagGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTCTCCTGTCATCCAGGAACAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGGGTCTGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCCCCGAACTGTCCGATATTTTATGCCTTGGGAGTCCCGGTCATGCTGAGCGCATTTCa                                                                                                                                                                                        >  1:273455/1‑239 (MQ=255)
                                                   agagGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCCCATTTCa                                                                                                                                                                                        >  2:394858/1‑239 (MQ=255)
                                                       gACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGt                                                                                                                                                                                                                >  2:531278/1‑211 (MQ=255)
                                                       gACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGt                                                                                                                                                                                                                <  1:531278/211‑1 (MQ=255)
                                                         cATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTTAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCCCCGAACTGTCCGATATTTTAACAATTGGGAGTCCCGGTCATGCTGAGCGCATTTAAAATGGa                                                                                                                                                                                  >  2:130575/1‑239 (MQ=255)
                                                         cATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGa                                                                                                                                                                                  >  1:792148/1‑239 (MQ=255)
                                                          aTTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGa                                                                                                                                                                                  <  1:443538/238‑1 (MQ=255)
                                                                   ggACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCg                                                                                                                                                                         >  1:177812/1‑238 (MQ=255)
                                                                         acCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGAcc                                                                                                                                                                  <  1:553323/239‑1 (MQ=255)
                                                                              gTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCg                                                                                                                                                                >  1:20858/1‑236 (MQ=255)
                                                                                   tACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGTTATTTTAAGCCTTGGAGTTCCCGGCACTGCTTAGGCAATTTAAACTGGAAAATAACCGAAGGCCCCGGCTGAAAg                                                                                                                                                        >  1:249533/1‑239 (MQ=255)
                                                                                   tACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGGGCGCATTTCAACTGGGAAAAAACCGACTGACCCGGCGGGAAg                                                                                                                                                        >  2:533071/1‑239 (MQ=255)
                                                                                    aCTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTTGGTGCCCACCGAACTGTCCGGTATTTTATGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACGGACGGACCCGGCTGGAAgg                                                                                                                                                       >  2:791625/1‑238 (MQ=255)
                                                                                    aCTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCCTTGTGAGTCCAGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAg                                                                                                                                                        >  1:887527/1‑238 (MQ=255)
                                                                                    aCTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGATTCCCGGTCATGCTGAGCGCAttt                                                                                                                                                                                          >  2:71612/1‑204 (MQ=255)
                                                                                    aCTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCAttt                                                                                                                                                                                          <  1:71612/204‑1 (MQ=255)
                                                                                     cTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTATCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCCTTGGGAGTCCCGGTCATGCTGACCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGt                                                                                                                                                       >  2:354592/1‑238 (MQ=255)
                                                                                       ggCTTACTTTCCTTTCATCCAGGTTCAACCCCTAAGGTTTTCCCTTTGTTGTGTTAGCCTTAGCGGTTGTCAGCGACCCTAATTTTTCCCGTTGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGGTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGa                                                                                                                                                    <  2:299272/239‑1 (MQ=255)
                                                                                         cTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCCTTGGGAGTCCCGGTCATGCTGAGCGCATTTTAACTGGAAAATAAACGACTGCCCCGGCTGGAAGTCGaaa                                                                                                                                                  >  2:638853/1‑238 (MQ=255)
                                                                                         cTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGGGCGCATTTCAACTGGAAAAAAACAGACTTACCCGGCCGGAAGACGAag                                                                                                                                                  >  2:761703/1‑239 (MQ=255)
                                                                                            aCTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGGAACCCAATTTGCTTCTGTGTGCCCCCCGAACTGTCCGATATTTTACGCATTGGGAGTCTCGGTCATCCTGAGCGCATTTCAACTGGAAAATAAATGCCTGACCCGGCTGGAAGGCGAagag                                                                                                                                                >  1:432846/1‑238 (MQ=255)
                                                                                               aTCCTGTCATCCAGGATCAACTCCTAAGGCTTTCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTcac                                                                                                                                            <  1:130575/239‑1 (MQ=255)
                                                                                                  cTGTCATCCCGGATCAACTCCTAAGGCTAGCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAAc                                                                                                                                         <  1:899588/239‑1 (MQ=255)
                                                                                                    gacATCCCGGATCGCCTCCTAAGGTTTTCCCTGTTTGCTGATAGCTTGAGCGGTTGTCAGGGGCCGCAATTTTTCCCGTCGCGCTGAGTCAGGCTATTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCGTCGAACTGTCCGATATTGTAAGCATGGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAAccc                                                                                                                                       <  2:273455/237‑1 (MQ=255)
                                                                                                      cttCCAGGATCACCTCCTAAGGCTTTCCTTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGGGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCt                                                                                                                                     <  2:579082/237‑1 (MQ=255)
                                                                                                      cATCCCGGATCAACTCCTAAGGCTATCCCTTTTTGCTGAAAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCt                                                                                                                                     <  1:297407/239‑1 (MQ=255)
                                                                                                        tCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTg                                                                                                                                   >  1:543775/1‑239 (MQ=255)
                                                                                                          cAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTa                                                                                                                                 >  1:466432/1‑239 (MQ=255)
                                                                                                            ggATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGGCTGACCCGGCTGGAAGTCGAGGAGTCAAAACCCCTTGTaaa                                                                                                                               >  1:880327/1‑239 (MQ=255)
                                                                                                             gATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAgcgc                                                                                                                                                                                              >  1:463994/1‑175 (MQ=255)
                                                                                                             gATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTAAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAgcgc                                                                                                                                                                                              <  2:463994/175‑1 (MQ=255)
                                                                                                              aTCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGGCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAAt                                                                                                                              <  1:486606/238‑1 (MQ=255)
                                                                                                                  accccTCAGGTTTTCCCTTTTTGCTGTGAGCCTGGCGGGTTGTCAGGGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGGGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGt                                                                                                                         <  2:466432/236‑1 (MQ=255)
                                                                                                                  aCTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGt                                                                                                                         <  1:640557/239‑1 (MQ=255)
                                                                                                                       tAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGAATGGa                                                                                                                    >  1:571247/1‑239 (MQ=255)
                                                                                                                        aaGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAAAGGCACAACCCCTTGGAAATCCAGAAAGGAt                                                                                                                   >  2:299740/1‑239 (MQ=255)
                                                                                                                         aGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAATCGAAGAGTCACCACCCCTTGTAAATGCAGTATGGttt                                                                                                                  >  1:718062/1‑239 (MQ=255)
                                                                                                                              aTCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGTGCGCATTTCAACTGGTAAATAACCGACTGCCCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCt                                                                                                             >  1:628341/1‑239 (MQ=255)
                                                                                                                              aTCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAAACGACTGACCCGGCTGGAAGGCGAAGAGTCACAACCCCTTTTAAATGCAGTAAGGAATGAtat                                                                                                             >  2:851941/1‑237 (MQ=255)
                                                                                                                              aTCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCAAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCCCAACCCCTTGTAAATGCAGAATGGATTGATCt                                                                                                             >  1:133549/1‑239 (MQ=255)
                                                                                                                               tCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTgg                                                                                                                                                                                                                    <  2:700027/135‑1 (MQ=255)
                                                                                                                               tCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTgg                                                                                                                                                                                                                    >  1:700027/1‑135 (MQ=255)
                                                                                                                               tCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGCCGAAGAGTCACAACCCCTTGTAAATGCAGGATGGATTgttatt                                                                                                            >  2:236232/1‑234 (MQ=255)
                                                                                                                                cccTTTTTGCTGATAGCCTTTGGGGTTGTCAGCGACCTCAATTTTTACCGTCCGGTTGAGTCAGGCTGTTTAATGGTCTGAACCCCAATTTGTTTCTGTGTGCCAACCGAACTGTCCGATATTTTAAGCATTGGGATTCCCGGTCATGCTGAGCACATTTCAACTGGAAAGTAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTg                                                                                                           <  1:135124/239‑1 (MQ=255)
                                                                                                                                     tttgcctcttcccTTACTGGTTGGCAGTGGCCCCAATTGTTTCCTGCGCGCGTGTGTAGTCTGTTTAATGGACTGAACCCCAATTTTTTGCTGGTGGCCACCCCGACTGTCCGATATTTTATGCATTGGGAGTCCCGGTCATGCTGACCGCATTTCAACGGGAAAATAAGCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGaa                                                                                                      <  2:466064/228‑1 (MQ=255)
                                                                                                                                     tttGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTTCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGGAAATAACCGACTGACCAGGCTGGAAGTCGGAAAGCCACAACCCCTTGAAACAGCAGTAGGGATTGACCttttgaa                                                                                                       >  2:575794/1‑233 (MQ=255)
                                                                                                                                     tttGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCCCAACCCCTTGTAAATGCAGTATGGAt                                                                                                                   >  1:346410/1‑225 (MQ=255)
                                                                                                                                           gATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGCCTGACCCGGCTGGAAGTCGAAGAGCCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGgac                                                                                                >  1:160879/1‑239 (MQ=255)
                                                                                                                                                    aGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAACTAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATTTTTTCGAACCGGAGGCCGAagag                                                                                       >  2:768637/1‑239 (MQ=255)
                                                                                                                                                                 gACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGGTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTa                                                                           <  1:236232/238‑1 (MQ=255)
                                                                                                                                                                 caCCTTAATTTTTCCCGTCGCGCTGAGGCAGGCTTTTTAATGGTCCGAAACCCAATTTGTTTCTGTGGGCCCACCGAACTGGCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTGCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTa                                                                           <  2:269225/237‑1 (MQ=255)
                                                                                                                                                                   ccTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGTAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGCCTGCGGGTACaa                                                                        >  1:566407/1‑239 (MQ=255)
                                                                                                                                                                    ccaaaattttCCCGTTGGGCTGAGTCAGGCTTTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAAt                                                                       <  2:792148/233‑1 (MQ=255)
                                                                                                                                                                          ttttCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGGCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAAc                                                                 <  1:127481/239‑1 (MQ=255)
                                                                                                                                                                                 gTCGCGCTGAGTCAGGCTGTTTAATGGACTGAAACCCAATTTGTTTCTGTGTGCCCACCAAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAGTAACCGACTGACACGGCTGGAGGTCGACGAGTCACAACCCCTTGTCAATGCAGTATGGATTGATCTTTTCCAACCGGAAGACGACGAACGACTGCGCGCAAAACCTAAAATTggacc                                                          >  2:151187/1‑236 (MQ=255)
                                                                                                                                                                                                     tttAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCAAAAGCCTGGCAAaccgcccg                                       >  2:399562/1‑239 (MQ=255)
                                                                                                                                                                                                        aattgttttaaaCCCAATTTTTTTCTGTGTGCCCCCCAAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGaa                                    <  2:133549/230‑1 (MQ=255)
                                                                                                                                                                                                               tGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAAAGCCTGGCAACCCGCCCGGAACTGGAAg                             >  1:655341/1‑239 (MQ=255)
                                                                                                                                                                                                                     cacaTTTGTTTCTGTGTGCCAACCGAACTGTCCGAAAGTATAAGCATTGGGAGTCCCGGTCATGCTGAGCACATTTCCACTGGAAAATAACCGACTGCCCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTTCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCg                       <  1:411982/236‑1 (MQ=255)
                                                                                                                                                                                                                       aaTTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGaa                     >  1:721493/1‑239 (MQ=255)
                                                                                                                                                                                                                               tCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACACTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAGGACATCGAAGCATCgg              >  1:718852/1‑238 (MQ=255)
                                                                                                                                                                                                                               tCTGGGTCCCCCCCGAGATGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGGCGAAGAGTCACACCCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGc             <  2:20858/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tgtgCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGCTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCAGCCCGGAACTGGAAGACATCGAAGCATCggaccgt         >  1:828578/1‑234 (MQ=255)
                                                                                                                                                                                                                                          aCCAAACTGTCCGATATTTTAAGCAGTGGGATGCCCGGTCATGCTGAGCGCATTTCAACCGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTg  <  1:151187/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           
CGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTG  >  NC_000913/4001153‑4001624

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: