Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913540,5580TG26.1% 194.6 / 19.1 88intergenic (+50/‑7)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base T (12/53);  new base G (23/0);  total (35/53)
Fisher's exact test for biased strand distribution p-value = 9.59e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGC  >  NC_000913/540341‑540792
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
cTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑CTC‑TCCGGTGAAGGGGGGTTTAGTTTCAATTTTGCa                                                                                                                                                                                                                        >  1:373843/1‑239 (MQ=255)
   aCCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGGCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGCC‑ATC‑TCCGGTT‑AGGGGTGTTTTGTTTGAAttttttcagt                                                                                                                                                                                                                     >  2:352632/1‑233 (MQ=255)
        tGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTCTTCCGGTT‑AGGGGTGGTTCTGTTTCATTTTTCCGTCAGcc                                                                                                                                                                                                                 >  2:239287/1‑238 (MQ=255)
           cGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                                                                                                                                                                              <  1:239287/238‑1 (MQ=255)
                                     gCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGGA‑TGC‑TCCCGGTCAGGTGTGGTTATTTTTAAATTTGCGCTCAGCCCCG‑AAAACCCGATATCACGAATTCACTGGt                                                                                                                                                                                  >  2:155411/1‑239 (MQ=255)
                                       ccGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGtt                                                                                                                                                                                 <  1:961543/238‑1 (MQ=255)
                                          tATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCt                                                                                                                                                                                                      <  1:1158532/214‑1 (MQ=255)
                                          tATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGGTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTTAGTGGtttt                                                                                                                                                                               >  1:2113173/1‑237 (MQ=255)
                                                ggCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAg                                                                                                                                                                                                                       <  1:2006751/192‑1 (MQ=255)
                                                ggCCGTACCATTGGCGCGCGTATACCGAAACCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAtt                                                                                                                                                                       <  2:2113173/239‑1 (MQ=255)
                                                    gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTCATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCAGTTT‑TCCGGTT‑TGGTGTGGGTTTGTTTAATTTTTTAGTCGGGCGCGCAAAGCC‑GTTAttcgga                                                                                                                                                                                            >  2:1158532/1‑209 (MQ=255)
                                                     tACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGGTTAGGGGTGGTTATGTTCAATTTTGCAGTCAGGCGCG‑AAAGCCTGGTATCAGGATTTTCGTGGtttttttttttttttt                                                                                                                                                                   >  1:87545/1‑238 (MQ=255)
                                                     tACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCCGTTAAGGGGGGTTTATGTTCAATTTTGCAGTTAGCCCCG‑AAAAACTGTTATCAGGATTTCCGGGGTTTTTCTTTTTTTTTTg                                                                                                                                                                  >  1:1192130/1‑239 (MQ=255)
                                                      aCCATTGGCGCGCGTATCCCGAAAACCATCTTACGTGGTGATGTGATTTACGACAGTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGc                                                                                                                                                                 <  1:1244139/239‑1 (MQ=255)
                                                        cATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGGTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAtt                                                                                                                                                                       <  1:564526/231‑1 (MQ=255)
                                                          ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑CTC‑TCCCGTGAAGGTGTGTTTATGTTTAATTTTTCAGTCAGCcgcgc                                                                                                                                                                                                             >  2:2006751/1‑191 (MQ=255)
                                                           tGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac                                                                                                                                                            <  1:1614958/239‑1 (MQ=255)
                                                           tGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac                                                                                                                                                            <  1:1080911/239‑1 (MQ=255)
                                                               gcgcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACTGtt                                                                                                                                                        <  1:388727/239‑1 (MQ=255)
                                                               gcgcgTATCACGAAAACCATCGTACGTGGTGATGTTATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCATCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTGCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                                        <  1:433849/239‑1 (MQ=255)
                                                                  cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTT‑TCCGGGTATGGGGTGGGTATGTTTAATTTTTTAGTGTGCCGCG‑GAAGCCTGTTATCAGGAGTTCTGTGGGTTTTTTTTTTTTTTTGTAacacc                                                                                                                                                           >  2:564526/1‑232 (MQ=255)
                                                                     acaaCGAACACCATCTTACGTGGTTATGTGATTTACGACATTGAACATGGCTTCCCTGTTGCGCCGAAAGGTCAAGTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTTCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                                                                                                                   <  2:1074984/235‑1 (MQ=255)
                                                                     aTCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                                                                                                                   <  2:1881923/238‑1 (MQ=255)
                                                                     aTCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                                                                                                                   <  1:137700/238‑1 (MQ=255)
                                                                                      aCGGGGTGATGTGATTTTCGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTc                                                                                                                                 <  1:122167/239‑1 (MQ=255)
                                                                                        gtggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTg                                                                                                                               <  1:1216775/239‑1 (MQ=255)
                                                                                          ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGGGTGGTTATGTTCAATTTTGCAGTCAGCCGCG‑CAAGCCTGTTATTAGGGTTTTAGGGGGTTTTCTTTTTTTTTTTCAACCCGGTTTTGGTTTCTCCTCAGCCACTTTTTGCt                                                                                                                             >  2:1216775/1‑239 (MQ=255)
                                                                                            tgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                                                                                                                                                                              <  1:535729/157‑1 (MQ=255)
                                                                                            tgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTAt                                                                                                                                                                                                 <  1:1810181/169‑1 (MQ=255)
                                                                                                gtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCC‑TTC‑TCCTCGTTAGGGGTGTTTTTTTTCAATTTTGCCGTCAGCCCCC‑AAAGAGCGTTTTTCGGGTTTTTGTGGTTTTTTTTTTTTTTTTTCTACACGGGTGTGGGTTCTCCCACCCTCATTTTTTCtttttta                                                                                                                       >  2:168824/1‑236 (MQ=255)
                                                                                                   atTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑GAAGc                                                                                                                                                                                                        >  1:1898461/1‑155 (MQ=255)
                                                                                                   atTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGc                                                                                                                                                                                                        <  2:1898461/155‑1 (MQ=255)
                                                                                                      tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TTCGGTTGAGGTGTGTTTATGTTCAATTTTCAGGTCAGCCGCC‑GAAGACTTGTTTTTGGGTTTTAGTAGGTTTTTTTTATTTTTTGCAACACGGTTTTCGTTTCTCCTTCGCTGCTTTTTGCGTTTTAGGTGca                                                                                                                  >  2:164854/1‑237 (MQ=255)
                                                                                                      tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGGTACGGGGGGGTTATGTTTTATTTTTCAGGCAGCCACG‑AAAGCCTggt                                                                                                                                                                                                   >  2:535729/1‑155 (MQ=255)
                                                                                                      tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCCGGTAAGGGGTGTTTTTGTTCAATTTTTCAGTCAGCCGCG‑CAAGCCTGTTATCAggtatttt                                                                                                                                                                                       >  2:1810181/1‑163 (MQ=255)
                                                                                                               gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCC‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                        <  2:1196293/239‑1 (MQ=255)
                                                                                                               gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCa                                                                                                               <  1:433330/232‑1 (MQ=255)
                                                                                                               gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                        <  1:1542610/239‑1 (MQ=255)
                                                                                                               gAAAAGGGCTTCCCTGGTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                        <  1:1155712/239‑1 (MQ=255)
                                                                                                                aaCAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTtcc                                                                                                                                                <  1:1001146/198‑1 (MQ=255)
                                                                                                                       cTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTa                                                                                                 <  1:669490/238‑1 (MQ=255)
                                                                                                                          cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                        <  1:993648/228‑1 (MQ=255)
                                                                                                                          cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGGTTATGTTCAATTTTTCAGTCAGCCGCG‑AAAGCCTGGTATCCGGATTTCAGTGGTTTTTTTTTATTTTTTGCAACACGGTTGGGGTTCCTCCCACGCTa                                                                                                                                      >  2:1001146/1‑198 (MQ=255)
                                                                                                                             tGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATc                                                                                                                                                                                                <  1:461578/137‑1 (MQ=255)
                                                                                                                             tGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATc                                                                                                                                                                                                >  2:461578/1‑137 (MQ=255)
                                                                                                                              gTTCCGCCGAAAGTCCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACg                                                                                         <  2:373843/239‑1 (MQ=255)
                                                                                                                                    ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTTAGGGGTGTTTATTTTTAATTTTTCTGTCCGCCGCG‑CAAGCCTGTTTTCCGGGTTTTTGTTGGTTTTTTTTTTTTTTTTCTACACTGTTTTGGTTTCTCCTCCCCTTCTTTTTTCTTTTTTGTTTTCCCAAAATATACTTCTTTCGCTccttctcttt                                                                                   >  2:993648/1‑230 (MQ=255)
                                                                                                                                       aaaGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTc                                                                                                                                                                                        >  1:560652/1‑135 (MQ=255)
                                                                                                                                       aaaGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTc                                                                                                                                                                                        <  2:560652/135‑1 (MQ=255)
                                                                                                                                              aaTTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca                                                                                                                                                             <  1:789673/155‑1 (MQ=255)
                                                                                                                                                      cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTTCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCtttt                                                                             <  1:310659/227‑1 (MQ=255)
                                                                                                                                                      cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                 <  1:155411/239‑1 (MQ=255)
                                                                                                                                                      cTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                 <  1:1930245/239‑1 (MQ=255)
                                                                                                                                                        tAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCCGTTTAGGGGTGTTTTTTTTCAATTTTTCAGTCAGCCGCG‑GAAGCCTGTTTTCAGGGTTTTTGTGGTTTTTTTTTTTTTTTTTCAACCAGGTTGTGGt                                                                                                                                                   >  2:789673/1‑155 (MQ=255)
                                                                                                                                                           aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGc                                                            <  1:426216/239‑1 (MQ=255)
                                                                                                                                                              tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt                                                                                                                                                        <  1:841252/144‑1 (MQ=255)
                                                                                                                                                               cagcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTT‑TCCCGTTTAGGGGGGTTTTTGTTCTATTTTTCAGTCTGCCCCG‑GCAACCCGTTTTCCGGGTTTTCGTGGTTTTTTTTTTTTTTTTTCAACACGGTTGTGGTTCCTCCTACCCCACTTTTTTCTTTTTCTTTTCCGCCAAGTATCCTGCTTTCGCTCACCCAATATGCTTTTCTTGTTTCCCCACCGGCCt                                                          >  2:649056/1‑237 (MQ=255)
                                                                                                                                                                agcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑ATC‑TCCCGTTAAGGGGGGTTTATGTTCCATTTTGCAGGCCGCCGCG‑AAAGCCTGTTTTCAGGGTTTCTGTGGTTTTTTTTTATTTTTTTCAACACGGGTTGGGTTTCTCCTACGCTCCTTTCTGTTTTTTAGTTTCCGCACAGGAGCCTGTCTACGCTCACCCAATATGCTTTTCTTTCTCCCGcccagccctcg                                                        >  2:1494723/1‑228 (MQ=255)
                                                                                                                                                                  cagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTAt                                                                                                                                                                                                 <  2:682266/99‑1 (MQ=255)
                                                                                                                                                                     tactCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCt                                                                                             <  1:291931/193‑1 (MQ=255)
                                                                                                                                                                      aaTCTGGCCCCTGCAGTGCCCGGCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCCCGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCt                                                  <  2:1192130/238‑1 (MQ=255)
                                                                                                                                                                       aTCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGcc                                                                                                                  <  1:473803/173‑1 (MQ=255)
                                                                                                                                                                       aTCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAACCTGTTATCAGGATTTCAGTGGTTTTTCTTTAtt                                                                                                                                                                       <  2:1072967/120‑1 (MQ=255)
                                                                                                                                                                        tCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGGGGGGTTTTGTTCAATTTTGGAGTCCGCCGCG‑AAAGCCTGTTTTTAGGATTTTAGTGGGTTTTCTTTATTTTTTTCAACACGGTTTTGGGttctc                                                                                                                                              >  2:841252/1‑144 (MQ=255)
                                                                                                                                                                            gCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCT‑TTT‑CCCGGGTAAGGGGGGTTTATTTTTAATTTTGCAGGCACCCCCG‑AAAGCCCGTTTTCACGATTTCa                                                                                                                                                                                       >  1:682266/1‑99 (MQ=255)
                                                                                                                                                                            gCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑a                                                                                                                                                                                                            >  1:1483528/1‑78 (MQ=255)
                                                                                                                                                                            gCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑a                                                                                                                                                                                                            <  2:1483528/78‑1 (MQ=255)
                                                                                                                                                                              cccTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                        <  1:1784842/176‑1 (MQ=255)
                                                                                                                                                                               ccTGCAATGCCCGTCCTTGCGGC‑GGGCC‑TTC‑TCCGGGTAAGGGGTGTTTATTTTCAATTTTGCAGTCTGCCGCG‑GAAGCCTGTTATCTGGGGTTCAGTGGTTTTTCTTTATTTTTTTCAACACCGGTGGGGGTCCTCCTCCGCTCCTTTCTTCTTTTTTGTTTCCGCCAAGTAACCTGCTTTCTCTCACGCCatat                                                                                   >  2:291931/1‑196 (MQ=255)
                                                                                                                                                                                 tGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCt                                                                           <  1:357228/202‑1 (MQ=255)
                                                                                                                                                                                 tGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGCTTAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTTAGTGGTTTTTTTTTTTTTTTTGCAACACGGTTGTGGGTCCTCCTACGCTACTTTCTTTTTTTTAGTTGCCGCCCAGTAGc                                                                                                        >  2:473803/1‑173 (MQ=255)
                                                                                                                                                                                 tGCAATGCCCGTCCTTGCGGC‑GGGCA‑ATC‑TCCGGGTAAGGTGTGTTTTTGTTCAATTTTTCATTCAGCCGCG‑AAAGCCTGTTATTTGGATTTCTGAGGttttttttttttttttttcaaa                                                                                                                                                             >  1:1072967/1‑114 (MQ=255)
                                                                                                                                                                                   cAATGCCCGTCCTTGCGGC‑GGGCC‑TTG‑TCCGGGTTATGGGTGGTTTTGTTCTATTTTTCTGTCCACCGCC‑ACAGCCCTTTTTCAGGATTTCAGTGGGTTTTCTTTATTTTTTTCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTTCTTTTCCGTTTCCGCAAAATAGCCTGCTTACGCTCACGCAATATGCTTTTTTTTCCTCCTCCCTGGCCTGCTTCTCTCAAGCGgttt                                        >  2:31503/1‑235 (MQ=255)
                                                                                                                                                                                           gTCCTTGCGGC‑GGGCA‑TTC‑TCCGGGTAAGGGGTGTTTATGTTCAATTTTGCAGTCCGCCGCG‑AAAGACTGTTTTCAGGGTTTTAGTGGGTTTTTTTTTTTTTTTTCAACACGGGTGTGGTTCCTCCTACCCTACTTTTTGCTTTTTAGTTTCCGCAAAGTAGCCTGCTTTCGCTCCCGCAATATGCTTTTCTTTCTACCCCa                                                                 >  2:357228/1‑202 (MQ=255)
                                                                                                                                                                                            tCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATGTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                           <  1:624676/239‑1 (MQ=255)
                                                                                                                                                                                             cctgtgcgc‑ggGAA‑TTC‑TCCGGTTAAGGTGTGTTTATTTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGTTTGTGGTTCCTCCTACGCTCCTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                           <  2:938375/231‑1 (MQ=255)
                                                                                                                                                                                             ccTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                           <  1:168824/238‑1 (MQ=255)
                                                                                                                                                                                              cTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                         <  1:817608/239‑1 (MQ=255)
                                                                                                                                                                                              cTTGCGGC‑GGGCA‑TTC‑GCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                         <  1:587560/239‑1 (MQ=255)
                                                                                                                                                                                                        ggCA‑TTC‑TCCGGTTTAGGGGTGTTTATGTTCAATTTTGCAGTCCGCCCCG‑GAAGCCTGTTATCAGGATTTTTGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCCCTTGCCGCCAAGTAGCCTGCTTACGCTCCCGCACTATGCTTTTCTTGCTACCGCCCTCGCCTGCCTCGCTCACGCGGTTTGTGGTCCtcg                              >  2:2003509/1‑225 (MQ=255)
                                                                                                                                                                                                           a‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGaa              <  1:1827261/238‑1 (MQ=255)
                                                                                                                                                                                                                    ggTTAAGGTGTTTTTATGTTCAATTTTGCAGTCAGCCNNG‑aa                                                                                                                                                                                                           <  2:1510885/42‑1 (MQ=255)
                                                                                                                                                                                                                        aaGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGTTTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATg                                                                                  >  2:1285166/1‑159 (MQ=255)
                                                                                                                                                                                                                        aaGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATg                                                                                  <  1:1285166/159‑1 (MQ=255)
                                                                                                                                                                                                                        aaGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGc  <  1:352632/239‑1 (MQ=255)
                                                                                                                                                                                                                         aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                           >  2:37313/1‑213 (MQ=255)
                                                                                                                                                                                                                         aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                           <  1:37313/213‑1 (MQ=255)
                                                                                                                                                                                                                         aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTgg   >  1:1201722/1‑237 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
CTTACCAATGACGATCTGGAATATCGCCACAAAGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGC‑GGGCA‑TTC‑TCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCG‑AAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGC  >  NC_000913/540341‑540792

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑