Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,085 C→A 9.5% 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,0850CA 9.5% 231.5 / 2.9 96Q183K (CAA→AAA) waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base C (55/30);  new base A (2/9);  total (57/39)
Fisher's exact test for biased strand distribution p-value = 6.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

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

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