Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,175 T→A 7.9% intergenic (+21/+43) emrB → / ← luxS multidrug efflux system protein/S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,1750TA 7.9% 298.9 / 14.3 89intergenic (+21/+43)emrB/luxSmultidrug efflux system protein/S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base T (41/40);  new base A (2/5);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-03

Only 100 of 101 total aligned reads displayed.

ACTGAGTCGGTAAACCCGTTCAACCCGAATGCCCAGGCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCA  >  NC_000913/2813912‑2814446
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                               
aCTGAGTCGGTAAACCCGTTCAACCCGAATGCCCAGGCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtt                                                                                                                                                                                                                                                                      >  2:66470/1‑275 (MQ=255)
   gAGTCGGTAAACCCGTTCAACCCGAATGCCCAGGCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCttttt                                                                                                                                                                                                                                                                   <  2:482869/275‑1 (MQ=255)
   gAGTCGGGAAGCCCGGTCAACCCGAATGCCCAGGCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCttttt                                                                                                                                                                                                                                                                   <  2:46552/275‑1 (MQ=255)
                  ttCAACCCGAATGCCCAGGCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGa                                                                                                                                                                                                                                                     <  2:30695/274‑1 (MQ=255)
                           aaTGCCCAGGCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                                                                                                                                                                                           <  4:66469/275‑1 (MQ=255)
                            tagCCCAGGGGATTGTCAAGGCNNCTGGGNGGGCGTGGGAGGGGGCAACGGAGTGCATAAGGTGGGATTGCCCAGCAGTTCACCAGTCAGGGGCTGATTAGTTCCGCCAATGGGATCTTCTGGATGTAAGCTGGGATATTCCTCTTCCTGCGGGGGCTGGTGTGGTGTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                                                                                                                                                                                           <  3:3278/273‑1 (MQ=255)
                              gCCCAGGCGATGTACAGTCAACGGGAAGGGCTTGGGATGACGCAACAGCATGAATCATGCTGGATTGACCAGCAGATCACCAATCAGTGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTg                                                                                                                                                                                                                                        <  1:1699/275‑1 (MQ=255)
                                   ggCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACCAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAg                                                                                                                                                                                                                                   <  2:102289/275‑1 (MQ=255)
                                                       gAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACt                                                                                                                                                                                                                <  2:218978/275‑1 (MQ=255)
                                                                      aCGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTCATTGTGAAGATCGTTTACTGACTATATGTCCAGTTCCTGCCACTTCTCTTTCGGCAGt                                                                                                                                                                                                  >  2:335335/1‑274 (MQ=255)
                                                                      aCGCAACAGCAGGCACCAGGCTGGATGGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTGAGCCGTGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGTGGTGCGCCGTAAGTACAACTAAGCCAGGCGATTTGTACTGGCTTTTTGCCATTAATTGTGAAGATAGTTGACTGAATAGATGGGCCGTTCATGCAAATTCCCTTTCGGccggg                                                                                                                                                                                                 >  4:205516/1‑271 (MQ=255)
                                                                           aacgcagGCATCAGGCTGGAATGCCCAGCAGATCACCAAGCAGGGGCTGATTATTTCCGCCAAGGGGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCGGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAg                                                                                                                                                                                             <  4:480058/271‑1 (MQ=255)
                                                                                     caaggtgggATTGCCCAGCAGATGAGAAGTAAGGGGGTGATCAGTGCCGCAAATGGAAGCGTCTGGAGGCGAGCCGGGATATGCCGCGTCCGGCTGGGGCTGGTGGGGTTTGCTAAACCGCCATAGTGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATGTGAACTGGCTTTTGTCAATTAATTGTGAAGATAGTTTACTGACTAGGTGTGCAGTTCCTGCAAGTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                                                                                                                                                                   <  4:130576/267‑1 (MQ=255)
                                                                                                         gATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACg                                                                                                                                                               <  2:332126/274‑1 (MQ=255)
                                                                                                          aTCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCCGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTTCAGTTCCTGCAACTTCTCTTTCGGCAGTTCCAGTTCTTCGTTTCTGTTTATGCGTACTTCACGt                                                                                                                                                              >  4:214032/1‑274 (MQ=255)
                                                                                                            cACCAATCATTGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCATCCGGGATATCCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTACACCGCCATTTTGCGTAGGTGGCGGCGGAGGCGGTGCGCTCTAAGTACAACTACGCCAGTTCATTTGAACTTGCTTTTTCTAATTAATTGTGACGATATTTTCCTGCCTTTATGCGCATTTCCGGCAACTTTCCTTTCGGGCGTGCCATTTCTCGTTTGCTTTTGTTTCGTAACTCACTTTc                                                                                                                                                            >  4:294928/1‑274 (MQ=255)
                                                                                                            cACCAATCAGGGGCTGATTATTTCCGCCAATGAGACCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGGTCTTTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTCCTGTTTTTGTGTACGTccttttc                                                                                                                                                            >  4:402459/1‑268 (MQ=255)
                                                                                                                     tgggCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTg                                                                                                                                                                                                                                        <  1:145583/188‑1 (MQ=255)
                                                                                                                      gggCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGATGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGGNAGTTCGTTTGAACTGGCTTTTTTCAATTAATTTTGa                                                                                                                                                                                                                                                     >  3:139067/1‑175 (MQ=255)
                                                                                                                        gCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGa                                                                                                                                                                                                                                                                              <  6:921700/148‑1 (MQ=255)
                                                                                                                         cTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                                                                                                                                                 <  2:205517/272‑1 (MQ=255)
                                                                                                                         cTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGNTCATTTGAACTGGCTGTTTTCAATTAATTGTGAAGCTAGTTTACTTACTAGAGGTGCAGTTCATGCAACTTCCTTTTCTGCAGTGCTTGTTCTTCATTTCtttc                                                                                                                                                                              >  3:110435/1‑240 (MQ=255)
                                                                                                                         cTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGACGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                                                                                                                                                 <  1:40049/272‑1 (MQ=255)
                                                                                                                            attattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTg                                                                                                                                                                                                                                                                          <  5:780534/148‑1 (MQ=255)
                                                                                                                            attattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTg                                                                                                                                                                                                                                                                          <  8:286467/148‑1 (MQ=255)
                                                                                                                              tattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGa                                                                                                                                                                                                                                                                              >  6:808075/1‑142 (MQ=255)
                                                                                                                              tattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGa                                                                                                                                                                                                                                                                              <  5:808075/142‑1 (MQ=255)
                                                                                                                              tattTCCGCCAACGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGc                                                                                                                                                                                                                                                                        <  5:555565/148‑1 (MQ=255)
                                                                                                                                      ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTTCCAGTTCTTCGTTGCTGTTGATGCGTACGTCCCGTTCCAGAATGCTACGCGCAATATCCTgcgc                                                                                                                                 >  2:284189/1‑275 (MQ=255)
                                                                                                                                      ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcgc                                                                                                                                 >  2:191322/1‑275 (MQ=255)
                                                                                                                                      ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTCCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCCCAATATCCTgccc                                                                                                                                 >  4:363389/1‑273 (MQ=255)
                                                                                                                                           gagaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGa                                                                                                                                                                                                                                                                              >  6:150284/1‑129 (MQ=255)
                                                                                                                                           gagaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGa                                                                                                                                                                                                                                                                              <  5:150284/129‑1 (MQ=255)
                                                                                                                                           gagaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaa                                                                                                                                                                                                                                                           >  7:934486/1‑148 (MQ=255)
                                                                                                                                             gaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtt                                                                                                                                                                                                                                                                      >  8:569993/1‑135 (MQ=255)
                                                                                                                                             gaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtt                                                                                                                                                                                                                                                                      <  7:569993/135‑1 (MQ=255)
                                                                                                                                              aTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 <  8:169737/123‑1 (MQ=255)
                                                                                                                                              aTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGAGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 >  7:169737/1‑123 (MQ=255)
                                                                                                                                               tcttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgt                                                                                                                                                                                                                                                       >  8:1082427/1‑148 (MQ=255)
                                                                                                                                                   ctGGATGTCAGACGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGAGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAg                                                                                                                    <  2:131555/275‑1 (MQ=255)
                                                                                                                                                           cAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGCTAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGTTGCGTACGTCACTTTCCAGAATGCTACTCGCAATATCCTGCGCTTCCTGcca                                                                                                                        >  1:21941/1‑261 (MQ=255)
                                                                                                                                                                     tttCCTCGTGCTGCGGGGGGTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGGGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTGCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGGGCca                                                                                                  <  4:80484/274‑1 (MQ=255)
                                                                                                                                                                       tCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 >  5:780190/1‑98 (MQ=255)
                                                                                                                                                                       tCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 <  6:780190/98‑1 (MQ=255)
                                                                                                                                                                           cGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAg                                                                                                                                                                                                                           <  7:1082427/148‑1 (MQ=255)
                                                                                                                                                                            gTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGNTAGTTTACTGACTAGATGTGCAGNTCCTNCAACTTCTCTTTCGGCAGTGCCCGTTCTTCGtt                                                                                                                                                                                    >  4:540783/1‑186 (MQ=255)
                                                                                                                                                                                tgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTg                                                                                                                                                                                <  2:540784/186‑1 (MQ=255)
                                                                                                                                                                                 gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 >  8:141710/1‑88 (MQ=255)
                                                                                                                                                                                 gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 >  6:399155/1‑88 (MQ=255)
                                                                                                                                                                                 gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 <  7:141710/88‑1 (MQ=255)
                                                                                                                                                                                 gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                                                                                                                                                 <  5:399155/88‑1 (MQ=255)
                                                                                                                                                                                 gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTgg                                                                                                                                                                                                                                                                         >  6:1026718/1‑96 (MQ=255)
                                                                                                                                                                                 gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTgg                                                                                                                                                                                                                                                                         <  5:1026718/96‑1 (MQ=255)
                                                                                                                                                                                        cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttt                                                                                                                                                                                                                                                                    >  5:464876/1‑94 (MQ=255)
                                                                                                                                                                                        cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttt                                                                                                                                                                                                                                                                    <  6:464876/94‑1 (MQ=255)
                                                                                                                                                                                          ggTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCACCTTCTCTTTCGGCAGTGTCAGTTCTTCGTTGCTGTTGATGCGTAAGTTACGTTCCAGACTGCTACGCGCAATATCCTGCTCTTCCTCCATCGCCTGCATCTGGTCACTGCCACACTGGTAGACgttccgttcc                                                                             >  1:141244/1‑275 (MQ=255)
                                                                                                                                                                                                gtgtGCTAAACCGCCATTTGGCGCAGTGGGCGGCGGAGGCGGTGCCCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGTAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGAt                                                                        <  2:117730/271‑1 (MQ=255)
                                                                                                                                                                                                     cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTTCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTTAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATTCGTATGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGAtctgatt                                                                  >  2:486020/1‑275 (MQ=255)
                                                                                                                                                                                                     cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTAc                                                                                                                                                                     >  4:369621/1‑176 (MQ=255)
                                                                                                                                                                                                        aaCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTg                                                                                                                                                                                                                                        <  5:657232/106‑1 (MQ=255)
                                                                                                                                                                                                        aaCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTg                                                                                                                                                                                                                                        >  6:657232/1‑106 (MQ=255)
                                                                                                                                                                                                         aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCa                                                                                                                                                                 <  2:369622/176‑1 (MQ=255)
                                                                                                                                                                                                         aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTg                                                                                                                                                                                                                                                                          >  7:488778/1‑71 (MQ=39)
                                                                                                                                                                                                         aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTg                                                                                                                                                                                                                                                                          <  8:488778/71‑1 (MQ=39)
                                                                                                                                                                                                          ccGCCCTTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATCGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATTCGTACGTCACGTTCCAGTATGCTACGCGCAATTCCCTGCGCTTCCTGCAACGACTGCATCTGGTACGTGCCACCATGGTAAGCGTTAATTTCCGGATtccgattctgat                                                             >  4:182014/1‑275 (MQ=255)
                                                                                                                                                                                                              cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                                                                                                                                                         <  6:267260/131‑1 (MQ=255)
                                                                                                                                                                                                              cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                                                                                                                                                         >  5:267260/1‑131 (MQ=255)
                                                                                                                                                                                                                                      gCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCt                                                                                                                                                                                                                      >  6:1087018/1‑94 (MQ=255)
                                                                                                                                                                                                                                      gCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCt                                                                                                                                                                                                                      <  5:1087018/94‑1 (MQ=255)
                                                                                                                                                                                                                                           gcgcACTAAGTACAACTAAGCCAGTCCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtg                                                                                                                                                                                                                                >  2:568109/1‑79 (MQ=255)
                                                                                                                                                                                                                                               aCTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                                                                                                                                                                                            <  4:568108/79‑1 (MQ=255)
                                                                                                                                                                                                                                                cTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCa                        <  4:191321/274‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGa                                                                                                                                                                                                                                                     <  5:530200/51‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGa                                                                                                                                                                                                                                                     >  6:530200/1‑51 (MQ=255)
                                                                                                                                                                                                                                                   aGTACAACTAGGCCAGTTCATTTGAACGTGCTTTTTTCAATTAATTGTGAAGATAGTGTACTGACTAGTTGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGc                     <  2:294929/274‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTc                                                                                                                                                            >  2:94950/1‑129 (MQ=255)
                                                                                                                                                                                                                                                             aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTNCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCACCGCGTGCATCTGGTAAGGGCCACACTGGTAGACNTTCAGTTCCGGGATCTGATTCTGATCCTGCCCTTTCACCACGTCTNCNATTTCCGCTTTCCCTGCATCACCAacaa          >  3:48991/1‑274 (MQ=255)
                                                                                                                                                                                                                                                             aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGCATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGAGCTGATCCTTATCCTACACTTTCACCACGCCTCCCATTGCCCCTTTCCAGCCATCAacaaca           >  4:339606/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATc                                                                                                              >  3:80588/1‑173 (MQ=255)
                                                                                                                                                                                                                                                               ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGt                                                                                                      >  4:166193/1‑181 (MQ=255)
                                                                                                                                                                                                                                                               ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACgc        >  2:450121/1‑275 (MQ=255)
                                                                                                                                                                                                                                                               ccAGTTCATTTGAACTGGCTTCTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAGGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTAATCCTGCACTTTCAGCAGGCCTTCCATTGCCGCTTCCCAGGCATACGCaaacccc        >  4:136931/1‑270 (MQ=255)
                                                                                                                                                                                                                                                                 aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAg                                                                                                                                                                                                   >  8:378157/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                 aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAg                                                                                                                                                                                                   <  7:378157/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                                                                                                                                                  >  5:434508/1‑119 (MQ=255)
                                                                                                                                                                                                                                                                 aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTc                                                                                                                                                                  <  6:434508/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGa                                                                                                                                                        <  4:94949/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAGCTTCTCTTTCGGCAGTGCCAGGTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACgctgc     >  2:115857/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCg                                                                            >  4:292108/1‑204 (MQ=255)
                                                                                                                                                                                                                                                                  gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACgctgc     >  4:146530/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                                                                                                                                                                                              >  7:799535/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                   ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                                                                                                                                                                                              <  8:799535/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGt                                                                                                          <  1:80588/173‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAAATTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGTTTCCTGCAACGAGTGCATCTGGTAAGTGCca                                                                                                  <  2:166194/181‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCa  <  2:214033/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aTTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGAt                                                                        <  2:292109/204‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGa                                                                                                                                                                                                                                  >  7:810032/1‑48 (MQ=255)
                                                                                                                                                                                                                                                                        ttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGa                                                                                                                                                                                                                                  <  8:810032/48‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGtt                                                                                                                                                             <  5:161686/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGtt                                                                                                                                                             >  6:161686/1‑117 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                               
ACTGAGTCGGTAAACCCGTTCAACCCGAATGCCCAGGCGATGTACAGTCAACTGGAAGGGCTTGGGATGACGCAACAGCAGGCATCAGGCTGGATTGCCCAGCAGATCACCAATCAGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCA  >  NC_000913/2813912‑2814446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: