Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,590,430 T→C intergenic (‑429/+424) ydeT ← / ← hipA fimbrial usher domain‑containing protein YdeT/serine/threonine‑protein kinase toxin HipA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,590,4300TC98.9% 320.4 / ‑3.9 87intergenic (‑429/+424)ydeT/hipAfimbrial usher domain‑containing protein YdeT/serine/threonine‑protein kinase toxin HipA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (47/39);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

CGAACAAGGGAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAGGGTTTATTAAAGCTGTCCGACAGACTTTCTTAACCTGGTAATATCTCTCTTGGGTGACTCTCCAAACATCCGCGAATATTCCCGACTGAAATGGGATAGCTTTCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAATTATGAA  >  NC_000913/1590203‑1590674
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gtaggcgatggatagtttgctggcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGGAACCGACagcatagctggct                                                                                                                                                                                                                               >  2:796416‑M2/198‑245 (MQ=255)
  aggcgatggatagtttgctggcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCtgt                                                                                                                                                                                                                             >  2:591358‑M2/196‑251 (MQ=255)
  aggcgatggatagtttgctggcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgatgcgtcgtttgagtaaacagtacgaacagataaacggtgattataatcaacctggctgtgagtaggtatagttcccatggccgCGAATATCCCCGACGGAAATGGGATAGCTTCCGGAACCGACAACATAAACTGCtgt                                                                                                                                                                                                                             <  2:315566‑M2/56‑1 (MQ=255)
            aagtttgctggcaaccccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCGGCTGTGGTGACATCg                                                                                                                                                                                                                   <  2:48730‑M2/66‑1 (MQ=255)
                      gcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGt                                                                                                                                                                                                                  <  2:25335‑M2/67‑1 (MQ=255)
                      gcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGt                                                                                                                                                                                                                  >  1:25335‑M2/176‑242 (MQ=255)
                       ggTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAGGGTTTATTAAAGCTGTCCGACAGACTTTCTTAACCTGGTAATATCTCTCTTGGGTGACTCTCCAAACATCCGCGAATATTCCCGACTGAAATGGGATAGCTTTCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcatt                                                                                                                                                                                                        <  2:652725/251‑1 (MQ=255)
                       caacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCTTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcatt                                                                                                                                                                                                        <  1:481122‑M2/77‑1 (MQ=255)
                       caacaccggttaacagaatgcgcttaccagaaagaaaacccatagctgtaatccttatttttgttgcttgttgtgcctgaaaatcaggctattcgtcgttttagtaaacagtacggacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTATCCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcatt                                                                                                                                                                                                        <  2:620219‑M2/77‑1 (MQ=255)
                        aacaccggttaccagaatgcgcttaccggaacgaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGCCAGCATAGGCTGCTGTGGTGACATCGTAATGttcattc                                                                                                                                                                                                       >  2:368471‑M2/174‑251 (MQ=255)
                        aacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcattc                                                                                                                                                                                                       <  1:793081‑M2/78‑1 (MQ=255)
                        aacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcattc                                                                                                                                                                                                       >  1:212498‑M2/174‑251 (MQ=255)
                             cggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCAt                                                                                                                                                                                                  <  1:19596‑M2/83‑1 (MQ=255)
                                   ccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCg                                                                                                                                                                                            >  1:295290‑M2/163‑251 (MQ=255)
                                           cgcttaccggaaagaaaacccatagctttaatccttattgtttatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCAGTCAGCATCTGTCGCCTGACtt                                                                                                                                                                                    <  2:934887‑M2/97‑1 (MQ=255)
                                                  cggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTg                                                                                                                                                                               >  1:320891‑M2/148‑249 (MQ=255)
                                                     aaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACg                                                                                                                                                                          >  2:945239‑M2/145‑251 (MQ=255)
                                                         aaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGc                                                                                                                                                                      >  1:837624‑M2/141‑251 (MQ=255)
                                                          aaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGcc                                                                                                                                                                     >  2:534921‑M2/140‑251 (MQ=255)
                                                           aacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGcc                                                                                                                                                                     >  1:453170‑M2/139‑250 (MQ=255)
                                                           aacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGcc                                                                                                                                                                     >  1:326312‑M2/139‑250 (MQ=255)
                                                            acccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCAt                                                                                                                                                                   >  2:214227‑M2/138‑251 (MQ=255)
                                                               catagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTc                                                                                                                                                                <  1:214227‑M2/117‑1 (MQ=255)
                                                                  agctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCTTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGa                                                                                                                                                             <  1:591358‑M2/120‑1 (MQ=255)
                                                                     tttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATAc                                                                                                                                                          >  2:409721‑M2/129‑251 (MQ=255)
                                                                       taatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCg                                                                                                                                                        >  2:220942‑M2/127‑251 (MQ=255)
                                                                          tccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCgg                                                                                                                                                       >  1:782665‑M2/124‑249 (MQ=255)
                                                                              tattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtatgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGCCGTAGCCATTTCTGATACCGGAGTgg                                                                                                                                                  >  2:85737‑M2/120‑250 (MQ=255)
                                                                              tattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTg                                                                                                                                                   <  2:33947‑M2/130‑1 (MQ=255)
                                                                                 tgttgatgaatgtagtgcctgacacacagccgattcgacgttttagtaaccagtacgaacagctaaacggttattataattaccctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGTACTc                                                                                                                                              >  2:888609‑M2/117‑251 (MQ=255)
                                                                                 tgtagatgctggttgtgcctgaaaatcaggcgtttcgtcgtttgagtaaccagtacgaacagataaacggttataataatcaacctggctgtgaggagctatagttgccaggtccgCGCATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACt                                                                                                                                               <  2:35810‑M2/134‑1 (MQ=255)
                                                                                        gcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgggagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTccc                                                                                                                                        <  2:657905‑M2/141‑1 (MQ=255)
                                                                                         cttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGc                                                                                                                                      >  2:758507‑M2/109‑251 (MQ=255)
                                                                                            gttgtgactgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCtt                                                                                                                                    <  1:796416‑M2/145‑1 (MQ=255)
                                                                                               gtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATc                                                                                                                                >  1:526177‑M2/103‑251 (MQ=255)
                                                                                                 gccttaaaatcaggcgattcgtcgttgtagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCg                                                                                                                              <  2:326312‑M2/151‑1 (MQ=255)
                                                                                                 gcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTc                                                                                                                                                                                   >  1:193672‑M2/101‑198 (MQ=255)
                                                                                                 gcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTc                                                                                                                                                                                   <  2:193672‑M2/98‑1 (MQ=255)
                                                                                                 gcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATcc                                                                                                                               <  1:220942‑M2/150‑1 (MQ=255)
                                                                                                     gaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTcc                                                                                                                           <  1:368471‑M2/154‑1 (MQ=255)
                                                                                                     gaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCt                                                                                                                          >  2:507774‑M2/97‑251 (MQ=255)
                                                                                                         atcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTaa                                                                                                                       >  2:214910‑M2/93‑250 (MQ=255)
                                                                                                           caggcgattcgtcgttttagaaaacagtacgaacagctaaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCCTTCCCGCTTGATCCGTCCTTAAAt                                                                                                                     >  2:482868‑M2/91‑250 (MQ=255)
                                                                                                               cgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgggagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTa                                                                                                                <  2:212498‑M2/165‑1 (MQ=255)
                                                                                                                     gtcgttttagtaaaccgtacgaacagatacacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCAGCTTGATCCGTCCTTAAATAAGTAACAGc                                                                                                           >  2:565878‑M2/81‑250 (MQ=255)
                                                                                                                       cgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCAc                                                                                                         <  2:837624‑M2/172‑1 (MQ=255)
                                                                                                                         ttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCt                                                                                                                                                                                                <  2:210811‑M2/85‑1 (MQ=255)
                                                                                                                         ttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCt                                                                                                                                                                                                >  1:210811‑M2/77‑161 (MQ=255)
                                                                                                                                 aacagtccgaacagataaacggttattctaatcaacctcgctgtgagtagctatagttgccagctccgCGAATATTCCCGACTGAAGTGGGATAGCTTCCGTAACCGACAGCTTAGGCTGCTGTGGTGACATCGTAATGTTCCTTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGAAACagcaccctctcgcgtt                                                                                              >  2:117109‑M2/69‑245 (MQ=255)
                                                                                                                                    agtacgaacagataaacggttattatactcaacctggctgtgactagctatagttgccaggtacgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATc                                                                                                                                                                                                 >  1:359235‑M2/66‑149 (MQ=255)
                                                                                                                                    agtacgaacagataaacggtgattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAAGGGGATAGCGTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATGCAGCATc                                                                                                                                                                                                 <  2:359235‑M2/84‑1 (MQ=255)
                                                                                                                                     gtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCg                                                                                          <  1:117109‑M2/187‑1 (MQ=255)
                                                                                                                                        cgaaccgataaacgcttcttctaatccacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTAt                                                                                       >  1:286228‑M2/62‑251 (MQ=255)
                                                                                                                                         gaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATg                                                                                      >  1:198576‑M2/61‑251 (MQ=255)
                                                                                                                                          aacagataaacggttattataatcagcctggctgtgagtagctatagttgccaggtccgCGAATATGCCCGACTGAAATGGGATAGCTTCCGTAACCTACAGCATAGGCTGCTGTGGTGTCATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCTCTTGTTCCGTCCTTAAATAAGTAACAGCACCCTCTTGTGTTTTCTTAtga                                                                                     <  2:426697‑M2/192‑2 (MQ=255)
                                                                                                                                           acagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCa                                                                                    <  2:198576‑M2/193‑1 (MQ=255)
                                                                                                                                                taaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTc                                                                                                                                              <  2:727843‑M2/135‑1 (MQ=255)
                                                                                                                                                taaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTc                                                                                                                                              >  1:727843‑M2/54‑188 (MQ=255)
                                                                                                                                                 aaacggttcgtataatcaacctggatgtgcgtagctatcgttgccaggcccgCGAATATTACCGACCGAAATGGGATAGCTTCCGTAACAGACAGCATAGGATGATGTGGTGACATCGTACTGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCt                                                                              >  2:866987‑M2/53‑251 (MQ=255)
                                                                                                                                                 aaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCATCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAAGAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTc                                                                               <  1:888609‑M2/198‑1 (MQ=255)
                                                                                                                                                              aatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGAGCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGacac                                                                 <  1:482868‑M2/212‑1 (MQ=255)
                                                                                                                                                               atcaaccagactgagaatagctatagttgcccggtccgCGAATATCCCCTACTGAAATGGGATAGCTTCCGTAACCGACAGAAAAGGCTGCTGCGGTGACATCGTAATGTTCATTCCGCATATGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACt                                                                >  2:62910‑M2/39‑251 (MQ=255)
                                                                                                                                                                  aacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTgcg                                                             >  1:491957‑M2/36‑251 (MQ=255)
                                                                                                                                                                    cctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGct                                                           <  1:214910‑M2/218‑1 (MQ=255)
                                                                                                                                                                      tggctgtgagtagctatagttgccaggtccgCGAATACTCCCGACTGCACTCGGATAGCTTCCGTAACCGACAGCCTAGGCTGCTGTGGCGAAATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGctct                                                         >  1:929041‑M2/32‑251 (MQ=255)
                                                                                                                                                                                 agctatagttgccaggtccgCGAATATTCCCGACGGAAATGGGATAGCTTCCGTAACCGACAGCAGAGGCTGCTGTGGTGACATCGTAATGTTCAGTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATGTCTGATGCCGGAGTGTACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCTCCCTCTGGTGTTTTCGGATTCAGTGTCTATGCTACAGACACTGCGCGCTTATCATGGGGc                                              <  2:929041‑M2/231‑1 (MQ=255)
                                                                                                                                                                                    tatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAgg                                           >  1:9376‑M2/18‑251 (MQ=255)
                                                                                                                                                                                          tgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATa                                      >  2:136218‑M2/12‑250 (MQ=255)
                                                                                                                                                                                           gccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAAc                                    <  2:453170‑M2/241‑1 (MQ=255)
                                                                                                                                                                                            ccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAAcc                                   <  1:945239‑M2/242‑1 (MQ=255)
                                                                                                                                                                                               ggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAgg                                >  2:697725‑M2/7‑251 (MQ=255)
                                                                                                                                                                                                  ccGCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCt                              >  2:55395/1‑250 (MQ=255)
                                                                                                                                                                                                  ccGCGAAGATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACGTCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGTACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCGCCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACTCTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTa                             <  1:466501/251‑1 (MQ=255)
                                                                                                                                                                                                   cgcgAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCTGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTgcgc                                                            <  2:908934/219‑1 (MQ=255)
                                                                                                                                                                                                   cgcgAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCTGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTgcgc                                                            >  1:908934/1‑219 (MQ=255)
                                                                                                                                                                                                   cgcgAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGCGGTGACATCATAATGTTCATTAAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTAt                            >  2:928618/1‑251 (MQ=255)
                                                                                                                                                                                                     cgAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTAttt                          <  2:491957/251‑1 (MQ=255)
                                                                                                                                                                                                      gAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTtc                         >  1:262102/1‑251 (MQ=255)
                                                                                                                                                                                                       aaTATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTtct                        <  1:241458/251‑1 (MQ=255)
                                                                                                                                                                                                          atTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTcc                      <  1:565878/250‑1 (MQ=255)
                                                                                                                                                                                                                 aCTGAAATGGGATAGATTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGcac                >  1:729219/1‑249 (MQ=255)
                                                                                                                                                                                                                       aTGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAAt        <  1:758507/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggATAGCTTCCGTAACCGACAGCCTAGGCTGCTGTGGTGCCATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCCGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACCGGCCAATTAt     >  2:549526/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAATTATg    >  2:365655/1‑250 (MQ=255)
                                                                                                                                                                                                                              aGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAAt        <  1:218178/244‑1 (MQ=255)
                                                                                                                                                                                                                              aGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAAt        >  2:218178/1‑244 (MQ=255)
                                                                                                                                                                                                                              aGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAATTATGaa  >  2:112109/1‑250 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CGAACAAGGGAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAGGGTTTATTAAAGCTGTCCGACAGACTTTCTTAACCTGGTAATATCTCTCTTGGGTGACTCTCCAAACATCCGCGAATATTCCCGACTGAAATGGGATAGCTTTCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAATTATGAA  >  NC_000913/1590203‑1590674

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