Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,550,108 T→C 5.6% Q789R (CAG→CGG)  malP ← maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,550,1080TC 5.6% 240.9 / 4.3 71Q789R (CAG→CGG) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (28/39);  new base C (2/2);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

AGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACC  >  NC_000913/3549839‑3550392
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            
agcTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCg                                                                                                                                                                                                                                                                           >  2:169226/1‑289 (MQ=255)
   tatcTTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGAt                                                                                                                                                                                                                                                                         <  1:580211/288‑1 (MQ=255)
            cgcTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGACGTTGGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATAGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAACCGCGAATAGGGAGATACGAGCCaaa                                                                                                                                                                                                                                                               >  2:294691/1‑289 (MQ=255)
              cTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACa                                                                                                                                                                                                                                                             <  1:294691/289‑1 (MQ=255)
                          accaTCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAgcggg                                                                                                                                                                                                                                                 <  1:207123/289‑1 (MQ=255)
                           ccaTCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAgcgggc                                                                                                                                                                                                                                                >  2:554905/1‑289 (MQ=255)
                                              gCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCcgcg                                                                                                                                                                                                                             <  1:324413/289‑1 (MQ=255)
                                                ttgTATAAACCAGGACTTGCGGGGCTGGGGTTCCCGCCACTTGCGTGTCGGTACGTTGATGCGTCGCGTAAAGAAAACGCCGTTTGGTTTCGGTGCTAATCGGCTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGTTCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCcgcgc                                                                                                                                                                                                                            <  2:351740/286‑1 (MQ=255)
                                                tGGTATAAACCATGACATTCGTGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCcgcgc                                                                                                                                                                                                                            <  1:584101/288‑1 (MQ=255)
                                                tGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCcgcgc                                                                                                                                                                                                                            <  1:212637/288‑1 (MQ=255)
                                                            tGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGcct                                                                                                                                                                                                                >  1:578201/1‑288 (MQ=255)
                                                            tGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGcct                                                                                                                                                                                                                >  1:354544/1‑288 (MQ=255)
                                                             gACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGcctc                                                                                                                                                                                                               >  2:373490/1‑288 (MQ=255)
                                                               cATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGcctcct                                                                                                                                                                                                             <  1:509522/288‑1 (MQ=255)
                                                                  tCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTTGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAACGCGAACCGAACTAAACATACCCCAGCCGGCGGTATTGAGGAACGCCGCGCTAGTGCCGGCCTGCTggg                                                                                                                                                                                                          >  2:126/1‑287 (MQ=255)
                                                                   cGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTcg                                                                                                                                                                                                        <  1:169226/289‑1 (MQ=255)
                                                                       aCTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCgg                                                                                                                                                                                                     <  2:578201/288‑1 (MQ=255)
                                                                         tGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTAc                                                                                                                                                                                                  >  1:163964/1‑289 (MQ=255)
                                                                           gCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCa                                                                                                                                                                                                                                                      >  2:192131/1‑235 (MQ=255)
                                                                           gCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCa                                                                                                                                                                                                                                                      <  1:192131/235‑1 (MQ=255)
                                                                               tACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGcac                                                                                                                                                                                             >  1:349562/1‑288 (MQ=255)
                                                                               tACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGcac                                                                                                                                                                                             <  1:568047/288‑1 (MQ=255)
                                                                                     cACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCCAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCCCGAATAGAGCGATCCGAGCTAAACATACCCCAGCGGGCGGCATTGAGGATCGCCCCCCGGATCCAGGCCTCCTGGTCGCGGTaaacgcctccccc                                                                                                                                                                                       >  2:320685/1‑276 (MQ=255)
                                                                                      aCTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGCCCAGGCCTCCTGGTCGCGGTACAGCACATCCa                                                                                                                                                                                        >  1:100025/1‑286 (MQ=255)
                                                                                       cTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCt                                                                                                                                                                                     >  1:426335/1‑288 (MQ=255)
                                                                                        tttCGGGGCGGTACGTTGATGCTGCGCGTCAAGCAAAGGCGGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGc                                                                                                                                                                                   <  2:387853/289‑1 (MQ=255)
                                                                                          tCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACACCCACCTGCtt                                                                                                                                                                                 >  1:519439/1‑289 (MQ=255)
                                                                                                 ggTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGcc                                                                                                                                                                          <  1:373490/289‑1 (MQ=255)
                                                                                                 ggTACGTGGATGCAGCGCGTGAAGGAAACGCCGTTTGGTTTGTGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGc                                                                                                                                                                           <  2:519439/288‑1 (MQ=255)
                                                                                                  gTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGAGcc                                                                                                                                                                          <  2:349562/288‑1 (MQ=255)
                                                                                                   tACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCct                                                                                                                                                                         <  1:554905/288‑1 (MQ=255)
                                                                                                                          aaaCGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGc                                                                                                                                                   <  2:466687/287‑1 (MQ=255)
                                                                                                                           aaCGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCAt                                                                                                                                                <  1:474518/289‑1 (MQ=255)
                                                                                                                             cGCCGTTGGGTTTCGCGCCTAATCGACTGCGGTTTACCGTGGGCCTTGATGTAATTGGGCATAATCGCGGCCGCAGGCGGGGCAGTAGCCAGACGTTTGGTATCCATCGAGCGTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCa                                                                                                                                              <  2:455294/289‑1 (MQ=255)
                                                                                                                              tcggTTTGGTTTCGTCGCTAATCGACTGCGGTTTACCGTCGGCATTGATGTAATTGGAGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTAGTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCAc                                                                                                                                             <  1:560453/286‑1 (MQ=255)
                                                                                                                                                          cGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCg                                                                                                                 >  1:232206/1‑289 (MQ=255)
                                                                                                                                                            gTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGcc                                                                                                                                                                          <  1:581560/230‑1 (MQ=255)
                                                                                                                                                            gTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGcc                                                                                                                                                                          >  2:581560/1‑230 (MQ=255)
                                                                                                                                                             tttACAGTGGGCATTGGTGTAATTGGGGGTAATCCCCGCCGCCTGCTCGGCATTATCCAGACGTTTGCTTTCCTGCGAGCTTCCTTAGCGTTTTGCGTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGGGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGGACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCGCTGTTTGCCGATa                                                                                                              <  2:298407/289‑1 (MQ=255)
                                                                                                                                                                aCCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTcg                                                                                                                                                                                                        >  1:234940/1‑196 (MQ=255)
                                                                                                                                                                aCCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTcg                                                                                                                                                                                                        <  2:234940/196‑1 (MQ=255)
                                                                                                                                                                  cGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACtgtg                                                                                                         <  2:163964/289‑1 (MQ=255)
                                                                                                                                                                     gggCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTg                                                                                                                        >  1:159767/1‑271 (MQ=255)
                                                                                                                                                                     gggCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTg                                                                                                                        <  2:159767/271‑1 (MQ=255)
                                                                                                                                                                      ggCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCa                                                                                                                                                 >  2:84888/1‑245 (MQ=255)
                                                                                                                                                                      ggCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCa                                                                                                                                                 <  1:84888/245‑1 (MQ=255)
                                                                                                                                                                                  aTTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGc                                                                                                                   >  2:231148/1‑263 (MQ=255)
                                                                                                                                                                                  aTTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGc                                                                                                                   <  1:231148/263‑1 (MQ=255)
                                                                                                                                                                                        gCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCAc                                                                                                                                             <  2:169369/231‑1 (MQ=255)
                                                                                                                                                                                        gCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCAc                                                                                                                                             >  1:169369/1‑231 (MQ=255)
                                                                                                                                                                                        gCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCAt                                                                                     >  1:140194/1‑287 (MQ=255)
                                                                                                                                                                                                 cgccgcCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCg                                                                          >  2:519953/1‑289 (MQ=255)
                                                                                                                                                                                                     gccAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGc                                                                      <  2:426335/289‑1 (MQ=255)
                                                                                                                                                                                                               cATTAGCCGGACGTTTGCTTTCCGTCGGTCTTCCTGACGCTTTTTGCCGCCTGATACGGGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTAc                                                             <  2:80935/288‑1 (MQ=255)
                                                                                                                                                                                                                          cGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGGGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGc                                                  <  2:140194/288‑1 (MQ=255)
                                                                                                                                                                                                                              tGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCttt                                             >  1:141266/1‑289 (MQ=255)
                                                                                                                                                                                                                              tGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCCGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTcg                                                                                                                                                                                                        >  2:208656/1‑134 (MQ=255)
                                                                                                                                                                                                                              tGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCCGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTcg                                                                                                                                                                                                        <  1:208656/134‑1 (MQ=255)
                                                                                                                                                                                                                                 tttCCATCGAGCGTCCTTAGCGTTTGGCCTGCCAGATACGGGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCaa                                          <  2:232206/289‑1 (MQ=255)
                                                                                                                                                                                                                                         gAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGcc                                                                                                                                                                          >  1:263970/1‑153 (MQ=255)
                                                                                                                                                                                                                                         gAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGcc                                                                                                                                                                          <  2:263970/153‑1 (MQ=255)
                                                                                                                                                                                                                                                          ttGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCTCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCCCTGTATTTACCGCTTTCCAGCTCTTTCAATCCTGCGTCAAGCCCCTTAtcttt                  >  2:298675/1‑289 (MQ=255)
                                                                                                                                                                                                                                                          ttGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGc                                                                                       <  2:376290/220‑1 (MQ=255)
                                                                                                                                                                                                                                                          ttGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGc                                                                                       >  1:376290/1‑220 (MQ=255)
                                                                                                                                                                                                                                                          ttGCCTGCCAGATACGAGCCCGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTAtcttt                  >  2:44176/1‑289 (MQ=255)
                                                                                                                                                                                                                                                           tgcctgccAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTAtctt                   >  2:20297/1‑287 (MQ=255)
                                                                                                                                                                                                                                                              ctgccAGATACGAGCCCGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTttcttc              <  1:44176/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                gccAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCg            <  1:320685/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCCTGGTAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCAcc  <  2:141266/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCAcc  <  1:298675/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCAcc  <  1:519953/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            
AGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGTTACCGCCACTTTCGTGGCGGTACGTTGATGCATCGCGTCAAGCAAACGCCGTTTGGTTTCGGCGCTAATCGACTGCGGTTTACCGTGGGCATTGATGTAATTGGGGCTAATCCCCGCCGCCAGCGCGGCATTATCCAGACGTTTGCTTTCCATCGAGCTTCCTTAGCGTTTTGCCTGCCAGATACGAGCCTGATAATCGCGAATAGAGCGATCCGAGCTAAACATACCGCAGCGGGCGGTATTGAGGATCGCCGCGCGAGTCCAGGCCTCCTGGTCGCGGTACAGCACATCCACCTGCTTTTGTGCCTCTACATAGGCTGCGAAATCCGCCATCACCAGATACGGATCGCCGCCCTGTTTGCCGATACTGTGCAGCATCTGGTCGAAGGCATGCTTATCGCCGTCGCTGTATTTACCGCTTTCCAGCTCTTTCAATACTGCGTCCAGCACCTTATCTTTCTTCCGCCATTTCACC  >  NC_000913/3549839‑3550392

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