Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 750,845 C→T S442N (AGC→AAC)  ybgQ ← putative fimbrial usher protein YbgQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913750,8450CT98.7% 246.4 / NA 76S442N (AGC→AAC) ybgQputative fimbrial usher protein YbgQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/41);  minor base A (1/0);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCAT  >  NC_000913/750598‑751088
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
gTCACCAATATCAACATTAAAAACTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGGAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGATGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAAAGACGATCTGAGAAGCGATAAGAGGCGAGCGAAATCGtgttgt                                                                                                                                                                                                                                                  >  1:497467/1‑251 (MQ=255)
     cAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGGGGAGGCGTCCACAATCCACAGGGCTTGATGTAGCAGGGTGGCGTACAGATTGAGGTTTAGTGGGGTAATCGGCTGGCCCACGGATAAGCTGATCGTCTGCTTTTCGTCCTGCGCCTCGCTGTCGTTGTATTTGTGATCCCGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGACATCGTGATGAtggtg                                                                                                                                                                                                                                             >  1:598431/1‑251 (MQ=255)
     cAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGAGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGGAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGACCCAAGGATAAGCTGATCGTCTGTTTTTCGACCTGCGCATCGCTGGAGTTGTATTTGTGATCCAGGTAGTTGGCGTAGCTGTGAACATGACGATCTGGGAAGCGATAAGCGGCGTGCGAAATCGTGTTGTtggtg                                                                                                                                                                                                                                             >  1:654315/1‑251 (MQ=255)
           cAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCAttc                                                                                                                                                                                                                                       <  2:225031/251‑3 (MQ=255)
            aaCATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCAttc                                                                                                                                                                                                                                       >  2:287446/1‑248 (MQ=255)
                  aaaaCCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGt                                                                                                                                                                                                                                 <  1:229488/250‑1 (MQ=255)
                      ccTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCt                                                                                                                                                                                                                            <  1:118383/251‑1 (MQ=255)
                              gATGTTGGCGGTCGTGGAGGCGTCTTCATTCCACCAGGTTTGATGTAGCAGGTTGGCTTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTaa                                                                                                                                                                                                                    <  1:79573/251‑1 (MQ=255)
                                       ggTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGc                                                                                                                                                                                                           <  1:327254/251‑1 (MQ=255)
                                           gTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTa                                                                                                                                                                                                        >  2:611902/1‑250 (MQ=255)
                                            tGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTa                                                                                                                                                                                                        >  1:515207/1‑249 (MQ=255)
                                             ggAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGc                                                                                                                                                                                                     >  1:216907/1‑251 (MQ=255)
                                               aGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGccc                                                                                                                                                                                                   >  1:632536/1‑251 (MQ=255)
                                                ggCGTATGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGACGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGTCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAACTCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAACACGGTAGCTTAAGcccc                                                                                                                                                                                                  >  2:268643/1‑251 (MQ=255)
                                                       gCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCg                                                                                                                                                                                             <  1:525306/249‑1 (MQ=255)
                                                        cATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGAGTGAGGTTTAGTGGGGTAATCGGTGGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTGGTTTTCTGTAATTAAAACGGTAGCTTAAGCTCCGCTCGTca                                                                                                                                                                                          <  2:616363/251‑2 (MQ=255)
                                                            ccaccaGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCctgct                                                                                                                                                                                      <  1:332366/251‑1 (MQ=255)
                                                              accaGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATGCACTTGTTTGCTGTAATTAAAACTGTAGCTTAAGCCCCGCTCGTCCTGCTGa                                                                                                                                                                                    <  1:74249/251‑1 (MQ=255)
                                                              aacaGGTTTGATGTAGCAGGATGGCGTAAAGATTGAGGTTTAGTAGGGTAATCGGATGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCctgctgg                                                                                                                                                                                    >  2:466719/3‑250 (MQ=255)
                                                                caGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGt                                                                                                                                                                                  >  1:402257/1‑251 (MQ=255)
                                                                     ttGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCa                                                                                                                                                                              >  1:338967/1‑250 (MQ=255)
                                                                        aTGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAAt                                                                                                                                                                          >  2:520928/1‑251 (MQ=255)
                                                                                   tGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCACAATGGCTACTGGcc                                                                                                                                                               >  2:484689/1‑251 (MQ=255)
                                                                                         aaaGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTg                                                                                                                                                         <  1:95763/251‑1 (MQ=255)
                                                                                                   gTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAgg                                                                                                                                               <  2:220694/251‑1 (MQ=255)
                                                                                                       aggggggTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGGTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAg                                                                                                                                           <  2:233028/248‑1 (MQ=255)
                                                                                                       aGTGGGGTAATCTGTTGGCCCACGGATAAGCTGGTCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAg                                                                                                                                           <  2:659753/251‑1 (MQ=255)
                                                                                                        gTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCTTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGt                                                                                                                                          <  2:417124/251‑1 (MQ=255)
                                                                                                        gTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGt                                                                                                                                          <  1:370172/251‑1 (MQ=255)
                                                                                                          ggggTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTgc                                                                                                                                        <  1:558012/251‑1 (MQ=255)
                                                                                                          ggggTAATCGGGTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACGTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTg                                                                                                                                         <  1:98582/250‑1 (MQ=255)
                                                                                                           gggTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTgcg                                                                                                                                       >  1:506142/1‑251 (MQ=255)
                                                                                                             gTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTCCTGAGTATCAAAATGGCTACTGGCCCAGGCGACATCAAAGGACAGTGCGcc                                                                                                                                     >  2:184364/1‑251 (MQ=255)
                                                                                                                aTCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAg                                                                                                                                  <  1:559096/251‑1 (MQ=255)
                                                                                                                   ggTTGGCCCACGGATAATCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGCCATCAAAGGACAGTGCGCCAAGccc                                                                                                                               >  2:95012/1‑250 (MQ=255)
                                                                                                                         cccACGGATAAGCTGAGCGTCTGTTGTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGAGCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAAGTAAAACGGTAGCTTAAGCCCCGCTCGGCCTGCTGAGTATTAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCat                                                                                                                         <  1:85180/251‑1 (MQ=255)
                                                                                                                              ggATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATAtttt                                                                                                                    <  1:611902/251‑1 (MQ=255)
                                                                                                                                  aaGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATGGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTTGCATTCACTTGTTTGCTGTAATTAAAACGGTTGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGcc                                                                                                                 <  2:546829/250‑1 (MQ=255)
                                                                                                                                        aTCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAc                                                                                                           >  1:349634/1‑250 (MQ=255)
                                                                                                                                          cGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAccc                                                                                                         <  1:184364/250‑1 (MQ=255)
                                                                                                                                                     cGTCCTGCGCATCGCTGTCGTTGTACTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGaa                                                                                             <  1:35827/251‑1 (MQ=255)
                                                                                                                                                     cGTCATGCGCATCGCTGTCGTTGTATGTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCCAGTGAAATCGTGTTGTTGGTGGCATTCCCTTGTTTGCTGTAATTAAACCGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAACGGCTACTGGCCCAGGTGACATCAAAGTACAGTGCCACAAGCCCCAGCATATTTTGCCCGATCCCCATTGCGGCAGa                                                                                              <  2:539761/250‑1 (MQ=255)
                                                                                                                                                        ccTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATgg                                                                                          >  1:12319/1‑251 (MQ=255)
                                                                                                                                                           gcgcATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTAACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGAGCCAAGCCACAGCATATTTTGCCCGATACCCATTGAGGCAGACTGGt                                                                                         >  1:287546/1‑249 (MQ=255)
                                                                                                                                                              cATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTATATGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAGTCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAACGGACAGTGCGCAAAGCCACAGCGTATTTTGCACGATACCCATTGCGGCAGAATGGTAGtc                                                                                     >  1:34878/1‑250 (MQ=255)
                                                                                                                                                                   cTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGtcatcatc                                                                               >  1:158160/1‑251 (MQ=255)
                                                                                                                                                                   cTGTCGTTGTATTTGTGATACAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGtcatcatc                                                                               >  1:354420/1‑251 (MQ=255)
                                                                                                                                                                     gTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGtcatcatc                                                                               >  2:54124/1‑249 (MQ=255)
                                                                                                                                                                          tGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATc                                                                        >  1:450454/1‑251 (MQ=255)
                                                                                                                                                                            tATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATc                                                                        >  1:204545/1‑249 (MQ=255)
                                                                                                                                                                             aTTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATcag                                                                      <  2:506142/250‑1 (MQ=255)
                                                                                                                                                                                   gATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGgccg                                                               <  2:402257/251‑1 (MQ=255)
                                                                                                                                                                                   gATCCAGATAGTGGGCGTAGCTGTGAAAATGACGATCTGAGAAGAGATAAGCGGCGAGCGAAATCGTGTGGTTGGTTGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTACTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGAGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGgccg                                                               <  2:505787/251‑1 (MQ=255)
                                                                                                                                                                                        aGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGt                                                           >  1:300293/1‑250 (MQ=255)
                                                                                                                                                                                                ggCGTAGCTGTGAAAAGGCAGATCTGAGAAGGGATAAGAGTATAGCTAAATCGTGTTGTTGGGGGCATTCTCTTGTTTGATGTAAATAAAGCGGTAGCTTAAGCCCCTCACGTCCGGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGTGCCAAGCCACAGGATATTTTGCCCGATACCCAGTGCGGCAGAATGGTAGTGATCATCAGAAAGGAGCAGGCTGGCGTACAGCGAgg                                                  <  2:59357/251‑1 (MQ=255)
                                                                                                                                                                                                  cGTAGCTGTGCCAATGACGATCTGAGCAGCGATAAGCGGCGAGAGAAATCGAGTTGTTGATGGCATTCACTTGTTTGCTGTAATTAAAAAGGGAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAACTGGATACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAGCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGAGAnggg                                                >  2:586441/1‑249 (MQ=255)
                                                                                                                                                                                                     aGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGttt                                             <  2:547591/251‑1 (MQ=255)
                                                                                                                                                                                                         gtgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCAAAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTg                                            >  1:6420/1‑248 (MQ=255)
                                                                                                                                                                                                          tgAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATAACCATTGCGGCCGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGATTgag                                          >  1:646453/1‑249 (MQ=255)
                                                                                                                                                                                                                tGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATccc                                  >  2:510817/1‑251 (MQ=255)
                                                                                                                                                                                                                 gACGATCTGAGAAGCGATAAGCGGCGAGTGAAATCGTGTTGTTGGTGGCATTCACTTGTGTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTTGCCCAGGTTACATCAAAGGACAGTGCGCCAAGCCACAGTATATTTTGCCCGATTCCCATTGCGGCAGAATGGTATTCATCATCAGAAATCAGCAGGCCGCTGTACAGCGAGGTGTTTGAGAGCATtccc                                 <  2:497467/251‑1 (MQ=255)
                                                                                                                                                                                                                 acgcacttccgatcggATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATcccc                                 <  2:562119/236‑1 (MQ=255)
                                                                                                                                                                                                                         gagaAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGa                            <  2:528062/248‑1 (MQ=255)
                                                                                                                                                                                                                         gagaAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAAc                         <  2:543957/251‑1 (MQ=255)
                                                                                                                                                                                                                         gagaAGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAAc                         <  2:355935/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaGCGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTc                      <  1:22317/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCTCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCAtt                   <  2:216907/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCAtt                   >  1:570540/1‑251 (MQ=255)
                                                                                                                                                                                                                                gATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATcccc                                 >  1:562119/1‑236 (MQ=255)
                                                                                                                                                                                                                                gATAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACGTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTg                  <  2:570960/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tAAGCGGCGAGCGAAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCt                >  1:538868/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gcgagcgaAATCGTTTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCTGGTGACATCAAAGGACAGTGCTCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAACCTTCATTGCTaaaaa           <  1:24843/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gagcgaAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTg                  <  1:355524/241‑1 (MQ=255)
                                                                                                                                                                                                                                          gagcgaAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTg                  >  2:355524/1‑241 (MQ=255)
                                                                                                                                                                                                                                           agcgaAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGt       <  2:570540/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cgaAATCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAACCGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGNttt     >  2:592509/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aaTCGTGTTGTTGGTGGCATTCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCAt  <  1:520928/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
GTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGTCGTTGTATTTGTGATCCAGATAGTTGGCGTAGCTGTGAAAATGACGATCTGAGAAGCGATAAGCGGCGAGCGAAATCGTGCTGTTGGTGGCATCCACTTGTTTGCTGTAATTAAAACGGTAGCTTAAGCCCCGCTCGTCCTGCTGAGTATCAAAATGGCTACTGGCCCAGGTGACATCAAAGGACAGTGCGCCAAGCCACAGCATATTTTGCCCGATACCCATTGCGGCAGAATGGTAGTCATCATCAGAAATCAGCAGGCCGCCGTACAGCGAGGTGTTTGAGAGCATCCCCCAGGAAACTTCATTGCTAAAAAAGGTTTCAT  >  NC_000913/750598‑751088

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