Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,939,335 A→T 4.4% V426D (GTC→GAC)  hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,939,3350AT 4.4% 329.3 / 2.9 91V426D (GTC→GAC) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base A (45/41);  new base T (2/2);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GATATCGAACTGGCAAGCTTTATGACGCCACCATTAACCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTA  >  NC_000913/3939048‑3939620
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                             
gATATCGAACTGGCAAGCTTTATGACGCCACCATTAACCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGAc                                                                                                                                                                                                                                                                                              <  2:550675/289‑1 (MQ=255)
   atCGAACTGGCAAGCTTTATGACGCCACCATTAACCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGACCATTGCTGCTGAATCACCAGTAATAATGCCCATTGTGATAAACGTAAAGTGGAATTGTTCGaaat                                                                                                                                                                                                                                                                                            >  1:514904/1‑285 (MQ=255)
        aCTGGCAAGCTTTATGACGCCACCATTAACCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTtg                                                                                                                                                                                                                                                                                      >  2:182286/1‑289 (MQ=255)
                           ccaccaTTAACCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATa                                                                                                                                                                                                                                                                   >  2:71600/1‑289 (MQ=255)
                            caccaTTAACCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATaa                                                                                                                                                                                                                                                                  >  1:688465/1‑289 (MQ=255)
                                    aCCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAg                                                                                                                                                                                                                                                          <  2:226737/289‑1 (MQ=255)
                                         tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTTATCAATAGGATAAACCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACg                                                                                                                                                                                                                                                      >  2:154942/1‑288 (MQ=255)
                                         tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACg                                                                                                                                                                                                                                                      >  2:9903/1‑288 (MQ=255)
                                         tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACg                                                                                                                                                                                                                                                      >  2:123066/1‑288 (MQ=255)
                                         tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACg                                                                                                                                                                                                                                                      >  1:616989/1‑288 (MQ=255)
                                         tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACg                                                                                                                                                                                                                                                      >  1:144158/1‑288 (MQ=255)
                                          aTCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGc                                                                                                                                                                                                                                                    <  1:97388/289‑1 (MQ=255)
                                           tCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGcc                                                                                                                                                                                                                                                   <  2:213765/289‑1 (MQ=255)
                                              accaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGc                                                                                                                                                                                                                                                 >  2:708319/1‑288 (MQ=255)
                                                caATCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCGTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCCGCTGAAGCAACAGTAATAATGCCCATGGAGATAAACGTATAGTGGAATTGTTAGACAGTCGTTGCGCCTTCCATTCCTTCAAAAACGGCAAGGACGGCCGCACt                                                                                                                                                                                                                                              >  1:604088/1‑289 (MQ=255)
                                                caACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACt                                                                                                                                                                                                                                              >  1:654673/1‑289 (MQ=255)
                                                caACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTAATCCGATTCTGAGGGCACGCGGTTCGGCTTAGAGATACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTTAGAAGCATGAACATTGCTGCGGAAGGAACAGTAATAATGCCCATTGTGATAAACGTATAATGGAATTGTTCGACAGTCGTTGTGTATTCCATTCCTTGATAAACGCGAGGGACGGCAGAACt                                                                                                                                                                                                                                              >  1:75993/1‑289 (MQ=255)
                                                       aaGCTGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGc                                                                                                                                                                                                                                        <  2:240711/288‑1 (MQ=255)
                                                       aaGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGc                                                                                                                                                                                                                                        <  2:466632/288‑1 (MQ=255)
                                                                            cGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCa                                                                                                                                                                                                                                       >  2:686713/1‑268 (MQ=255)
                                                                            cGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCa                                                                                                                                                                                                                                       <  1:686713/268‑1 (MQ=255)
                                                                                  aTGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGTTTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGTCAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTtgctgc                                                                                                                                                                                                            >  2:697763/1‑289 (MQ=255)
                                                                                      aCTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTg                                                                                                                                                                                                         <  2:466210/288‑1 (MQ=255)
                                                                                       cTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGAc                                                                                                                                                                                                       <  2:353410/289‑1 (MQ=255)
                                                                                       cTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGAc                                                                                                                                                                                                       <  2:59232/289‑1 (MQ=255)
                                                                                       cTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTAGTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGAc                                                                                                                                                                                                       <  1:154942/289‑1 (MQ=255)
                                                                                          atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACCGCAACGCCTAAACTAATCCACTGTTGCTGCGTGACCg                                                                                                                                                                                                     >  1:505668/1‑288 (MQ=255)
                                                                                          atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCg                                                                                                                                                                                                     >  2:503201/1‑288 (MQ=255)
                                                                                            ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCa                                                                                                                                                                                                  >  2:274701/1‑289 (MQ=255)
                                                                                               tcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGaa                                                                                                                                                                                               <  2:616989/289‑1 (MQ=255)
                                                                                               tcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGaa                                                                                                                                                                                               >  1:330795/1‑289 (MQ=255)
                                                                                                cGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAac                                                                                                                                                                                              <  2:516285/289‑1 (MQ=255)
                                                                                                       cccAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGc                                                                                                                                                                                                                                        >  2:546390/1‑240 (MQ=255)
                                                                                                       cccAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGc                                                                                                                                                                                                                                        <  1:546390/240‑1 (MQ=255)
                                                                                                          aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCg                                                                                                                                                                                    <  2:688465/289‑1 (MQ=255)
                                                                                                            ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTAcc                                                                                                                                                                                     >  2:153586/1‑286 (MQ=255)
                                                                                                             cGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCgctg                                                                                                                                                                                 <  1:9903/289‑1 (MQ=255)
                                                                                                             cGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCgctg                                                                                                                                                                                 <  2:144158/289‑1 (MQ=255)
                                                                                                                         aacaacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCa                                                                                                                                                                     <  2:32130/289‑1 (MQ=255)
                                                                                                                         aacaacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCa                                                                                                                                                                     <  2:88116/289‑1 (MQ=255)
                                                                                                                                 ttACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAg                                                                                                                                                             <  2:330795/289‑1 (MQ=255)
                                                                                                                                 ttACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAg                                                                                                                                                             <  1:677608/289‑1 (MQ=255)
                                                                                                                                       cTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGc                                                                                                                                                       <  1:503201/289‑1 (MQ=255)
                                                                                                                                       attAGTCCGATGGTGCTGGAACGCGGTTCGGCTTAGATTTACGCTGTTTTTTGCTCAAATTATTTCCATCGGTAGTTTTCAGAAGCAGGAAAATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATATACGTGTAGCTGAAGTGTTCGTCAGTCGTTGTGCCTTCCATTCCTTCGTAAACGCGAAGGACGGCCGCCCTTAAAGTGACTCCTGAACTAATAGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGACAGATCGGc                                                                                                                                                       <  1:697763/288‑1 (MQ=255)
                                                                                                                                          aCTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCaa                                                                                                                                                     >  2:16682/1‑288 (MQ=255)
                                                                                                                                          aCTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCaa                                                                                                                                                     <  1:16682/288‑1 (MQ=255)
                                                                                                                                            tCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTTTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAgtg                                                                                                                                                  <  1:708319/289‑1 (MQ=255)
                                                                                                                                            tCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAgtg                                                                                                                                                  <  2:159808/289‑1 (MQ=255)
                                                                                                                                            tCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAgtg                                                                                                                                                  <  2:72379/289‑1 (MQ=255)
                                                                                                                                            tCCGATTCTGATGGAAAGCCGGTCAGCTTAGATTTACGCGGTCTTTTTGTCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACAGTGCTGCTGGAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGGTCGACAGTCGTTGTGCCTTCCATTCCTTCAGAAACTCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGATGCGTGACCGCCAGAACACTGTTACCGCTGCTGGGATTGTCATCGGTCAGATCGGCAAgtg                                                                                                                                                  <  1:501130/289‑1 (MQ=255)
                                                                                                                                                      aTGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTAt                                                                                                                                         >  2:690150/1‑288 (MQ=255)
                                                                                                                                                        ggAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCa                                                                                                                                      <  1:153586/289‑1 (MQ=255)
                                                                                                                                                                         aGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGt                                                                                                                      >  1:313036/1‑288 (MQ=255)
                                                                                                                                                                               aCGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACAt                                                                                                                >  1:403433/1‑288 (MQ=255)
                                                                                                                                                                                cGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCATTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTACACTAATCGACAGTTGCTGCGTGACTGCCAGAACACTGTTACCGCTGCTGGCATTTCTCTCGGTCAGTTCGGCAAGTGGGATGGTATTCATCGAGGTAAATTGCGTCAACATa                                                                                                               >  1:116852/1‑288 (MQ=255)
                                                                                                                                                                                 gCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGc                                                                                                             <  2:313036/289‑1 (MQ=255)
                                                                                                                                                                                      cTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATcc                                                                                                        <  2:505668/289‑1 (MQ=255)
                                                                                                                                                                                       ttttGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGTTGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATccc                                                                                                       <  1:182286/289‑1 (MQ=255)
                                                                                                                                                                                                        ccATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCa                                                                                      <  1:71600/289‑1 (MQ=255)
                                                                                                                                                                                                               ttgttTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGTCAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCa                                                                                >  2:522084/1‑288 (MQ=255)
                                                                                                                                                                                                               ttgttTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTATATAAACAACGGCAAGATCAGc                                                                              >  2:748109/1‑289 (MQ=255)
                                                                                                                                                                                                                tgttTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGc                                                                              <  2:654673/289‑1 (MQ=255)
                                                                                                                                                                                                                gttgttCAGAAGCAAGAACATTGGTGCTGAAGCAACAGTAATAATGCCCATTTGTATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGGCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGc                                                                              <  1:711336/285‑1 (MQ=255)
                                                                                                                                                                                                                 tgtttCAGAAGAATGAGCATTGCTGCGTAAGCAACAGTAATAAAGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCa                                                                             <  1:244533/287‑1 (MQ=255)
                                                                                                                                                                                                                             aTGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGc                                                                   <  2:227960/287‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGATGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAgtg                                                                                                                                                  <  1:382890/201‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGATGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAgtg                                                                                                                                                  >  2:382890/1‑201 (MQ=255)
                                                                                                                                                                                                                                     tgctgcTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCg                                                                                                                                                                                                                                                  <  1:504303/104‑1 (MQ=255)
                                                                                                                                                                                                                                     tgctgcTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCg                                                                                                                                                                                                                                                  >  2:504303/1‑104 (MQ=255)
                                                                                                                                                                                                                                       ctgctgAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCCTCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTAGACATAGCCATCCCTAATACAAACAACGGCCAGATCAGCATCCATATAGCCATTGCTGGTGa                                                       >  2:377489/1‑289 (MQ=255)
                                                                                                                                                                                                                                           tgAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCAGAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGc                                                   <  2:514904/289‑1 (MQ=255)
                                                                                                                                                                                                                                             aaGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCaa                                                 >  1:160243/1‑289 (MQ=255)
                                                                                                                                                                                                                                                    aGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGaa                                           >  1:83487/1‑288 (MQ=255)
                                                                                                                                                                                                                                                         taatGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGc                                     >  2:640333/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               ccATTTTGATAAACGTGTAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCAt                               <  2:604088/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                cATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCTATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCAtt                              >  1:567486/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                cATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCAtt                              >  2:207498/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                cATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCACCGCGGTAAACTGCGTCGAACTAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGCGACTGCAACGAGAATTGAGCTATCAtt                              >  1:96955/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                cATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGGGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCAtt                              >  2:646688/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  ttGTTATAGACGGATAGTGTGACTGTTATACAATCTTTGTGCCTTCCAATCCTGCATACACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAAGCGACAGTGGCTGCGTGACCGCCAGAACACTGTTACCGTTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAg                            <  1:174037/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ataGTGGAATTGTTCGACAGACGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGAAGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGCGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTCCTGGTAATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCAc                >  2:136120/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                 aGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGt             >  2:86283/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   tGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGa           >  2:235550/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                   tGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACAGTGa           >  1:540924/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                    ggATTTGTTCGCCAGTAACTGTGCCTTCCTTTCCCACTTAAACTCTAAGGACGGCCAAACTTACAGCAACCCCTAAACTAACCGAAGTTTGCTGCGTGACGGCCAGAACACTGTTACCGCTGCTGGCATTGTAATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCTTCGACATAGCCATCCCTATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGAt          <  1:748109/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     gAATTGTTCGTCAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGAt          <  1:522084/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      aaTTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATc         >  2:74414/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                          gTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATcccc      >  2:242354/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                                           ttCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACt   >  1:515916/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           ttCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATACCCATTGCTTGTGATTGCAAAGAGAACTGAGCGATCATTAGACCAATACTCCACGTGTACACCACt   >  2:726734/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                             cGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTa  >  2:287639/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                             cGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTa  >  2:421786/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                             
GATATCGAACTGGCAAGCTTTATGACGCCACCATTAACCACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTTCCATTCCTTCATAAACGCGAAGGACGGCCGCACTTACAGCAACGCCTAAACTAATCGACAGTTGCTGCGTGACCGCCAGAACACTGTTACCGCTGCTGGCATTGTCATCGGTCAGATCGGCAAGTGTGATGGTATTCATCGCGGTAAACTGCGTCGACATAGCCATCCCTAATATAAACAACGGCAAGATCAGCATCCATATAGCCATTGCTGGTGATTGCAAAGAGAACTGAGCGATCATTAGCCCAATAATCACCGTGATCCCCACTA  >  NC_000913/3939048‑3939620

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