Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,256,157 C→T Q165Q (CAG→CAA cyuA ← putative L‑cysteine desulfidase CyuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,1570CT97.6% 277.6 / ‑1.3 83coding (495/1311 nt)cyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/42);  minor base . (0/2);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATT  >  NC_000913/3255911‑3256390
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
ggATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTggg                                                                                                                                                                                                                                         >  2:406547/1‑251 (MQ=255)
ggATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCCCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTGCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTgg                                                                                                                                                                                                                                           <  2:175558/250‑1 (MQ=255)
    gAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTaa                                                                                                                                                                                                                                     <  1:308639/250‑1 (MQ=255)
         aagaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGACGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGGAAacacc                                                                                                                                                                                                                                >  2:56569/2‑250 (MQ=255)
          agaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCCCACACGCCTGTTGGGTAAacacccc                                                                                                                                                                                                                              >  2:237446/1‑249 (MQ=255)
           gaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGACTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAacaccac                                                                                                                                                                                                                              >  1:375214/1‑250 (MQ=255)
             tCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACCCACGCCTGTTGGGTAAacaccacacc                                                                                                                                                                                                                           >  2:696092/1‑251 (MQ=255)
                    ccAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGc                                                                                                                                                                                                                    <  1:32614/251‑1 (MQ=255)
                     cAGCAAACCGCGCTCGCCCTGTTTTTCCAGCGTCGCGCCAATATGCCGCCCCCACTTACCGCTCAACCCTTCCTGCGATAACGCACACTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCAGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCg                                                                                                                                                                                                                   >  2:307396/1‑251 (MQ=255)
                     cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCg                                                                                                                                                                                                                   >  1:468616/1‑251 (MQ=255)
                     cAGCAAACCGCCCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGACCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCGCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCg                                                                                                                                                                                                                   <  2:409589/251‑1 (MQ=255)
                       gCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACGTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGtc                                                                                                                                                                                                                 <  2:164134/251‑1 (MQ=255)
                        cAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGtct                                                                                                                                                                                                                <  1:305556/251‑1 (MQ=255)
                         aaaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGtctc                                                                                                                                                                                                               >  2:78328/1‑251 (MQ=255)
                         aaaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGtctc                                                                                                                                                                                                               <  1:570876/251‑1 (MQ=255)
                               cgcTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCCCACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATg                                                                                                                                                                                                         >  1:697109/1‑251 (MQ=255)
                                 cTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGc                                                                                                                                                                                                        >  1:314812/1‑250 (MQ=255)
                                 atCGCAATGTTTATCCATCGTCGCCCCAATATGCAGCCCCCACTTAACGCTCACTCCTTCCTGCGCTAAGGCTCAATTTAGCTTCGCGGAATCGAGCATAAAGCGGATCGCCGCAAACGGGACTTCATTGCCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCCCCGTCAGCGGAGAATCTTGTTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGCGCGTCTCGTTATGc                                                                                                                                                                                                        >  2:499479/2‑250 (MQ=255)
                                           ttttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACATTCCTCCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGAGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTAAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGt                                                                                                                                                                                             >  1:707602/1‑251 (MQ=255)
                                              tCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                                                                                                                                                                                           <  2:617634/250‑1 (MQ=255)
                                                cAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                                                                                                                                                                                        >  1:525016/1‑251 (MQ=255)
                                                 aGCGTCGCGCCAATATGCAGCCCCCACTTACCGCCCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGccc                                                                                                                                                                                       <  2:602570/251‑1 (MQ=255)
                                                       gcgcCAAAATGCAGCCCCCACTTACCCCTAAAACCTTCCTGCGATAACGCACAATTTAGCTGCGCGGAAGCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTGCGGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCGGCCACCCACGCCTGTTGGGTAAACCCCACCCCATCGTGCGTCTCGATATGCCCAATGTTGGTATGCCCGCcgac                                                                                                                                                                                 <  2:517887/251‑1 (MQ=255)
                                                           cAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATg                                                                                                                                                                             <  1:237446/251‑1 (MQ=255)
                                                             atatGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCTGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGt                                                                                                                                                                           <  2:691050/251‑1 (MQ=255)
                                                               atGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTg                                                                                                                                                                          >  1:106095/1‑250 (MQ=255)
                                                               atGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTCCATTGACGATCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGa                                                                                                                                                                         <  2:664756/251‑1 (MQ=255)
                                                                 gCAGCCCCCACTTCCCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACGTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACGCTTGCTCGCCCTCTGCCACACCCGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGTGaca                                                                                                                                                                       <  2:469696/251‑1 (MQ=255)
                                                                       cccACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                                                                                                                                                                                           <  2:110706/225‑1 (MQ=255)
                                                                       cccACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                                                                                                                                                                                           >  1:110706/1‑225 (MQ=255)
                                                                         cACTTACCGCTCAACCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCa                                                                                                                                                               <  2:535743/251‑1 (MQ=255)
                                                                                   tCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGCGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCg                                                                                                                                                     >  1:103298/1‑251 (MQ=255)
                                                                                   ccAAACCTCCCTGCAATAACGCGCAAGTTAGCTGCGCGGAATCGAGAGTAAAGCGGATCGCCTCAAACGGGACCGCAGTGACGAACTTCAGGATCGCAGCCAGCGTTGTTGTGGAACGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTCCCACACACGCCTGGTGGGTAATCTCTACTCCGTCGGTCGTCTCGATATGTACATTGTGGGGATGCCCGCCGTTGATGGTGGCACACGCCCACTTCTCAcc                                                                                                                                                      <  2:70110/249‑1 (MQ=255)
                                                                                     aacccTTCCTGCGATAACGCACAATGTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTGCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTGGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGtt                                                                                                                                                   <  2:478747/248‑1 (MQ=255)
                                                                                     aaaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGtt                                                                                                                                                   >  1:297930/1‑251 (MQ=255)
                                                                                                cGATAACGCACAATTTAGCTGCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACTCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTgg                                                                                                                                        <  1:434204/251‑1 (MQ=255)
                                                                                                    aaCGCACAATATAGCTTCGCGAAAGCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCGCAGCCATCGTCGGTCTGGAAAGCACCTTCAGCGGAGACTCTTGCGCGCCCGCTGCCACACACGCCTGTGGTGTAAACACCACACCATCGTTGGTCTCGATATTCGCAGTGGTGGTATGCCCGTCGAGGATGTTGACACACGCCCCCGGCTCTCCGTTCCAGACTTTGGcgcg                                                                                                                                    <  2:40845/251‑1 (MQ=255)
                                                                                                      cGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                                  >  1:405635/1‑251 (MQ=255)
                                                                                                         acaATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                                  >  1:196338/1‑248 (MQ=255)
                                                                                                            aTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTgaga                                                                                                                               >  2:520057/1‑248 (MQ=255)
                                                                                                            aTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGagaaga                                                                                                                            <  2:205995/251‑1 (MQ=255)
                                                                                                                  cTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGAttt                                                                                                                      >  1:121851/1‑251 (MQ=255)
                                                                                                                       cgtaaTCGAGAATACAGCGGATCGCCGCAAACCGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGCGACTCTTGCTCGCCCTCTGCCACACACCCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGCCACACGCCCACTTCTCACCGTTCCAGCCTTTGACGCGTGCGAAGTGGATTTCATCg                                                                                                                 >  2:506083/4‑251 (MQ=255)
                                                                                                                           aTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAg                                                                                                             <  1:620485/251‑1 (MQ=255)
                                                                                                                                  aTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTg                                                                                                      >  1:247196/1‑251 (MQ=255)
                                                                                                                                      aGCGGCTCGCCGCAAACGGGACTTCATTGCCCACCTTCAGGATCTCAGCCAGCGTCGTTCTGGCACGCACCGTCAGCGGAGACTCTTGCACGACCTCTGCCCCACCCGCCTGTTGGGAAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGAAAGGTTCCTGGATc                                                                                                  >  2:430394/1‑251 (MQ=255)
                                                                                                                                       gCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACCCTTGCTCGCCCCCTGCCACACACGCCTGTGGGGTAAACACCACACCATCGTGCGTCTCGATAGGCACAATGTTGGTATGCCCGTCGACGATGGTGACACACGCCCACATCTCACCGTTCCAGTCGTTGGCGCTTGAGAAGAGGATGTCATCGCAAGGTGCCTTGATCt                                                                                                 <  2:723749/251‑1 (MQ=255)
                                                                                                                                           aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAAc                                                                                             >  2:250329/1‑251 (MQ=255)
                                                                                                                                           aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAAc                                                                                             <  1:373368/251‑1 (MQ=255)
                                                                                                                                           aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAAc                                                                                             >  1:415693/1‑251 (MQ=255)
                                                                                                                                                   aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGAGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACtt                                                                                     >  1:365985/1‑251 (MQ=255)
                                                                                                                                                    aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACttt                                                                                    <  2:611677/251‑1 (MQ=255)
                                                                                                                                                     aCGGGACCACATTGACGAACTTCAGGCTCTCAGCCAGCGTGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATCTGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGTTTCCTGGATCTTAAGGGAGACTTTc                                                                                   >  2:447414/1‑251 (MQ=255)
                                                                                                                                                      ctggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTAGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTaa                                                                                              >  2:627508/3‑239 (MQ=255)
                                                                                                                                                      cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTGCCTGGATCTTaa                                                                                              <  1:627508/239‑1 (MQ=255)
                                                                                                                                                               atggacaaaCTTCAGGATCTCAGCCCGCGTCGTTCGGGAAAGCACCGTCAGCGGAGACTTTTGCTCGCCCTCTGCCCCACCCGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACATTGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCca                                                                         <  2:525016/244‑1 (MQ=255)
                                                                                                                                                                  gACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcagc                                                                       >  2:531436/1‑250 (MQ=255)
                                                                                                                                                                      aaCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGc                                                                  >  1:647505/1‑251 (MQ=255)
                                                                                                                                                                       aCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTATCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCt                                                                 >  2:637197/1‑251 (MQ=255)
                                                                                                                                                                          tCAGGATCTCAGCCATCGTCGTCCGGGAAAGCCCCGTCACCGGAGACTCTGGTCCCCCCTCGGCCACACACTCTTGTTGGGTAAACACCACACCATCGGGCGTCTCGATATGCACAAGGTTGGTATGCCCGCCGACGATGGTTCCACACGACCCTTTCGTACGGTTCCAGACTTTGGCGCGTTAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCtt                                                                <  1:474157/249‑1 (MQ=255)
                                                                                                                                                                          tCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCtt                                                                <  1:696092/249‑1 (MQ=255)
                                                                                                                                                                           cAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTgg                                                             >  2:261334/1‑251 (MQ=255)
                                                                                                                                                                             ggATCTCAGCCAGCGTCGTTCTGGAAAGCACCGCCAGCGGAGACCCTTGCCCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCa                                                           <  1:691507/251‑1 (MQ=255)
                                                                                                                                                                              gATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCAt                                                          >  1:104848/1‑251 (MQ=255)
                                                                                                                                                                                    ggccagcGCCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTGGGTTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGGGACACACGCCCACTTCTCTCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCaa                                                    <  2:468616/250‑1 (MQ=255)
                                                                                                                                                                                         gcgtcttgcTGGAACGCCCCGTCAGTGGAGCCCCTTGCCCGCCCGCTGCCACCCACGCCTGTTGGGTAACCACCCCACCTTCGTCCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGTTGACACTCGCCCCCTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGcc                                              <  2:297930/243‑1 (MQ=255)
                                                                                                                                                                                           gtctttCTGGAAAGCTCCGGCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTg                                            <  1:61265/247‑1 (MQ=255)
                                                                                                                                                                                              gtTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACAAACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTgcgc                                         >  1:645963/1‑251 (MQ=255)
                                                                                                                                                                                                tCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCCCCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCt                                        <  1:80494/250‑1 (MQ=255)
                                                                                                                                                                                                  gggAACGCACCGACAGCGGAGTCTCTTGCTCGCCCTCTGCCACCCGCGCCTGTGGGGTAAGCACCACACCATCGTGCGTCTCGAGATGCACAATGTTGGTATGCCCGCCGACGATGTTGACGCGCGCCCACTTCACACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACGGAGCCTTTCCCCGCCGCCATCAGTGCTTTGGCATCGGCATTGTCCTGCGCTGtt                                     <  1:499479/250‑1 (MQ=255)
                                                                                                                                                                                                      aaGCACCGTCAGCGGAGACTCTTGCTCTCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCCGCCTTTTGCGCGTGAGAAGAGGCTTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCCAGCCGTGCTTTGGCATCGGCAATTGCCTGCGCTTTTGctt                                 >  2:194720/1‑249 (MQ=255)
                                                                                                                                                                                                             gTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                            >  1:446164/1‑249 (MQ=255)
                                                                                                                                                                                                               cAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGTTAACCCCCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCa                        <  1:1780/251‑1 (MQ=255)
                                                                                                                                                                                                                 gCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACt                      <  1:227019/251‑1 (MQ=255)
                                                                                                                                                                                                                     agaCTTTTGCACGCCCTCCGCCACACACGCCTGTTGGTTAAACACCACCCCATCGTGCGTCCCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGGGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTAAGAAGAGGATTGCATCGCAAGGTTCCTGGATCTTAACGTAGACGTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAAGTGCCTGCGCGGTTGCGTCTTTCAGCACTTCCa                  <  2:63490/251‑1 (MQ=255)
                                                                                                                                                                                                                            tGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCgg            <  1:153022/250‑1 (MQ=255)
                                                                                                                                                                                                                            tGCTCGCCCCCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc              <  1:166955/248‑1 (MQ=255)
                                                                                                                                                                                                                             gCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGc           <  1:338911/250‑1 (MQ=255)
                                                                                                                                                                                                                             ccccGCCCTCGGCCACCCACGCCGGTGGGGTAAACACCACCCCACCGTGCGTCTCGATATGCACAAGGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCCTTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCg          <  2:697256/248‑1 (MQ=255)
                                                                                                                                                                                                                              cccGCCCTCCGCCACACACGCCTGTGGGGTAAACACCGCACCATCGTGCGTCTCTAGATGCACAATGTTGGTATTCCCGCCGACGATGGTGACACACGCCCACTTCTCACCTTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGt         <  1:15011/249‑1 (MQ=255)
                                                                                                                                                                                                                                   ccTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGCGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCCGCACTTCCACCCCGGCGTTGGc     >  2:75333/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ctGCCACACACGCTTGTGGGGAAAACCCCACCCCACCTGGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACCCGCCCACTTCTCGCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGGTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAt   <  2:365985/249‑1 (MQ=255)
                                                                                                                                                                                                                                      cggCCACACACGCCTGTTGGGTAAACACCACTCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAtt  <  1:520057/248‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATT  >  NC_000913/3255911‑3256390

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