Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,350 T→C 13.1% intergenic (+37/+3) metE → / ← ysgA 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,3500TC13.1% 179.4 / 3.2 84intergenic (+37/+3)metE/ysgA5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base T (28/43);  new base C (3/9);  total (31/53)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGC  >  NC_000913/4015122‑4015576
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
tATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGGGCAGGCGGCGCAGAACTTGCGTCGGGGGGAAAATCCAAACCGGGTGGTAATACCCCCCGGGCTTTtct                                                                                                                                                                                                                                     >  1:626761/1‑229 (MQ=255)
 aTGACTTTCACTCGCTAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCGGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACGTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                                                                                                                                                                     <  1:138068/228‑1 (MQ=255)
    aCATTCACTCGCCAAACGTACCGAGCGAGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCAtt                                                                                                                                                                                                                                 <  1:625194/229‑1 (MQ=255)
        tCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCCCGGCTGGCCAGAAACCCGCGCGCGACTGGCCAACCTGGGGCAGGCGGCGCAGAACTTGCGTCGGGGGGAAAAACCCAACCGGGTGGTAATAACCCCCGGGCTTTTCTCAttcaa                                                                                                                                                                                                                              >  2:698345/1‑225 (MQ=255)
           cgccccAAGCGTACCGAGCGGGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAACCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGa                                                                                                                                                                                                                          <  1:25570/225‑1 (MQ=255)
            tCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGGGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGGCTTTTCTCCTTACAGCGAc                                                                                                                                                                                                                         >  1:49380/1‑229 (MQ=255)
              gCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGAACTGGCGAACAGGGTGCAGGCGGCGCCGAACTTGCGCCGGGGGGAAAACCCAACACCGGGGGGAAAACCCCCCGGCCTTTTTCCCTTACaaggcatt                                                                                                                                                                                                                       >  2:369869/1‑222 (MQ=255)
                     gttccgggggTGGAATGGATTGAAGCCTTGATGAAAAAAGCGGTAAAACGCATTCCGGCAGGGCGCCTGGTGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCcca                                                                                                                                                                                                                <  2:414611/220‑1 (MQ=255)
                          gggggttgAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATa                                                                                                                                                                                                           <  1:151885/222‑1 (MQ=255)
                          gAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAACACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTCCCGCGGGGGTAAAAGCCCAACCGGGGGGGAAAACCCCCAGGTCTTTTCTCCTTACAGCGACCCGTTCccaccca                                                                                                                                                                                                            >  2:482112/1‑226 (MQ=255)
                          gAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGACCTTGCGTCGGGGGGAAAATCCAAACCGTGTGGTAATACCACCCGGGCTTTTTTCCTTACACCGACTTCTTCCCACCAt                                                                                                                                                                                                            >  1:149973/1‑228 (MQ=255)
                          gAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGGAAAATCCAAACCGGGGGGTGAGCCCCCCCGGCCTTTTCTCCtttccggcgcctcttccccccaa                                                                                                                                                                                                            >  2:413408/1‑207 (MQ=255)
                          gAGCGTGGAAGGGATTGAAGCCTGGGTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATa                                                                                                                                                                                                           <  1:281593/229‑1 (MQ=255)
                               tGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCCCAGAACTTGCGCCGGGGGTAAAATCCAAACCGGGTGGGAATAACCCCCGGCCTTTTTTCATTCCAGCGCCCTCTTCCCCCCCTACCGc                                                                                                                                                                                                       >  2:623902/1‑228 (MQ=255)
                               tGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGGGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCt                                                                                                                                                                                                      <  1:554075/229‑1 (MQ=255)
                                 gAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCCGGGGCAAAATCCAAACCCGCGGGGAAACCCACCCGGGCTTTTTCCCTTACACCGCCTTCCTCCCACCATACCGCCTa                                                                                                                                                                                                    >  2:243145/1‑229 (MQ=255)
                                   tgtgATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGGGCGCCTGGGGGTCAACCCGGACTGGGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                                                                                                                  <  2:626760/226‑1 (MQ=255)
                                   agggATTGAAGCCTGGCTGAAGAAAGCGGCAAAACGCATTCCGGGAGAGCGCCTGTGGGTCAACCCGGACTGGGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                                                                                                                  <  2:285010/227‑1 (MQ=255)
                                       attggaagctggcggaagaaAGCGGCAAAACGCAGTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGGCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCTGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCa                                                                                                                                                                                               <  2:207527/214‑1 (MQ=255)
                                       aTTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCa                                                                                                                                                                                               >  1:26910/1‑228 (MQ=255)
                                         tgaacgcTTGCTGAAGAAAGCGGCAAAACAGATTCCGGGAGAGCGCCTTTGGGTCAACCCGGACGGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGTGCGGCACTGGCGAACATGGGGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAAGACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATTCTGCTTAAACCAtt                                                                                                                                                                                             <  2:472591/222‑1 (MQ=255)
                                           aaGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGGGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGGGGGAATACCACCCGGGCTTTTCTCATTAAAAAGACTTCTTCCCACCAAAATCCTTAAACCCTTcc                                                                                                                                                                                           >  2:685354/1‑228 (MQ=255)
                                            aGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGGCTTTTCTACTTACAGCAACTTCTTCCCACCATACTCCTTACACCATTCCac                                                                                                                                                                                         >  1:339466/1‑228 (MQ=255)
                                                ttGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCCGCGGGGGTAAACTCCAAACCGGGCGGTAATCCCACCCGGTCTTTTCCCCTTACACCGAATTCCTCCCACCCCACTGCATAAAACATTTCCGCAt                                                                                                                                                                                     >  2:696585/1‑228 (MQ=255)
                                                ttGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAAAACTTGCGTCGGGGGTAAAATCCAAACCGGGGGGTAATCCCACCCGGTCTTTTCTCATTTCAGCGACTTCTTCCCCCCATACCCCTTAAACCATTCCAGCaa                                                                                                                                                                                     >  2:738266/1‑227 (MQ=255)
                                                         aaaGCGGCAAAACGCATTCCGGCAGAGCGCCTGGGGGTCAACCCGGACTGGGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                                                                                                                                                                           <  1:55090/229‑1 (MQ=255)
                                                           cgccgCAAACCCCATGCCGGAAGAGCGCCTGTGGGTCAACCCGGGCGGTGGCCGGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACTCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                                                                                                                                                                         <  2:322009/225‑1 (MQ=255)
                                                           aGCGGCAAAACGCCATCCGGCAGGGCGCCGGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCTAACATGGTGCAGGCGGCCCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                                                                                                                                                                         <  2:507153/229‑1 (MQ=255)
                                                                cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCCGAACTTGCGTCGGGGGTAAGATCCAAACCCGGGGGTAATACCACCCGGGCTTTTTTCTTTAAAACGGATTCTTACCCCCATAACTCTTAAACCCTTCCAGTCACACCTCCCAGCCAactt                                                                                                                                                                    >  2:693340/1‑229 (MQ=255)
                                                                cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGAAAAATCCAAACCGGGTGGTAATACCACCCGGCCTTTTCTCATTACAGCGACTTCTTCCCACCCTACTGCTTAAACCATTCCAGCATACCCTGCCAGCCAtctt                                                                                                                                                                    >  1:356861/1‑229 (MQ=255)
                                                                    aCGCATGCCGGGAGGGCGCCTTTGGGTCCACCCGGACTGTGGCCTGAAACCGCGGGGCGGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg                                                                                                                                                                 <  2:219660/228‑1 (MQ=255)
                                                                       cagtCCGGGAGGGCGCCTGTGGGTGAAGCCGGCCTGGGGCCCGAAAACGGGGGGCTGGCCAGAAACCCGCGCGGCACTGGCGAAAATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGa                                                                                                                                                             <  2:383190/226‑1 (MQ=255)
                                                                          tgcggcggGGCGCCTGGGGGTCAACCCGGACTGTGGGCCGAAAACGCGCGGCTGGCCAGAAACCCGGGGGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATg                                                                                                                                                           <  2:247256/226‑1 (MQ=255)
                                                                           ccGGCGGAGCGCCTGTGGGTCACCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGATAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg                                                                                                                                                         <  2:228273/229‑1 (MQ=255)
                                                                            cgggagggggcCTGTGGGGCCCCCCGGACTGGGGCCTGAAAACCGGCGGGCGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCgg                                                                                                                                                        <  2:615350/220‑1 (MQ=255)
                                                                            cGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg                                                                                                                                                         <  1:369869/228‑1 (MQ=255)
                                                                                    cgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTTGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGGGCCTTCTTCCCACCATAATGCTTTAACCATTCCCGCATACGCTGCCAGCCCTCTCCTGACGATGCGGCATGATaa                                                                                                                                                >  2:744789/1‑228 (MQ=255)
                                                                                    cgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGGAAAATCCAAACCGGGTGGGAATACCACCCGGTCTTTTCTCATTCAAGCGAATTCTTCCCACCATCATGCTTAAACCCTTCCACCATAACCTCCCCGCCAGCTTCTGCAAAGCCGGCAtaaaag                                                                                                                                                >  2:395329/1‑224 (MQ=255)
                                                                                       cTGGGGGTCAACCCGGACGGTGGCCCGAAAACGCGCGGCTGGCCAAAAACCCGGGCGGCACTGGCGAACAGGGTGCAGGCGGCGCAGGACTTGCGTCGGGGGTAAAAGCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTTCCAGCCATCTTCTGCAGATGCGGCATGATAGCTc                                                                                                                                             <  2:577053/229‑1 (MQ=255)
                                                                                                 ccccGGACTGGGGCCTGAAACCGCGGGGGTGGCCAGAAACCCGCGCGGCACTGGCGAAAATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAAt                                                                                                                                   <  2:339466/228‑1 (MQ=255)
                                                                                                 aCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGCAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTTTGCAGATGCGGCATGATAACTCGGGCGCTAAt                                                                                                                                   >  1:503130/1‑229 (MQ=255)
                                                                                                       aCTGTGGCCTGAAAACGCGCGGCTGTCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCCGGGGTAAAACCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTAAAGCGACCTCCTCCCACCTTACTGCCTAAACCATTCCCGCCAACGCTGCCTGCCACCCTCTGCAGATCCGCCACGAAAGCCCGTCCGCTAAtaccac                                                                                                                              >  2:132324/1‑223 (MQ=255)
                                                                                                        cgtgggCCTGATACCGGGCGGTGGGCCAGAAACCCGGTCGGCACTGGCGAACAGGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCACATGACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGt                                                                                                                             <  2:431134/224‑1 (MQ=255)
                                                                                                               cggAAAACGCGGGGCCGGCCAGAAACCCGGGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgca                                                                                                                     <  1:599063/227‑1 (MQ=255)
                                                                                                                       gcgcgGGTGGCTCGAAACCCGGGCGGGCCTGGCGGACAAGGTGCGGGGGGCGCAGAACTTGGCGCGGGGGTAAAAGCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGCCTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCgg                                                                                                              <  2:567609/228‑1 (MQ=255)
                                                                                                                        cgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                                                                                                                                                                           <  2:195211/166‑1 (MQ=255)
                                                                                                                        cgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                                                                                                                                                                           >  1:195211/1‑166 (MQ=255)
                                                                                                                        cgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGCCGATAATCAGCGTTGAATGCATGCCCGGc                                                                                                             >  1:647517/1‑228 (MQ=255)
                                                                                                                         gcgGCTGGCCAAAAACCCGCGGGGCACTGGCGAACATGGGGCGGGCGCCGCAGACCTTGCGTTGGGGGTAAAATCCAATCCGGGTGGGAATACCACCCGGGCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCg                                                                                                            <  2:618678/228‑1 (MQ=255)
                                                                                                                                 ccAAAACCCCGGGCGGCCCTGGCGAACATGGGGCAGGCGGAGCAGAACTTGCGTCGGGGGTAACATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCGTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGt                                                                                                     <  2:31839/227‑1 (MQ=255)
                                                                                                                                     aaaCCCGGGCGGCACTGGCGAACATGGTGAGGGCGGCGCAGAACTGGCCTCGGGGGTAAAATCCAAACCTGGTTGTAATACCTCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCAGTCGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGa                                                                                               <  1:744790/229‑1 (MQ=255)
                                                                                                                                          cgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATTCGGCATGATAGCTCGGGCGATAATCACCGTTGAATGCATGCCCGGCGTCCCGGTACACggaaat                                                                                           >  1:624462/1‑223 (MQ=255)
                                                                                                                                                  tcTGGCCAAAATGGTGCATGCGGCTCAAAACTGGCGTCGGGGGTAAAATCCAACCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGCCTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCttt                                                                                  <  2:49380/228‑1 (MQ=255)
                                                                                                                                                           cATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGCTAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTag                                                                          >  1:689774/1‑228 (MQ=255)
                                                                                                                                                                   aggctgtcccgcaactGCGTCGGGTGTAAAACCCAAACCGGGTGGTAATAACACCCGGTCTTTTCTGATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCa                                                                 <  2:485485/214‑1 (MQ=255)
                                                                                                                                                                   aGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGc                                                                  >  1:602678/1‑228 (MQ=255)
                                                                                                                                                                   aGGCGGCCGAGAACTTGCGTCGGGGGTAAAATCCAACCCGGGTGGTAATCCCCCCCGGTCTTTTCCCATTACAGCCACTTCCTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCa                                                                 <  2:26910/229‑1 (MQ=255)
                                                                                                                                                                    ggctggcgaGAACTTGCGTCGGGGATAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCa                                                                 <  1:685355/220‑1 (MQ=255)
                                                                                                                                                                         cgcAGAACTTGCCTCGGGGGTAAAATCCAAACCGGGTGGTAATTCCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGcc                                                            <  1:30966/228‑1 (MQ=255)
                                                                                                                                                                           cAGAACTTGCGTCTGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTgg                                                         <  1:395329/229‑1 (MQ=255)
                                                                                                                                                                                tttcgctgtggggGAAGACCCAAGCCGGGGGGGAATACCCCCCGGTCTTTTCCCATTTCGGGGCCTTTTTCCCGCCATCCTGCGTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCa                                                     <  2:356861/219‑1 (MQ=255)
                                                                                                                                                                                cTTGCTGCGGGGGTAAAATCCAAACCGGGGGGTAATACCCCCCGGTCTTTTCCCAGTACGTCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCGTTCGCATTAGCAGCCCGCAGCGCCTGGCGCa                                                     <  2:640006/228‑1 (MQ=255)
                                                                                                                                                                                          ggttAAACTCCAAACCGGGTAGGAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACg                                          <  1:698346/226‑1 (MQ=255)
                                                                                                                                                                                          ggttAAAATCCAAACCGGGTGTTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTCCCCACCATACTGCTTAAACCATTCCAGCATCCGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACg                                          <  1:696586/226‑1 (MQ=255)
                                                                                                                                                                                          gggTAAAATCCAAACCGGGTGGTAATTCCACCCGGTCTTTTCCCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACg                                          <  1:522326/229‑1 (MQ=255)
                                                                                                                                                                                               aaaTCCAAACCGGGTGGTAATACCACCAGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCtctc                                     <  1:302699/229‑1 (MQ=255)
                                                                                                                                                                                                    cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCg                                >  1:21833/1‑229 (MQ=255)
                                                                                                                                                                                                     caaCCGGGTGGTAATCCCACCCGGTCTTTTCCCATTCCAGCGACTTCTTCCCCCCATACCGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCg                                <  2:624461/227‑1 (MQ=255)
                                                                                                                                                                                                     aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCg                                <  1:568690/228‑1 (MQ=255)
                                                                                                                                                                                                     aaaCCGGGTGGTAATACCACCCGGTCTTTTCCCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCg                                <  1:630187/228‑1 (MQ=255)
                                                                                                                                                                                                      aaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCTGGCGATAATCAGCGCTGAATGCATGCCCGGCGTCCGGGGCCACGATAATCTCTGCTTTACCATTAGCAGCCCCCAGCGCCTGGCGCATGGTTTGaaactcctcctacgg                               >  2:490818/1‑215 (MQ=255)
                                                                                                                                                                                                        cccGTGTGGTAATTCCCCCCCGTCTTTTCTCATTACAGCGTCTTCTTCCCACAATACTGCTGAAACCATTCCAGCATCCGCTGCCAGCCATTTTCTGCAGTGGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGaa                             <  2:695782/228‑1 (MQ=255)
                                                                                                                                                                                                         ccGGGTGTTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCGTCCCCCCATACTGCTTAAACCATTCCAGCAGACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGaa                             <  2:273322/228‑1 (MQ=255)
                                                                                                                                                                                                         ccGGGTGGTAATTCCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCCCCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGaa                             <  1:175300/228‑1 (MQ=255)
                                                                                                                                                                                                             gtggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGACTGGCGCATGGTTTCAACGCCATCCTGCGGAATggt                         >  1:312028/1‑226 (MQ=255)
                                                                                                                                                                                                               ggtAATTCCCCCCGGTCTTTTCCTATTACCGGGCCTTCTTCCCACCAGACTGGTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGtt                      <  2:85006/229‑1 (MQ=255)
                                                                                                                                                                                                               ggtAATACCCACCGGTCTTTTCCCAATCCAGCGACTTCGTCCCACCATTCTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTg                        <  1:243145/227‑1 (MQ=255)
                                                                                                                                                                                                                 tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGtt                      >  1:85006/1‑227 (MQ=255)
                                                                                                                                                                                                                         cccgtgcTTTTCTCATTCCAGGGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATGCGCTGCCAGCCATCTTCTGCAGATGGGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGacc             <  2:463382/222‑1 (MQ=255)
                                                                                                                                                                                                                         cccGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCCGATGCGGCATGATAGCTCGGGCGCTAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCCCTGCTTTCGCATTAGCACCCCGCAGCCCCTGGCCCCTGGTTTCAACCCTCCCCTTCGGACTGCACTCATCCTGacc             >  2:73784/1‑228 (MQ=255)
                                                                                                                                                                                                                         cccGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGACTGGTTTCAACGCTATCCTGCGGAATGCTGTTATCCTGacc             >  1:237025/1‑228 (MQ=255)
                                                                                                                                                                                                                         cccGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCACCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGacc             >  1:297827/1‑228 (MQ=255)
                                                                                                                                                                                                                         cccGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTGAAACCATTCCAGGATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGacc             <  1:623903/228‑1 (MQ=255)
                                                                                                                                                                                                                         cccGGGCTTTTCCCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGacc             <  1:482112/228‑1 (MQ=255)
                                                                                                                                                                                                                              tCTTTTCTCATTACAGCGACTTCTCCCCAACAGACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAGCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCata       <  2:647516/229‑1 (MQ=255)
                                                                                                                                                                                                                               tttttCTCATTCCCGCGCCGTCTTCCCCCCAGACTGCTTAAACCCATCCAGCATACGCGGCCAGCCGTTTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCatat      <  2:13092/228‑1 (MQ=255)
                                                                                                                                                                                                                                   tctcATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGGGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGc  >  1:722936/1‑229 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   
TATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGC  >  NC_000913/4015122‑4015576

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