Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,577,406 T→C 4.5% Y185C (TAC→TGC)  mcrC ← 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,577,4060TC 4.5% 317.2 / 3.7 89Y185C (TAC→TGC) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base T (33/52);  new base C (2/2);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

TGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGA  >  NC_000913/4577118‑4577682
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                    
tGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTG‑AAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                                                                                                                                                                                                                                                                     <  2:424393/289‑1 (MQ=255)
 gATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                                                                                                                                                                                                                                                                     <  2:408841/289‑1 (MQ=255)
 gATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTa                                                                                                                                                                                                                                                                                     >  1:205871/1‑289 (MQ=255)
    aaaGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAgtg                                                                                                                                                                                                                                                                                  >  1:499498/1‑289 (MQ=255)
         ttttGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAGTCATAGAAACGGTAGTGTCCt                                                                                                                                                                                                                                                                              <  1:103186/288‑1 (MQ=255)
         ttttGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCt                                                                                                                                                                                                                                                                              <  1:157310/288‑1 (MQ=255)
            tGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGt                                                                                                                                                                                                                                                                          >  2:312290/1‑289 (MQ=255)
            tGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGt                                                                                                                                                                                                                                                                          >  1:340799/1‑289 (MQ=255)
                aaTGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTg                                                                                                                                                                                                                                                                      <  1:484999/289‑1 (MQ=255)
                  tGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGa                                                                                                                                                                                                                                                                     >  2:92587/1‑288 (MQ=255)
                  tGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAc                                                                                                                                                                                                                                                                    <  1:106544/289‑1 (MQ=255)
                  tGAAATGTTTCTGGTCCCATTTGTCGTGAAAAAATGCTCTTATAGTATTTGGGTGCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATTTCAGTTTCCATGCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGGGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTTTCCTTTGTTTTGAc                                                                                                                                                                                                                                                                    <  1:28998/289‑1 (MQ=255)
                   gAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTTGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAcc                                                                                                                                                                                                                                                                   <  1:318520/289‑1 (MQ=255)
                   gAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAcc                                                                                                                                                                                                                                                                   <  1:75946/289‑1 (MQ=255)
                   gAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGAcc                                                                                                                                                                                                                                                                   <  2:339038/289‑1 (MQ=255)
                        ttttCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGaa                                                                                                                                                                                                                                                              <  1:85183/289‑1 (MQ=255)
                        ttttCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGaa                                                                                                                                                                                                                                                              <  2:205871/289‑1 (MQ=255)
                                    aTTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTg                                                                                                                                                                                                                                                   >  2:173821/1‑288 (MQ=255)
                                          cgtGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGa                                                                                                                                                                                                                                            >  2:426958/1‑289 (MQ=255)
                                                 aaaTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCa                                                                                                                                                                                                                                     <  2:38132/289‑1 (MQ=255)
                                                 aaaTGCTCTTATAGTATTTGGCGGCAACGATAAATATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCa                                                                                                                                                                                                                                     <  1:198173/289‑1 (MQ=255)
                                                   aTGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGa                                                                                                                                                                                                                                   <  2:257028/289‑1 (MQ=255)
                                                             aGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAAc                                                                                                                                                                                                                         <  2:351910/289‑1 (MQ=255)
                                                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATACAGTAATGACATCTCTTTTTCGTTTCTTTCAAACTCATAGAAACGGTATTGTCCTTTGTTTTGACCTGGAATCGAATTATTGACGATGAATTTGCAGCCACTTATACCGact                                                                                                                                                                                                                     >  2:481117/1‑287 (MQ=255)
                                                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGaa                                                                                                                                                                                                                                           <  1:153161/267‑1 (MQ=255)
                                                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGaa                                                                                                                                                                                                                                           >  2:153161/1‑267 (MQ=255)
                                                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTg                                                                                                                                                                                                                                       >  2:434235/1‑271 (MQ=255)
                                                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTg                                                                                                                                                                                                                                       <  1:434235/271‑1 (MQ=255)
                                                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAAt                                                                                                                                                                                                                     >  2:387775/1‑289 (MQ=255)
                                                                 tttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATAACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAAt                                                                                                                                                                                                                     <  2:275749/289‑1 (MQ=255)
                                                                  ttGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATTGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAAtt                                                                                                                                                                                                                    <  1:116772/289‑1 (MQ=255)
                                                                   tGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAAttt                                                                                                                                                                                                                   >  2:437389/1‑289 (MQ=255)
                                                                   tGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAAttt                                                                                                                                                                                                                   >  2:466019/1‑289 (MQ=255)
                                                                        tCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTataat                                                                                                                                                                                                              >  1:162903/1‑289 (MQ=255)
                                                                         cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGa                                                                                                                                                                                                                                   >  2:517781/1‑267 (MQ=255)
                                                                         cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGa                                                                                                                                                                                                                                   <  1:517781/267‑1 (MQ=255)
                                                                         cAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTataat                                                                                                                                                                                                              <  2:218131/288‑1 (MQ=255)
                                                                                       ttttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                                                                                               <  1:466019/289‑1 (MQ=255)
                                                                                       ttttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                                                                                               <  2:340799/289‑1 (MQ=255)
                                                                                        tttCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTcc                                                                                                                                                                                              <  1:387775/289‑1 (MQ=255)
                                                                                             gatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                                                                                               >  2:40056/1‑283 (MQ=255)
                                                                                             gatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTc                                                                                                                                                                                               <  1:40056/283‑1 (MQ=255)
                                                                                             gatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCAt                                                                                                                                                                                         >  1:333620/1‑289 (MQ=255)
                                                                                             gatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGCAGTGTCCTTTGTTTTGAcc                                                                                                                                                                                                                                                                   <  2:347505/215‑1 (MQ=255)
                                                                                             gatgaGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGCAGTGTCCTTTGTTTTGAcc                                                                                                                                                                                                                                                                   >  1:347505/1‑215 (MQ=255)
                                                                                                    gAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                                                                                                   <  1:312290/288‑1 (MQ=255)
                                                                                                            aTGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAt                                                                                                                                                                          >  1:498457/1‑289 (MQ=255)
                                                                                                                cAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgctg                                                                                                                                                                      <  2:312905/289‑1 (MQ=255)
                                                                                                                               gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGaa                                                                                                                                                        <  2:499498/288‑1 (MQ=255)
                                                                                                                               gTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGaa                                                                                                                                                        <  1:173821/288‑1 (MQ=255)
                                                                                                                                tAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGaa                                                                                                                                                        <  2:498457/287‑1 (MQ=255)
                                                                                                                                tAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAag                                                                                                                                                       <  1:196702/288‑1 (MQ=255)
                                                                                                                                          aGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                                                                                                   <  1:131212/250‑1 (MQ=255)
                                                                                                                                          aGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGAt                                                                                                                                                                                   >  2:131212/1‑250 (MQ=255)
                                                                                                                                           gTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGc                                                                                                                                            >  1:363174/1‑288 (MQ=255)
                                                                                                                                              aCTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCgg                                                                                                                                        >  1:200982/1‑289 (MQ=255)
                                                                                                                                                            cTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaaa                                                                                                                          >  2:305037/1‑289 (MQ=255)
                                                                                                                                                             tCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAg                                                                                                                         >  2:391969/1‑289 (MQ=255)
                                                                                                                                                              cGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGa                                                                                                                                                                                                                       <  1:449006/194‑1 (MQ=255)
                                                                                                                                                              cGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGa                                                                                                                                                                                                                       >  2:449006/1‑194 (MQ=255)
                                                                                                                                                              cGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGt                                                                                                                        >  1:207725/1‑289 (MQ=255)
                                                                                                                                                                    aTCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACg                                                                                                                   <  1:426958/288‑1 (MQ=255)
                                                                                                                                                                    aTCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACg                                                                                                                   <  2:54881/288‑1 (MQ=255)
                                                                                                                                                                      cccATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACg                                                                                                                   <  2:267842/286‑1 (MQ=255)
                                                                                                                                                                      cccATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACg                                                                                                                   >  1:267842/1‑286 (MQ=255)
                                                                                                                                                                           tttAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAt                                                                                                           <  2:162903/289‑1 (MQ=255)
                                                                                                                                                                               aaTAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATctct                                                                                                       <  1:437389/289‑1 (MQ=255)
                                                                                                                                                                                             gtgtTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa                                                                                         <  2:324778/289‑1 (MQ=255)
                                                                                                                                                                                             gtgtTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa                                                                                         <  2:419752/289‑1 (MQ=255)
                                                                                                                                                                                                  tGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACttt                                                                                    <  2:333620/289‑1 (MQ=255)
                                                                                                                                                                                                        cggtAATTCACGACGGCAAAATTCATAAAGAAAATTTTGATAAAGTAATGACATGTCTTTTTCGTTTCTGTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGc                                                                              <  1:159877/286‑1 (MQ=255)
                                                                                                                                                                                                        cGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCAt                                                                                >  1:42086/1‑287 (MQ=255)
                                                                                                                                                                                                        cGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCAt                                                                                >  2:243576/1‑287 (MQ=255)
                                                                                                                                                                                                        cGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATg                                                                               <  2:200982/288‑1 (MQ=255)
                                                                                                                                                                                                          ttAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCtt                                                                            >  2:183314/1‑289 (MQ=255)
                                                                                                                                                                                                                   gacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAatat                                                                   <  1:92587/289‑1 (MQ=255)
                                                                                                                                                                                                                                                         gACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATACCGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCGCGGTCATTTTCTATAAAGTGATCGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCgtgt                             <  1:393947/289‑1 (MQ=255)
                                                                                                                                                                                                                                                         gACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCgtgt                             <  2:42086/289‑1 (MQ=255)
                                                                                                                                                                                                                                                               tctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAGTTGAACGACTTTAATCCCTGATGGTTGATTTTAACTTTCCTTGCTTATTTAATTTGCCAATTGTCCTTTTTAT‑ATTTGGTTAGCCAGCGTTTccttctt                      >  1:275690/1‑283 (MQ=255)
                                                                                                                                                                                                                                                               tctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTAAt                       >  1:405737/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               tctTTTTCGTTTCTTTCAAAATCATAGAAACGGCAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                                                                                                                                                                                                                                 <  2:135246/87‑1 (MQ=255)
                                                                                                                                                                                                                                                               tctTTTTCGTTTCTTTCAAAATCATAGAAACGGCAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGaca                                                                                                                                                                                                                                 >  1:135246/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                      cGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCat                <  1:183314/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAatat                                                                   <  2:299369/235‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAatat                                                                   >  1:299369/1‑235 (MQ=255)
                                                                                                                                                                                                                                                                           ctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCaa           <  2:363174/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  2:207725/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              ggTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTataat                                                                                                                                                                                                              <  1:78676/75‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              ggTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTataat                                                                                                                                                                                                              >  2:78676/1‑75 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                    
TGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTCTTATAGTATTTGGCGTCAACGATAAGTATTTTTTCTGATGAGCGAATGGTGATGTCAGTTTCCATTCGAGGTAACAAATTAAGTGACTGATCCGATATACTCGATGCATCCCATTTTAAATAAGAGCGGGTTGTGTTTGCAGACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGA  >  NC_000913/4577118‑4577682

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