Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,094 G→C 5.5% intergenic (+61/+167) glxK → / ← allE glycerate kinase II/S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,0940GC 5.5% 975.3 / 35.4 280intergenic (+61/+167)glxK/allEglycerate kinase II/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (137/127);  new base C (7/9);  total (144/136)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

Only 100 of 282 total aligned reads displayed.

CGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCT  >  NC_000913/542822‑543340
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
tggaggaGTTTTATGTAGCGGTGTGTGGGATGGCTGGCGTATTTGGTGATGGTGTGTGTGTGGGGGACCAGTATGGCATTGACGCGGTATTCAGCATTTTGCCTCGGCTGGCACCTTTGGCCGAAGTGCTCGCCAGCGGTGAAACCAAGCTCTTCAACAGTGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGc                                                                                                                                                                                                                                                       <  4:253440/268‑1 (MQ=255)
      aGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATg                                                                                                                                                                                                                                                 <  2:486244/274‑1 (MQ=255)
                                           gggTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCCCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAACCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTTCAGTACCCCAGCCGGCAGCGGCATGAGTAACCTATCCATCCTGAAAATTCCGCAACTTCtaga                                                                                                                                                                                                           >  4:195368/1‑273 (MQ=255)
                                                                           gCATTGACGCGGTATTCAGCATTTTGCCTCGGCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcg                                                                                                                                                                            <  2:338795/274‑1 (MQ=255)
                                                                                      gTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCa                                                                                                                                                                  <  2:195371/273‑1 (MQ=255)
                                                                                                   tGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAACTTCTCTATCTTGCAGAGATCATGTATGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaacg                                                                                                                                               >  1:99517/1‑277 (MQ=255)
                                                                                                         gTCTTGCACCTTTAGCCGAAGTGCGCGCCAGGGGGGAAAACAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGcc                                                                                                                                             <  4:399775/275‑1 (MQ=255)
                                                                                                                    ttAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaa                                                                                                                                                 >  2:399779/1‑260 (MQ=255)
                                                                                                                                cTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGc                                                                                                                                                                                                                                                    >  5:1556560/1‑149 (MQ=255)
                                                                                                                                  cGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                                                                                                                                                                         >  2:435438/1‑174 (MQ=255)
                                                                                                                                    ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                                                                                                                                <  6:553546/149‑1 (MQ=255)
                                                                                                                                    ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                                                                                                                                <  5:924832/149‑1 (MQ=255)
                                                                                                                                    ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                                                                                                                                <  7:664365/149‑1 (MQ=255)
                                                                                                                                    ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                                                                                                                                <  5:752870/149‑1 (MQ=255)
                                                                                                                                         ggTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAcc                                                                                                                                                                                                                                           >  7:496176/1‑149 (MQ=255)
                                                                                                                                                 cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcctt                                                                                                                                                                                                                                   >  5:1297564/1‑149 (MQ=255)
                                                                                                                                                  aagctctTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTAATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                                                                                                                                                                  >  8:1130566/4‑149 (MQ=255)
                                                                                                                                                   aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                                                                                                                                                                 <  8:62057/149‑1 (MQ=255)
                                                                                                                                                    tctcTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCa                                                                                                                                                                                                                                          <  8:58629/139‑1 (MQ=255)
                                                                                                                                                       ctTCAAAAGGGCGCGAAAAATTGCCTGGGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAg                                                                                                <  4:127996/274‑1 (MQ=255)
                                                                                                                                                        ttCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGaaa                                                                                                                                                                                                                            >  7:341528/1‑149 (MQ=255)
                                                                                                                                                          cAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAACCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCCCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAACGAAGCGCATCATGCAATGTTACAAAAAAAGCAACGGATAAAAACGTTGCAAAATCCAACATTACTAAAATAAATAAGGcag                                                                                            >  4:304178/1‑273 (MQ=255)
                                                                                                                                                          cAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGa                                                                                            >  2:547456/1‑275 (MQ=255)
                                                                                                                                                                  cgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtatata                                                                                                                                                                                                         >  4:41373/1‑158 (MQ=255)
                                                                                                                                                                  cgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACtt                                                                                >  1:155969/1‑279 (MQ=255)
                                                                                                                                                                     gAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTc                                                                                                                                                                                                               >  8:1017329/1‑149 (MQ=255)
                                                                                                                                                                     gAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAg                                                                                                                                                                                                                                               <  7:1330484/117‑1 (MQ=255)
                                                                                                                                                                     gAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAg                                                                                                                                                                                                                                               >  8:1330484/1‑117 (MQ=255)
                                                                                                                                                                      aaaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGc                                                                                                                                                                                                     <  2:41373/158‑1 (MQ=255)
                                                                                                                                                                       aaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                                                                                                                             <  6:1168336/149‑1 (MQ=255)
                                                                                                                                                                             gCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGGTTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATAtt                                                                                                                                                                                                       <  6:628109/149‑1 (MQ=255)
                                                                                                                                                                                  cgcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCACAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCaga                                                                                                                                                                                                  <  6:1052673/149‑1 (MQ=255)
                                                                                                                                                                                   gcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCagag                                                                                                                                                                                                 >  6:268388/1‑149 (MQ=255)
                                                                                                                                                                                   gcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGAATTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCagag                                                                                                                                                                                                 >  6:1719664/1‑149 (MQ=255)
                                                                                                                                                                                      aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                                                                                                                              >  5:54144/1‑149 (MQ=255)
                                                                                                                                                                                      aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                                                                                                                              >  5:150670/1‑149 (MQ=255)
                                                                                                                                                                                      aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                                                                                                                                >  8:844895/1‑99 (MQ=255)
                                                                                                                                                                                        tAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                                                                                                                                                            <  6:1099779/149‑1 (MQ=255)
                                                                                                                                                                                        tAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                                                                                                                                                            <  5:1422576/149‑1 (MQ=255)
                                                                                                                                                                                            aTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAg                                                                                                                                                                                        <  7:306669/149‑1 (MQ=255)
                                                                                                                                                                                              aGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGc                                                                                                                                                                                      <  6:905079/149‑1 (MQ=255)
                                                                                                                                                                                               ggTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc                                                                                                                                                                                                                         >  4:784358/1‑113 (MQ=255)
                                                                                                                                                                                                gTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt                                                                                                                                                                                    >  8:126534/1‑149 (MQ=255)
                                                                                                                                                                                                 tCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCg                                                  >  1:187317/1‑278 (MQ=255)
                                                                                                                                                                                                  cAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGa                                                                                                                                                                                  <  5:300248/149‑1 (MQ=255)
                                                                                                                                                                                                  cAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGa                                                                                                                                                                                  >  7:268176/1‑149 (MQ=255)
                                                                                                                                                                                                    gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                                                                                                                                                                                               >  4:481677/1‑134 (MQ=255)
                                                                                                                                                                                                     ggAAGTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATGTAGATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTTAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTTATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAATATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGa                                                 <  4:547452/275‑1 (MQ=255)
                                                                                                                                                                                                          tAAAAACTAACCCTTACGTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTAAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACa                                                                                                                 <  1:119260/206‑1 (MQ=255)
                                                                                                                                                                                                            aaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaagc                                                                                                                                                                        >  6:1363530/1‑149 (MQ=255)
                                                                                                                                                                                                               aCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCa                                                                                                                                                                     <  7:1261661/149‑1 (MQ=255)
                                                                                                                                                                                                                   aCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCt                                                                                                                                                                                                              <  1:23022/104‑1 (MQ=255)
                                                                                                                                                                                                                   aCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCt                                                                                                                                                                                                              >  3:23022/1‑104 (MQ=255)
                                                                                                                                                                                                                     ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                                                                                                                              <  7:736143/118‑1 (MQ=255)
                                                                                                                                                                                                                     ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                                                                                               >  7:1031981/1‑149 (MQ=255)
                                                                                                                                                                                                                      cTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTa                                                                                                                                                                                                                                             >  8:1475347/1‑70 (MQ=255)
                                                                                                                                                                                                                      cTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTa                                                                                                                                                                                                                                             <  7:1475347/70‑1 (MQ=255)
                                                                                                                                                                                                                       ttACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa                                                                                                                                                             >  8:319057/1‑149 (MQ=255)
                                                                                                                                                                                                                             ttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTcacc                         <  2:523257/275‑1 (MQ=255)
                                                                                                                                                                                                                             ttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCAGGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTcacc                         <  2:304182/275‑1 (MQ=255)
                                                                                                                                                                                                                                 ctcttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaa                                                                                                                                                   >  8:696739/6‑149 (MQ=255)
                                                                                                                                                                                                                                 aGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaa                                                                                                                                                   >  8:1679616/1‑149 (MQ=255)
                                                                                                                                                                                                                                     gATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAg                                                                                                                                               >  7:932881/1‑149 (MQ=255)
                                                                                                                                                                                                                                     gATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACg                      <  4:600241/270‑1 (MQ=255)
                                                                                                                                                                                                                                     gATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACg                      >  2:600245/1‑270 (MQ=255)
                                                                                                                                                                                                                                         cGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                                                                                                                                  <  4:27864/46‑1 (MQ=255)
                                                                                                                                                                                                                                          gATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACg                                                                                                                                          >  6:319071/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aTTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACgg                                                                                                                                         >  8:891189/1‑149 (MQ=255)
                                                                                                                                                                                                                                                tatCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATaa                                                                                                                                    >  6:645007/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                                                                                                                                    <  3:180047/81‑1 (MQ=255)
                                                                                                                                                                                                                                                     cATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGt                                                                                                                               >  5:163603/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     cATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGt                                                                                                                               <  6:150670/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     cATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCt  >  4:449780/1‑274 (MQ=255)
                                                                                                                                                                                                                                                      aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTg                                                                                                                              <  8:496176/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                                                                                                                              <  5:1143934/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                                                                                                                              <  5:639166/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTaaaaa                                                                                                       >  2:25942/1‑172 (MQ=255)
                                                                                                                                                                                                                                                        aaGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGc                                                                                                                            >  7:415707/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         agagTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTGGct                                                                                                                           <  7:857780/149‑2 (MQ=255)
                                                                                                                                                                                                                                                          gagTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGNATCAGGCAATGTTACAAAAAAAANNACGGTATANACCGTGGCAAAATCCAACATAGCTAAAAATANACAGGCGAg                                                                                           >  3:65388/1‑180 (MQ=255)
                                                                                                                                                                                                                                                          gagTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTGGctt                                                                                                                          <  6:1659147/149‑3 (MQ=255)
                                                                                                                                                                                                                                                             tGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAgaga                                                                                                                                                                                                <  1:49242/76‑1 (MQ=255)
                                                                                                                                                                                                                                                             tGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa                                                                                                                                                                                             >  2:140994/1‑79 (MQ=255)
                                                                                                                                                                                                                                                               cAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATcc                                                                                                                     <  8:1589911/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTa                                                                                                                                                                                         <  4:140992/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa                                                                                                                   <  6:754297/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa                                                                                                                   >  8:379934/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTg                                                                                                                                                                                   >  7:1627816/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                     tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa                                                                                                               >  8:110398/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     tCATGCCGCATCCGGCATGAGTACTATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa                                                                                                               >  5:841077/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                                                                                                                              >  8:749004/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                       aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAAt                                                                                                                       <  8:1524548/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAAt                                                                                                                       >  7:1524548/1‑139 (MQ=255)
                                                                                                                                                                                                                                                                       aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGc                                                                                                             >  5:472495/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAgaga                                                                                                                                                                                                <  2:715101/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCt                                                                                                            >  5:1306512/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         gCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTa                                                                                                           <  6:130163/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              aTGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAaat                                                                                                      >  7:1174361/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                gCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGc                                                                                                                                                                                                     <  2:473024/52‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAaataa                                                                                                    <  7:78884/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
CGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCT  >  NC_000913/542822‑543340

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: