Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,289,866 A→T 4.4% Q352L (CAG→CTG)  yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,289,8660AT 4.4% 329.5 / 2.2 89Q352L (CAG→CTG) yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base T (2/2);  total (46/43)
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.96e-01

GAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGC  >  NC_000913/3289580‑3290144
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
gaAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACCCCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCTAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAg                                                                                                                                                                                                                                                                                       <  2:129948/288‑1 (MQ=255)
 aaaGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAg                                                                                                                                                                                                                                                                                       >  2:512921/1‑287 (MQ=255)
      gcCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTAtt                                                                                                                                                                                                                                                                                <  1:512921/289‑1 (MQ=255)
        caccaGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTAttct                                                                                                                                                                                                                                                                              >  2:751279/1‑289 (MQ=255)
                  aGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTAccc                                                                                                                                                                                                                                                                     <  2:573388/288‑1 (MQ=255)
                         tttCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCa                                                                                                                                                                                                                                                             <  1:655080/289‑1 (MQ=255)
                         tttCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCa                                                                                                                                                                                                                                                             <  2:374799/289‑1 (MQ=255)
                                 tcccacTTTACTGGCGTGTAACTTGCTTCCGACGTTAATACAATGCCAAACAGCAAGCGCGTATTGGCCCCAACGGTAGGCTGTCTCGCAAACAGTAATGCACTCGTGACTATCAGCCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGCGGCCTTTGAGATTAACGATCTCTACCCCTCTGCCAACAGCAGCGATTTAGAAGTTACGATTGAAGAAAGTGACGGTACACAACCTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCTACCTg                                                                                                                                                                                                                                                     <  1:724709/284‑1 (MQ=255)
                                    acCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGc                                                                                                                                                                                                                                                  <  1:112860/289‑1 (MQ=255)
                                    acCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGc                                                                                                                                                                                                                                                  <  2:321970/289‑1 (MQ=255)
                                    acCTTTACTGGCGGGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGc                                                                                                                                                                                                                                                  <  2:333590/289‑1 (MQ=255)
                                      ctgtactgggtggcgattCGGTTCCGAAGGTAACAGGTTGCCAAACAGACAGCTCGAATTTGCCCAAACGATACGCGGTATCTCAACCAGGAGTGCATGCTTGACTATCAGGCAAAAGGGTTTTGTGATCGAGCAAAGCAACGTGTCAGCGGGTGCCTTTGAGTTTAACGATCTCTAGCCCTCTTCCTACAGCGGCGATTTAGAAGTCACGATTGAAAAAAGAGACGGTACACAACGTCGCTTTATACAGCCGTATTATTCATTACCCATGATGCAGCGACCTGGGCAt                                                                                                                                                                                                                                                <  1:25853/272‑1 (MQ=255)
                                        ttACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCt                                                                                                                                                                                                                                              <  2:291067/289‑1 (MQ=255)
                                          aCTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCACCAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTaa                                                                                                                                                                                                                                            >  1:409407/1‑289 (MQ=255)
                                          aCTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTGGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTaa                                                                                                                                                                                                                                            >  2:724709/1‑289 (MQ=255)
                                             ggCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGt                                                                                                                                                                                                                                         >  2:353160/1‑289 (MQ=255)
                                                gTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTgg                                                                                                                                                                                                                                                    <  1:415769/275‑1 (MQ=255)
                                                gTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTgg                                                                                                                                                                                                                                                    >  2:415769/1‑275 (MQ=255)
                                                   cAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAgc                                                                                                                                                                                                                                    >  1:358235/1‑288 (MQ=255)
                                                      cTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGa                                                                                                                                                                                                                                <  2:670898/289‑1 (MQ=255)
                                                           ttCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGAcc                                                                                                                                                                                                                              >  1:359825/1‑286 (MQ=255)
                                                              cgacgaTAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGa                                                                                                                                                                                                                        <  2:359825/289‑1 (MQ=255)
                                                                  gaTAATATGTTGCCAAACAGCCAGCGCGCATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCGACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCTGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGc                                                                                                                                                                                                                                                  <  2:498362/259‑1 (MQ=255)
                                                                        atGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATcgcg                                                                                                                                                                                                              <  2:332381/289‑1 (MQ=255)
                                                                            tGCCAATCAGCCAGCGGGGATTTGCACCAACGGTACGCGGTCTCGCAAACAGTAGTGCAATCGTGACGATCAGCCAACATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTCCCAACAGCGGCGCTCTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTAATACTGCCTTATTCTTCATTACCCATGATGCAGCGAGCTGGGCATCTAAAGtt                                                                                                                                                                                                                                        >  1:498362/1‑258 (MQ=255)
                                                                              ccAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATTCAGCGACCTGGGCATCTAAATTATAGCGCTACCGCTGGACGCTATccctcttat                                                                                                                                                                                                         >  1:413065/1‑280 (MQ=255)
                                                                                 aaCAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCaa                                                                                                                                                                                                     >  1:462793/1‑289 (MQ=255)
                                                                                  aCAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCaaa                                                                                                                                                                                                    <  2:298845/289‑1 (MQ=255)
                                                                                       cagcGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGt                                                                                                                                                                                                <  1:353160/288‑1 (MQ=255)
                                                                                                         aCGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAAt                                                                                                                                                                             >  2:577163/1‑289 (MQ=255)
                                                                                                                    tATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGc                                                                                                                                                                   >  1:77749/1‑288 (MQ=255)
                                                                                                                     aTCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCa                                                                                                                                                                 <  2:336189/289‑1 (MQ=255)
                                                                                                                       cGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACg                                                                                                                                                               <  2:358235/289‑1 (MQ=255)
                                                                                                                                tagtGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGACTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAatatat                                                                                                                                                      >  2:553592/1‑289 (MQ=255)
                                                                                                                                   tGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCGTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCAGGACGCTATCGCGCTGATGCAATCCGTGCTAGCAAGGGACCACAATTTGCTGAAGCCCAGGCAATATATGGt                                                                                                                                                   >  1:345777/1‑289 (MQ=255)
                                                                                                                                        tCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCa                                                                                                                                                                 <  2:539171/270‑1 (MQ=255)
                                                                                                                                        tCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCa                                                                                                                                                                 >  1:539171/1‑270 (MQ=255)
                                                                                                                                          gTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCa                                                                                                                                                                                         >  2:44591/1‑244 (MQ=255)
                                                                                                                                          gTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCa                                                                                                                                                                                         <  1:44591/244‑1 (MQ=255)
                                                                                                                                          gTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGaata                                                                                                                                            >  1:253300/1‑289 (MQ=255)
                                                                                                                                          gTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGaata                                                                                                                                            >  2:674611/1‑289 (MQ=255)
                                                                                                                                            gACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGaataa                                                                                                                                           <  2:413065/288‑1 (MQ=255)
                                                                                                                                            gACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGaataa                                                                                                                                           <  2:97768/288‑1 (MQ=255)
                                                                                                                                                  cAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTa                                                                                                                                                                                                                                             <  1:438418/184‑1 (MQ=255)
                                                                                                                                                  cAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTa                                                                                                                                                                                                                                             >  2:438418/1‑184 (MQ=255)
                                                                                                                                                   aGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTa                                                                                                                                   <  1:233875/289‑1 (MQ=255)
                                                                                                                                                                    gATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTg                                                                                                                  <  2:253300/289‑1 (MQ=255)
                                                                                                                                                                     aTCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTTTATGGCGGCCt                                                                                                                   >  1:7563/1‑287 (MQ=255)
                                                                                                                                                                      tCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCt                                                                                                                <  1:577163/289‑1 (MQ=255)
                                                                                                                                                                        tATCAAAGCAACGTGTCGGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCGAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGATTGATGGCGGccctcccg                                                                                                              >  1:418762/1‑283 (MQ=255)
                                                                                                                                                                             aaGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGt                                                                                                            >  1:612334/1‑286 (MQ=255)
                                                                                                                                                                             aaGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCt                                                                                                         <  2:77749/289‑1 (MQ=255)
                                                                                                                                                                              aGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTg                                                                                                        <  1:719442/289‑1 (MQ=255)
                                                                                                                                                                               gCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGa                                                                                                       >  1:56604/1‑289 (MQ=255)
                                                                                                                                                                                  aCGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACTTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGGTTTATGGCTGCCTGCCCGGTTCTGgaga                                                                                                    >  1:610438/1‑289 (MQ=255)
                                                                                                                                                                                        cAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATGTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGattatt                                                                                               <  2:462793/288‑1 (MQ=255)
                                                                                                                                                                                        cAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGattatt                                                                                               <  2:612334/288‑1 (MQ=255)
                                                                                                                                                                                           cGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATg                                                                                            >  1:335724/1‑288 (MQ=255)
                                                                                                                                                                                                            aaCGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCt                                                                                        <  2:53076/275‑1 (MQ=255)
                                                                                                                                                                                                            aaCGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCt                                                                                        >  1:53076/1‑275 (MQ=255)
                                                                                                                                                                                                            aaCGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCGTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATcg                                                                          <  1:751279/289‑1 (MQ=255)
                                                                                                                                                                                                                       ccccTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCAt                                                                                                                                                                                                                                                >  2:74944/1‑112 (MQ=255)
                                                                                                                                                                                                                       ccccTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCAt                                                                                                                                                                                                                                                <  1:74944/112‑1 (MQ=255)
                                                                                                                                                                                                                       ccccTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACtt                                                               <  2:409407/289‑1 (MQ=255)
                                                                                                                                                                                                                           tctTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAAGACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGgcg                                                           >  2:379274/1‑289 (MQ=255)
                                                                                                                                                                                                                            ctTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGgcg                                                           >  1:306756/1‑288 (MQ=255)
                                                                                                                                                                                                                             ttCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGgcgc                                                          >  1:542339/1‑288 (MQ=255)
                                                                                                                                                                                                                                     gcggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGgcgc                                                 >  2:350539/1‑289 (MQ=255)
                                                                                                                                                                                                                                      cggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGcgcg                                                >  1:91735/1‑289 (MQ=255)
                                                                                                                                                                                                                                       ggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCAtt                                                                                                                                                                                                                                                                         <  1:626988/71‑1 (MQ=255)
                                                                                                                                                                                                                                       ggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCAtt                                                                                                                                                                                                                                                                         >  2:626988/1‑71 (MQ=255)
                                                                                                                                                                                                                                       ggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGgc                                                   <  2:306756/285‑1 (MQ=255)
                                                                                                                                                                                                                                       ggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGcgcg                                                <  2:542339/288‑1 (MQ=255)
                                                                                                                                                                                                                                       ggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGt                                               >  1:428802/1‑289 (MQ=255)
                                                                                                                                                                                                                                       ggcgATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTAGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGt                                               >  2:118843/1‑289 (MQ=255)
                                                                                                                                                                                                                                           aTTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTc                                           >  1:777474/1‑289 (MQ=255)
                                                                                                                                                                                                                                             ttAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCTGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGa                                         <  2:629993/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGat                                    <  2:7563/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGat                                    <  2:34786/289‑1 (MQ=255)
                                                                                                                                                                                                                                                      aCGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCTGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGTGCATGTGGCCGGTTGTTGATGGTAATCa                                >  1:629993/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      aCGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCa                                >  2:277235/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      aCGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCa                                >  2:646909/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      aCGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCa                                >  2:313860/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                aCGTCGCTTTATCCATCCTTATTCTTCATTACCCAGGATGCAGCGGCCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAAc      <  1:350539/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATaa       >  1:219011/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                 cGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAAcc     >  2:316475/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                    cGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACCGAGCTGACACCCAATTCGATAACCAGc  >  2:699810/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                    cGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAg   >  2:456233/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                    cGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGc  >  1:498217/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
GAAAGCGCCACCAGCGGCAGTATTTTTTCCAGCTACACCTTTACTGGCGTGCAACTCGCTTCCGACGATAATATGTTGCCAAACAGCCAGCGCGGATTTGCCCCAACGGTACGCGGTATCGCAAACAGTAGTGCAATCGTGACTATCAGGCAAAATGGTTATGTGATCTATCAAAGCAACGTGTCAGCGGGTGCCTTTGAAATTAACGATCTCTACCCCTCTTCCAACAGCGGCGATTTAGAAGTCACGATTGAAGAAAGTGACGGTACACAACGTCGCTTTATCCAGCCTTATTCTTCATTACCCATGATGCAGCGACCTGGGCATCTAAAGTATAGCGCGACCGCTGGACGCTATCGCGCTGATGCAAACAGTGATAGCAAGGAACCCGAATTTGCTGAAGCCACGGCAATATATGGTTTGAATAATACTTTTACGCTGTATGGCGGCCTGCTCGGTTCTGAAGATTATTATGCGCTGGGGATCGGTATCGGCGGCACACTTGGCGCACTGGGCGCGTTGTCGATGGATATCAACAGAGCTGACACCCAATTCGATAACCAGC  >  NC_000913/3289580‑3290144

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