Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,086 C→A 9.3% Q183K (CAA→AAA)  waaA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,0860CA 9.3% 222.5 / 2.8 91Q183K (CAA→AAA) waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base C (51/30);  new base A (2/8);  total (53/38)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTA  >  NC_000913/3808859‑3809307
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                              
cGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTgcgc                                                                                                                                                                                                                                >  1:694343/1‑228 (MQ=255)
   gaaacaggCCCGATGCAATTATCCGGTTCCTGAATAACGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTAGGGCCTATCCTGATGGCGGCGCTACATAAACGTAAAATCCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCaaa                                                                                                                                                                                                                             <  2:409208/221‑1 (MQ=255)
   aTGATCTGCCCGGTGGACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCaaa                                                                                                                                                                                                                             <  1:683899/228‑1 (MQ=255)
         tGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTATGCCCGCTAGGCCGCAGGTTATGCCAACCTGGGTAAATTCGTCCGTCGCTTGCTGCCTCGTTTTACGCGGATTGCTGCGaaa                                                                                                                                                                                                                              >  2:683899/1‑221 (MQ=255)
           cccGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATAGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTACATTCGTCCGTCGCTTGCTGCGCCGGATAACGCTCATTGCTGCGAaaaacaaaaaa                                                                                                                                                                                                                     >  2:310881/1‑221 (MQ=255)
            ccGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCCAAGGTTAGGCCAAACTGGGTAAATTCGTCCTCCGCTTGCTGGGCCGTATTACGCTGTTTGCGGCGCAAAATgaagaaa                                                                                                                                                                                                                    >  2:459301/1‑227 (MQ=255)
             agacagaCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACCGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAagaaga                                                                                                                                                                                                                   <  1:240314/222‑1 (MQ=255)
                   acgcAACCGTTTCCTGAATAAAGGCGACCCTAAACGGGTGTTGATTATGGAAACCGAACGATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGc                                                                                                                                                                                                             <  1:392076/225‑1 (MQ=255)
                     tCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCCCTTGCTGCGTCGTATTACCCTGCTGGCTGCCCAAAATGTagaagacctgccc                                                                                                                                                                                                            >  1:409208/1‑220 (MQ=255)
                          cGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGATATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGGGCACGttt                                                                                                                                                                                                       >  1:396765/1‑227 (MQ=255)
                             ttCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGCATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCCCGTTTtgtg                                                                                                                                                                                                   >  1:482540/1‑228 (MQ=255)
                              tCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATAGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCCCTGCCCGCTCGGCCGCAGGGTATGCCAAACAGGGGAAATTCGTCCGTCGCTTGCTGCGTCGGTTTACGCTGCTTTGTGTGCAAAATGAAGCAAATGGGGCACGTTCtgac                                                                                                                                                                                                   >  2:240314/1‑225 (MQ=255)
                                        aGTCGCCCCTAAACTGGGTTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGgcgc                                                                                                                                                                                        <  1:466604/228‑1 (MQ=255)
                                        aGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGCCTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCCCTTGCTGCGTCGGATTACGCTGATTGCTGCGCAAAAAGATGAAGATGTGGCCCGTTTTGTGGCGCTgggggg                                                                                                                                                                                        >  2:674987/1‑225 (MQ=255)
                                         ccagACCCTAACTTGTGGTTGATTAAGGAAACCGAACTTTGGCCTAACCTGATTCCGGGCCTACATAAACGTAAAATTCCGCTGGTGATCGGTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTCGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCa                                                                                                                                                                                       <  2:290759/225‑1 (MQ=255)
                                             aCCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGGTGCGCTGGGCCCaaaaa                                                                                                                                                                                   >  2:68823/1‑228 (MQ=255)
                                                   aaCTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTCCATAAACGTAAAACTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGCAGAAGATGGTGCACGTTTTTTGGCGCTGGGCGCAAAAaataat                                                                                                                                                                              >  1:124005/1‑227 (MQ=255)
                                                        gtgtTGATTATGGAAACCGAACTATGGCCTAACCTTATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGCTATGCCAAACTGGGTAAATTCGTCCGGCGCTTGCTCCGTCGTCTTCCGCTGATTGCTGCGCAAAATGAAAAAAAGGGGGCACGTTTTGTGGCCCTCGGCGCAAAAAAAAAACAgggg                                                                                                                                                                        >  2:257444/1‑226 (MQ=255)
                                                        gtgtTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTgcgc                                                                                                                                                                                                                                >  1:423159/1‑172 (MQ=255)
                                                        gtgtTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGTCGCTGGGCGCCAAAAATAATCAGGTg                                                                                                                                                                        >  1:596288/1‑228 (MQ=255)
                                                        gtgtTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGGGCACGTTTTTTGGCGCTGGgcgc                                                                                                                                                                                        >  1:290759/1‑212 (MQ=255)
                                                         tgtTTATTTAGGAAACCCAACTAGGGCCTATCAGGATTGCGGCGCTGCATAAACGTAAAATTCCGGTGGTGATCGCGAACGCGCGACTCTCTTCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCCCGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGa                                                                                                                                                                       <  2:432925/228‑1 (MQ=255)
                                                         tgtTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCa                                                                                                                                                                                                                               <  2:423159/172‑1 (MQ=255)
                                                         tgtTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGGCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGAAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATAAGGTGa                                                                                                                                                                       >  1:370428/1‑228 (MQ=255)
                                                             gATTATGGAAACCGATCTTTGGCCTACCCTTATTGCGGCGCTTCATAACCGGAAAATTCCGCTTGTGGTCGCTAACGGGCGACTCTCTGCCCGCTTGGCCGCAGGTTATGCCCAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGt                                                                                                                                                                   <  2:366853/228‑1 (MQ=255)
                                                             gATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTTGCGTTGGGCGCCAAAAATAATAAGGTGAcc                                                                                                                                                                     >  1:702696/1‑226 (MQ=255)
                                                             gATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGCCGCAAAAAATAATCAGGTGACCGt                                                                                                                                                                   >  1:149398/1‑228 (MQ=255)
                                                             gATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGCCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGCAGATGTGGCCCGTTTTTTGGCGCTGGGCGCAAAAAAAAAACAGGGAACCGt                                                                                                                                                                   >  2:685852/1‑228 (MQ=255)
                                                                tATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTAc                                                                                                                                                                >  1:70236/1‑228 (MQ=255)
                                                                tATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGGGACCGTTAc                                                                                                                                                                >  1:162958/1‑228 (MQ=255)
                                                                  tGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGCCCTTCGCTTGCTGCGCCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGGGCACGTTTTGTGGCGCGGGGGGCCAAAAATAATCAGGTg                                                                                                                                                                        >  2:662823/1‑218 (MQ=255)
                                                                  tGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGGTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGAGCAAAATGAAGAAGATGGTGCACGTTTTTTGGCGCTGGGCGCACAAACATATCAGGGGAACGATACCg                                                                                                                                                              >  2:597333/1‑228 (MQ=255)
                                                                   ggAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGa                                                                                                                                                                       <  1:662823/218‑1 (MQ=255)
                                                                       aCCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGACGCAGATGGTGCCCGTTTTTTGGCGCTGGGCGCACAAAATAATCAGGGGACCGTTACCGGTAGc                                                                                                                                                         >  1:704124/1‑228 (MQ=255)
                                                                       aCCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTGGCTGCGTCGCATTACGCTAATTGCTGAGCAAAATGAAGAAGATGGTGCCCGTTTTGTGGCGCGGGGCGCAAAAAATAATAAGGAGACCGTTACCGGTAGc                                                                                                                                                         >  2:651473/1‑228 (MQ=255)
                                                                             cTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATGCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGaaa                                                                                                                                                   <  1:685852/228‑1 (MQ=255)
                                                                                tGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGACGAGGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTc                                                                                                                                                >  1:535599/1‑228 (MQ=255)
                                                                                tGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGTGCAAAAAATAATCAGGTGCCCGTTACCGGTACCCTGCAATTc                                                                                                                                                >  1:31163/1‑228 (MQ=255)
                                                                                       aCCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGGTTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGGGGCGCTGGGCGCAAAAAATAATAAGGGAGCCGGTACCGGTAGCCGGAAATTCGGTAttt                                                                                                                                         >  2:430720/1‑228 (MQ=255)
                                                                                       aCCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATATTCAGGTGACCGTTACCGGTAACCTGAAATTCGATCttt                                                                                                                                         >  1:34154/1‑228 (MQ=255)
                                                                                       aCCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATAtt                                                                                                                                          >  1:29754/1‑227 (MQ=255)
                                                                                            aTTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGGGCCCGTTTTTTGGCCCTGGGGGCAAAAAATAATAAGGTGACCGTACCCGGTAGCCGGAAATTCGTAATTTCttaa                                                                                                                                    >  2:594812/1‑225 (MQ=255)
                                                                                            aTTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGGGCACGTTTTGTGGCACTGGGCGCAGAAAATAATCAGGTGACCGTAACCGTTAGCCCTAAATTGGAATTTTCtata                                                                                                                                    >  2:182007/1‑228 (MQ=255)
                                                                                               ggggcgCTACATAAACGTAAAATTCCCGTGGGGACCGCTAACGCGCGATTCGCGGCTCGCTCGGTCGCAGGGTGTGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAGGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAAcg                                                                                                                                 <  2:535599/226‑1 (MQ=255)
                                                                                                 ggcgCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGGGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGAGACCGTTACCGGTAGCTCGAAATCCGATATTTCTGTAAccc                                                                                                                                >  1:439015/1‑225 (MQ=255)
                                                                                                 ggcgCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGGCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTTGGCGCAAAAAATAATAAGGTGACCGTTACCGGTAGCCTGAAATTCGAAATTTCTGTAAagc                                                                                                                                >  2:567186/1‑227 (MQ=255)
                                                                                                  gcgcTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAAACAGGGGACCGTTACCGGTAGCCTGAACTGCGAATTTTCTGTAAagccg                                                                                                                              >  2:112398/1‑228 (MQ=255)
                                                                                                  gcgcTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGGGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGGGACCGTTCCCGGTAACTTGAAATTCGAAAtttttaaaaaacca                                                                                                                              >  2:150753/1‑217 (MQ=255)
                                                                                                        taTAGACGTGAAATTCCCGTGTGGATCGTGCACGCGAGCCTCTCTGCCCGCGCGGCCGCAGGTTATGCCAAACTTGGTAAATTCTTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTg                                                                                                                        <  2:702696/227‑1 (MQ=255)
                                                                                                        cATAAACGTAAACTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTg                                                                                                                        <  1:676441/228‑1 (MQ=255)
                                                                                                        cATAAACGTAAAATCCCGCTGGTGATCCCTAACGCGCGACTCCCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTg                                                                                                                        <  2:70236/228‑1 (MQ=255)
                                                                                                                  aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAGTTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTGTGTGGCGCTGGGCGCAAAAAATAATCAGGGGACCGTTACCGGGAGCGCGAAAGTCGATATTTTTGTAAAGCCGCAGGTGGCTGCAAAAg                                                                                                              >  2:296168/1‑228 (MQ=255)
                                                                                                                  aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGTGGACCGTTACCGTTAGCCTGAAATTCGGTATTTCTGTAAGGCCCCAGTTGGCGGCTcaa                                                                                                               >  1:632368/1‑227 (MQ=255)
                                                                                                                  aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTaaa                                                                                                               >  1:518902/1‑227 (MQ=255)
                                                                                                                   aaTTCCCGTGGTGATCGGTAACGCGCGTCTCCCTGCCCGCTCGGCAGCAGGTTATGCCAAAACTGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTTGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTgg                                                                                                                       <  1:296168/218‑1 (MQ=255)
                                                                                                                                       acccgggaccttctGCCCGTCCGGGCGCAGGTTATGCCAAACGGGGTAAATTCGTCCGTCGCTTGCTGGGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGt                                                                                        <  1:667451/217‑1 (MQ=255)
                                                                                                                                       aCGGGCGACTCTTCGCCCGCTCGGCCGTCGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGt                                                                                        <  1:674987/228‑1 (MQ=255)
                                                                                                                                            cgtcccttcccccGCTGGGCCCGCGGTTATGCCAAACGGGGTAAATTCGCCCGTCGCTTGCTGGGTCGTATTGCGCTGATGGCGGCGAAAAATGAAGAAGGTGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTCCCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCATTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGca                                                                                   <  2:510934/218‑1 (MQ=255)
                                                                                                                                                 ttctGCCCGCCTGGCCGCAGGTTACGCCAAACTGGGAAAATTCGGCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACcaca                                                                              <  1:182007/227‑1 (MQ=255)
                                                                                                                                                 tacggCCCGCTCGGCCGCAGGTTAGGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACcaca                                                                              <  1:597333/224‑1 (MQ=255)
                                                                                                                                                 ccctGCCCGCCCGGCGCGAGGGTATGCCAAACGGGGTAAATTGGGCCATCGCTGGCTGCGTCGTATTACGCTGATGGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACcaca                                                                              <  2:396765/226‑1 (MQ=255)
                                                                                                                                                  tctGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCacac                                                                             >  1:43689/1‑228 (MQ=255)
                                                                                                                                                    tGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCg                                                                           >  1:631133/1‑228 (MQ=255)
                                                                                                                                                     cccccgtgggCCGATGGTTATGCCAAATGGGGATAGTGCGGCCGTGGCTTGCGGCGTCCTATTACTATGGTTGGTGGCAAAAATGAAGCAGAGGCTTCTCGTTTGGTGGCGCTGGGCGCAAAAAATAATCAGGTGGCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGc                                                                          <  2:149398/220‑1 (MQ=255)
                                                                                                                                                       cgccTCGGCCGCATGTTAGTCCAAACTGGGTAAATTCGGCCGTCGCTGGGTGCGTCGTTTTACGCTGATTGCGGCGCAAAATGAAGAAGTGGGGGCAGGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGGGACCGTTACCGGTAGCTTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGcc                                                                         <  2:694343/224‑1 (MQ=255)
                                                                                                                                                       ccGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGcc                                                                         <  1:493308/227‑1 (MQ=255)
                                                                                                                                                         gCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAATTAATCAGGTGACAGTTACCGGCAGCCTGCAATTCGATATTTCTGTAACGCCTCCGTTGGCTGCTAAAGCCCTGCCGCTGCGCCGCCCGTGGGCACCACcccgcccg                                                                       >  2:249095/1‑227 (MQ=255)
                                                                                                                                                          cTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGgcgc                                                                                                                                                                                        >  1:482643/1‑114 (MQ=255)
                                                                                                                                                          cTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACCCTGCCCCGCCAGTGGGCCCCCAACCCCCCgg                                                                      >  2:493308/1‑227 (MQ=255)
                                                                                                                                                          cTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTAGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCGGAAATTCGATATTTCTGTAACGCCGCAGCTGGCTGCTAAAGCCGTGACGCTGCGCGGCCAGTGGGCACCACACCGCCCGGt                                                                     >  1:212177/1‑228 (MQ=255)
                                                                                                                                                           tCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCa                                                                                                                                                                                       <  2:482643/114‑1 (MQ=255)
                                                                                                                                                            tggctgcgtgTTATGCCAAAGGGGGTAAAATCATGCGTGCGTTGTTGGGTCGTATTCGGGGGATGGGTGGCAAAAATGAAGAAGAGGGGGCACCTTTTGTGGTGGTGGGCGAAAAAAATAAGCAGGTGACCGTTACCGGTAGCTCGAAATCCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTAt                                                                   <  2:317765/219‑1 (MQ=255)
                                                                                                                                                            cGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGACGCTTGCCGCAAAAAATAATCAGGAGACCGGTACGGGTAGCCTCAAATTCGATATTTCTGTAACGACGCAGTTGGCGCCTAAAGCCGGAACGCTCCGCCCGCCAGAGCCACCACACCCCCCGGTac                                                                   >  2:165616/1‑227 (MQ=255)
                                                                                                                                                                 gggttgtATGCCAAACTGGGTAAATTCGTCCTTCTCTTGGTGCGTCGTATTACGCTTATTGCCGCCCAAAATGAAGAAGATGGTGGCGTTTTTGTGGCGATGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGAtt                                                              <  2:439015/222‑1 (MQ=255)
                                                                                                                                                                  cAGGTTATGGCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTCCGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGAtt                                                              <  1:132711/227‑1 (MQ=255)
                                                                                                                                                                         tGCCAAACTGGGGAAATTCATCCGGCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGGTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTACAGCCGTTACGCTGCGCCGCCAGTGGGAACCACACCGCACGGAATCGATTCCCACCAg                                                      >  1:538976/1‑228 (MQ=255)
                                                                                                                                                                               aCTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCCCCAGCCCTcc                                                >  1:592172/1‑227 (MQ=255)
                                                                                                                                                                               aCTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCa                                                >  1:668967/1‑228 (MQ=255)
                                                                                                                                                                               aCTGGGGAATTTTGCCCGTCGCTTGGCGCGTCTTTTTCCGGTGTTTGCGGCGAAAAAGGAAGAAGATGGTGCACGTTTTGTGGCGTTGGGCGCAAAAAATAATCAGGTGCCCGTGCCCGGTAGCCTGAATTTCGATATTTCTGTAACGCCGCAGTTGGCTCCTAAAGCCGTGCCGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCa                                                <  2:162958/228‑1 (MQ=255)
                                                                                                                                                                                 tGGGTAAATTTGTCCTTCCTTTGCCGCGTCGTATTACGCTGTGTGCTGCGAAAAATGCAGAAGGTGGTGCCCGGTTTGGGGGGCTGGGCGAAAAAAATAATCAGGTGACCGTTACCGGGAGCCTGAAATTCGATATCTCGGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACg                                              <  2:212177/228‑1 (MQ=255)
                                                                                                                                                                                 tGGGTAAATGTGTCGGTGGGTTGGTGGGTTGTTTTGCGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGCCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGGTATTTCTGTACCGCGGCAGTTGGCTGCTAAAGCCGTGACGATGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACg                                              <  2:124005/228‑1 (MQ=255)
                                                                                                                                                                                   ggaaaaTTCGTCCGTCGCTTGCTTCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCGGTTGTCTGCTCAAGCCGTGACGCTGAGCCGCCAGTGGGCACCACACCGCCCGGATAGGGTTGCCACCAGCACGCACCaa                                            >  1:689746/4‑228 (MQ=255)
                                                                                                                                                                                   ggTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGGTACCGGTAGCCTGCAATTCGATATTTCTGTAACGGCGCAGTTGGCGCGTAAAGCCGTGACGCTGCGCCGCCat                                                                                         >  2:150065/1‑182 (MQ=255)
                                                                                                                                                                                        aTTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGa                                       >  1:300585/1‑228 (MQ=255)
                                                                                                                                                                                        aTTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGa                                       <  1:68823/228‑1 (MQ=255)
                                                                                                                                                                                        aTTCGTCCGTCGCTTGCTGCGTCGTATTAAGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTAGATCTTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTTACGCTGCGCCGCCAGTGGGCCCAACACCGCCCTGTATGGCTTGCCACCAGCACTCACGAAGGCGa                                       >  1:611937/1‑228 (MQ=255)
                                                                                                                                                                                         gtCGTGCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGaa                                      <  1:512474/227‑1 (MQ=255)
                                                                                                                                                                                           cgtccggcGTTTGCTGCGTCGTATTGCGCTGATTGCTGCGAAAAATGAGGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAaga                                    <  1:158904/228‑1 (MQ=255)
                                                                                                                                                                                             tccgtcGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGAATGGATTGCCACCAGCACTCCCGAAGGCGAagag                                   >  1:341945/1‑227 (MQ=255)
                                                                                                                                                                                             gccgttGCGTGCTGCGTTGTATTACGCTGGTTGCTGCGAAAAATGAAGAAGATGGTGGACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagaga                                  <  1:310881/227‑1 (MQ=255)
                                                                                                                                                                                                             ttgttttACGTGGATTGCTGCGAAAAATGAAGAAGATGGTGCCCGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgccg                   <  2:43689/222‑1 (MQ=255)
                                                                                                                                                                                                                  ttACGCTGATTGCTGCGCAAAATGAAGATGATGGTGCCTGTGTTGTGGCGATGGGCGCAAAAAATAATCAGGGGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAt              <  1:504213/227‑1 (MQ=255)
                                                                                                                                                                                                                          aTTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATGGCCGCACATCAGGCAtt      >  1:140307/1‑227 (MQ=255)
                                                                                                                                                                                                                             gctgcGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCTGTGGGCACCACACGGCCCGGTATGGATTGCCCCCAGCCACCACGAAGGGGAAGAGAGGCTGGGGAGCGCCGCCCATCAGGCACTGttc  >  2:504213/1‑227 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                              
CGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTA  >  NC_000913/3808859‑3809307

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