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

GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATT  >  NC_000913/3255911‑3256390
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
ggATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAAAGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTgg‑g                                                                                                                                                                                                                                         >  2:406547/1‑251 (MQ=255)
ggATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCCCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTGCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTgg                                                                                                                                                                                                                                           <  2:175558/250‑1 (MQ=255)
    gAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑aa                                                                                                                                                                                                                                     <  1:308639/250‑1 (MQ=255)
         aagaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGACGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GG‑AAacacc                                                                                                                                                                                                                                >  2:56569/2‑250 (MQ=255)
          agaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCCCACACGCCTGTTGG‑GT‑AAacacccc                                                                                                                                                                                                                              >  2:237446/1‑249 (MQ=255)
           gaTCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGACTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAacaccac                                                                                                                                                                                                                              >  1:375214/1‑250 (MQ=255)
             tCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACCCACGCCTGTTGG‑GT‑AAacaccacacc                                                                                                                                                                                                                           >  2:696092/1‑251 (MQ=255)
                    ccAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGc                                                                                                                                                                                                                    <  1:32614/251‑1 (MQ=255)
                     cAGCAAACCGCGCTCGCCCTGTTTTTCCAGCGTCGCGCCAATATGCCGCCCCCACTTACCGCTCAACCCTTCCTGCGATAACGCACACTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCAGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCg                                                                                                                                                                                                                   >  2:307396/1‑251 (MQ=255)
                     cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCg                                                                                                                                                                                                                   >  1:468616/1‑251 (MQ=255)
                     cAGCAAACCGCCCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGACCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCGCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCg                                                                                                                                                                                                                   <  2:409589/251‑1 (MQ=255)
                       gCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACGTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGtc                                                                                                                                                                                                                 <  2:164134/251‑1 (MQ=255)
                        cAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGtct                                                                                                                                                                                                                <  1:305556/251‑1 (MQ=255)
                         aaaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGtctc                                                                                                                                                                                                               >  2:78328/1‑251 (MQ=255)
                         aaaCCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGtctc                                                                                                                                                                                                               <  1:570876/251‑1 (MQ=255)
                               cgcTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCCCACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATg                                                                                                                                                                                                         >  1:697109/1‑251 (MQ=255)
                                 cTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGc                                                                                                                                                                                                        >  1:314812/1‑250 (MQ=255)
                                 atCGCAATGTTTATCCATCGTCGCCCCAATATGCAGCCCCCACTTAACGCTCACTCCTTCCTGCGCTAAGGCTCAATTTAGCTTCGCGGAATCGAGCATAAAGCGGATCGCCGCAAACGGGACTTCATTGCCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCCCCGTCAGCGGAGAATCTTGTTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGCGCGTCTCGTTATGc                                                                                                                                                                                                        >  2:499479/2‑250 (MQ=255)
                                           ttttCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACATTCCTCCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGAGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTAAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGt                                                                                                                                                                                             >  1:707602/1‑251 (MQ=255)
                                              tCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                                                                                                                                                                                           <  2:617634/250‑1 (MQ=255)
                                                cAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGcc                                                                                                                                                                                        >  1:525016/1‑251 (MQ=255)
                                                 aGCGTCGCGCCAATATGCAGCCCCCACTTACCGCCCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGccc                                                                                                                                                                                       <  2:602570/251‑1 (MQ=255)
                                                       gcgcCAAAATGCAGCCCCCACTTACCCCTAAAACCTTCCTGCGATAACGCACAATTTAGCTGCGCGGAAGCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTGCGGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCGGCCACCCACGCCTGTTGG‑GT‑AAACCCCACCCCATCGTGCGTCTCGATATGCCCAATGTTGGTATGCCCGCcgac                                                                                                                                                                                 <  2:517887/251‑1 (MQ=255)
                                                           cAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATg                                                                                                                                                                             <  1:237446/251‑1 (MQ=255)
                                                             atatGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCTGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGt                                                                                                                                                                           <  2:691050/251‑1 (MQ=255)
                                                               atGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTg                                                                                                                                                                          >  1:106095/1‑250 (MQ=255)
                                                               atGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTCCATTGACGATCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGa                                                                                                                                                                         <  2:664756/251‑1 (MQ=255)
                                                                 gCAGCCCCCACTTCCCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACGTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACGCTTGCTCGCCCTCTGCCACACCCGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGTGaca                                                                                                                                                                       <  2:469696/251‑1 (MQ=255)
                                                                       cccACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                                                                                                                                                                                           <  2:110706/225‑1 (MQ=255)
                                                                       cccACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTAt                                                                                                                                                                                           >  1:110706/1‑225 (MQ=255)
                                                                         cACTTACCGCTCAACCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCa                                                                                                                                                               <  2:535743/251‑1 (MQ=255)
                                                                                   tCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGCGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCg                                                                                                                                                     >  1:103298/1‑251 (MQ=255)
                                                                                   ccAAACCTCCCTGCAATAACGCGCAAGTTAGCTGCGCGGAATCGAGAGTAAAGCGGATCGCCTCAAACGGGACCGCAGTGACGAACTTCAGGATCGCAGCCAGCGTTGTTGTGGAACGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTCCCACACACGCCTGGTGG‑GT‑AATCTCTACTCCGTCGGTCGTCTCGATATGTACATTGTGGGGATGCCCGCCGTTGATGGTGGCACACGCCCACTTCTCAcc                                                                                                                                                      <  2:70110/249‑1 (MQ=255)
                                                                                     aacccTTCCTGCGATAACGCACAATGTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTGCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTGGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGtt                                                                                                                                                   <  2:478747/248‑1 (MQ=255)
                                                                                     aaaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGtt                                                                                                                                                   >  1:297930/1‑251 (MQ=255)
                                                                                                cGATAACGCACAATTTAGCTGCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACTCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTgg                                                                                                                                        <  1:434204/251‑1 (MQ=255)
                                                                                                    aaCGCACAATATAGCTTCGCGAAAGCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCGCAGCCATCGTCGGTCTGGAAAGCACCTTCAGCGGAGACTCTTGCGCGCCCGCTGCCACACACGCCTG‑TGGTGT‑AAACACCACACCATCGTTGGTCTCGATATTCGCAGTGGTGGTATGCCCGTCGAGGATGTTGACACACGCCCCCGGCTCTCCGTTCCAGACTTTGGcgcg                                                                                                                                    <  2:40845/251‑1 (MQ=255)
                                                                                                      cGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                                  >  1:405635/1‑251 (MQ=255)
                                                                                                         acaATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTg                                                                                                                                  >  1:196338/1‑248 (MQ=255)
                                                                                                            aTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTgaga                                                                                                                               >  2:520057/1‑248 (MQ=255)
                                                                                                            aTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGagaaga                                                                                                                            <  2:205995/251‑1 (MQ=255)
                                                                                                                  cTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGAttt                                                                                                                      >  1:121851/1‑251 (MQ=255)
                                                                                                                       cgtaaTCGAGAATACAGCGGATCGCCGCAAACCGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGCGACTCTTGCTCGCCCTCTGCCACACACCCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGCCACACGCCCACTTCTCACCGTTCCAGCCTTTGACGCGTGCGAAGTGGATTTCATCg                                                                                                                 >  2:506083/4‑251 (MQ=255)
                                                                                                                           aTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAg                                                                                                             <  1:620485/251‑1 (MQ=255)
                                                                                                                                  aTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTg                                                                                                      >  1:247196/1‑251 (MQ=255)
                                                                                                                                      aGCGGCTCGCCGCAAACGGGACTTCATTGCCCACCTTCAGGATCTCAGCCAGCGTCGTTCTGGCACGCACCGTCAGCGGAGACTCTTGCACGACCTCTGCCCCACCCGCCTGTTGG‑GA‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGAAAGGTTCCTGGATc                                                                                                  >  2:430394/1‑251 (MQ=255)
                                                                                                                                       gCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACCCTTGCTCGCCCCCTGCCACACACGCCTGTGGG‑GT‑AAACACCACACCATCGTGCGTCTCGATAGGCACAATGTTGGTATGCCCGTCGACGATGGTGACACACGCCCACATCTCACCGTTCCAGTCGTTGGCGCTTGAGAAGAGGATGTCATCGCAAGGTGCCTTGATCt                                                                                                 <  2:723749/251‑1 (MQ=255)
                                                                                                                                           aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAAc                                                                                             >  2:250329/1‑251 (MQ=255)
                                                                                                                                           aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAAc                                                                                             <  1:373368/251‑1 (MQ=255)
                                                                                                                                           aTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAAc                                                                                             >  1:415693/1‑251 (MQ=255)
                                                                                                                                                   aaaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGAGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACtt                                                                                     >  1:365985/1‑251 (MQ=255)
                                                                                                                                                    aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACttt                                                                                    <  2:611677/251‑1 (MQ=255)
                                                                                                                                                     aCGGGACCACATTGACGAACTTCAGGCTCTCAGCCAGCGTGGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATCTGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGTTTCCTGGATCTTAAGGGAGACTTTc                                                                                   >  2:447414/1‑251 (MQ=255)
                                                                                                                                                      ctggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTAGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTaa                                                                                              >  2:627508/3‑239 (MQ=255)
                                                                                                                                                      cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTGCCTGGATCTTaa                                                                                              <  1:627508/239‑1 (MQ=255)
                                                                                                                                                               atggacaaaCTTCAGGATCTCAGCCCGCGTCGTTCGGGAAAGCACCGTCAGCGGAGACTTTTGCTCGCCCTCTGCCCCACCCGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACATTGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCca                                                                         <  2:525016/244‑1 (MQ=255)
                                                                                                                                                                  gACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcagc                                                                       >  2:531436/1‑250 (MQ=255)
                                                                                                                                                                      aaCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGc                                                                  >  1:647505/1‑251 (MQ=255)
                                                                                                                                                                       aCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTATCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCt                                                                 >  2:637197/1‑251 (MQ=255)
                                                                                                                                                                          tCAGGATCTCAGCCATCGTCGTCCGGGAAAGCCCCGTCACCGGAGACTCTGGTCCCCCCTCGGCCACACACTCTTGTTGG‑GT‑AAACACCACACCATCGGGCGTCTCGATATGCACAAGGTTGGTATGCCCGCCGACGATGGTTCCACACGACCCTTTCGTACGGTTCCAGACTTTGGCGCGTTAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCtt                                                                <  1:474157/249‑1 (MQ=255)
                                                                                                                                                                          tCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCtt                                                                <  1:696092/249‑1 (MQ=255)
                                                                                                                                                                           cAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTgg                                                             >  2:261334/1‑251 (MQ=255)
                                                                                                                                                                             ggATCTCAGCCAGCGTCGTTCTGGAAAGCACCGCCAGCGGAGACCCTTGCCCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCa                                                           <  1:691507/251‑1 (MQ=255)
                                                                                                                                                                              gATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCAt                                                          >  1:104848/1‑251 (MQ=255)
                                                                                                                                                                                    ggccagcGCCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTG‑TGG‑GTTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGGGACACACGCCCACTTCTCTCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCaa                                                    <  2:468616/250‑1 (MQ=255)
                                                                                                                                                                                         gcgtcttgcTGGAACGCCCCGTCAGTGGAGCCCCTTGCCCGCCCGCTGCCACCCACGCCTGTTGG‑GT‑AACCACCCCACCTTCGTCCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGTTGACACTCGCCCCCTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑cc                                              <  2:297930/243‑1 (MQ=255)
                                                                                                                                                                                           gtctttCTGGAAAGCTCCGGCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTg                                            <  1:61265/247‑1 (MQ=255)
                                                                                                                                                                                              gtTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACAAACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTgcgc                                         >  1:645963/1‑251 (MQ=255)
                                                                                                                                                                                                tCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCCCCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCt                                        <  1:80494/250‑1 (MQ=255)
                                                                                                                                                                                                  gggAACGCACCGACAGCGGAGTCTCTTGCTCGCCCTCTGCCACCCGCGCCTGTGGG‑GT‑AAGCACCACACCATCGTGCGTCTCGAGATGCACAATGTTGGTATGCCCGCCGACGATGTTGACGCGCGCCCACTTCACACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACGGAGCCTTTCCCCGCCGCCATCAGTGCTTTGGCATCGGCA‑TTGTCCTGCGCTGtt                                     <  1:499479/250‑1 (MQ=255)
                                                                                                                                                                                                      aaGCACCGTCAGCGGAGACTCTTGCTCTCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCCGCCTTTTGCGCGTGAGAAGAGGCTTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCCAGCCGTGCTTTGGCATCGGCAATTG‑CCTGCGCTTTTGctt                                 >  2:194720/1‑249 (MQ=255)
                                                                                                                                                                                                             gTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTc                            >  1:446164/1‑249 (MQ=255)
                                                                                                                                                                                                               cAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑TT‑AACCCCCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCa                        <  1:1780/251‑1 (MQ=255)
                                                                                                                                                                                                                 gCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACt                      <  1:227019/251‑1 (MQ=255)
                                                                                                                                                                                                                     agaCTTTTGCACGCCCTCCGCCACACACGCCTGTTGG‑TT‑AAACACCACCCCATCGTGCGTCCCGATATGCACAATGTTGGTAGGCCCGCCGACGATGGGGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTAAGAAGAGGATTGCATCGCAAGGTTCCTGGATCTTAACGTAGACGTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAAGTG‑CCTGCGCGGTTGCGTCTTTCAGCACTTCCa                  <  2:63490/251‑1 (MQ=255)
                                                                                                                                                                                                                            tGCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCgg            <  1:153022/250‑1 (MQ=255)
                                                                                                                                                                                                                            tGCTCGCCCCCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc              <  1:166955/248‑1 (MQ=255)
                                                                                                                                                                                                                             gCTCGCCCTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGc           <  1:338911/250‑1 (MQ=255)
                                                                                                                                                                                                                             ccccGCCCTCGGCCACCCACGCCGGTGGG‑GT‑AAACACCACCCCACCGTGCGTCTCGATATGCACAAGGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCCTTTGGCGCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCg          <  2:697256/248‑1 (MQ=255)
                                                                                                                                                                                                                              cccGCCCTCCGCCACACACGCCTGTGGG‑GT‑AAACACCGCACCATCGTGCGTCTCTAGATGCACAATGTTGGTATTCCCGCCGACGATGGTGACACACGCCCACTTCTCACCTTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGt         <  1:15011/249‑1 (MQ=255)
                                                                                                                                                                                                                                   ccTCTGCCACACACGCCTGTTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGCGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCCGCACTTCCACCCCGGCGTTGGc     >  2:75333/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ctGCCACACACGCTTGTGGG‑GA‑AAACCCCACCCCACCTGGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACCCGCCCACTTCTCGCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGGTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAt   <  2:365985/249‑1 (MQ=255)
                                                                                                                                                                                                                                      cggCCACACACGCCTGTTGG‑GT‑AAACACCACTCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAtt  <  1:520057/248‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGG‑GT‑AAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTG‑CCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATT  >  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: ‑