Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,185,1990TG28.2% 172.0 / 10.6 78T35P (ACC→CCC) thiMhydroxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (9/47);  new base G (22/0);  total (31/47)
Fisher's exact test for biased strand distribution p-value = 1.42e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAA  >  NC_000913/2184955‑2185444
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
aGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGG‑GTTGcac                                                                                                                                                                                                                                                    >  1:229395/1‑251 (MQ=255)
 gTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGcacc                                                                                                                                                                                                                                                   <  2:69835/251‑1 (MQ=255)
       gACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGCTCAAGCGTCCAGGGTGGTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑tgtggcacacccc                                                                                                                                                                                                                                              >  2:30902/1‑238 (MQ=255)
          aaaaaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTTAAGG‑GTTGGCCCAcatacttt                                                                                                                                                                                                                                          >  2:264770/1‑246 (MQ=255)
           aaaaTGGCGGCGATAATTGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTGcccccctcatttg                                                                                                                                                                                                                                          >  1:92468/1‑239 (MQ=255)
            aaaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGt                                                                                                                                                                                                                                        <  2:223795/251‑1 (MQ=255)
             aaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTTGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGGTTCGATAACCATCGCTGGCGATGCAACGAGCGCCAGCAAGGTATTGGCGGGAAAGG‑GTTGtaccaccccatttggt                                                                                                                                                                                                                                       >  1:315630/1‑241 (MQ=255)
                ggcggcGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGG‑GTTGTACCACATCAtttgtattg                                                                                                                                                                                                                                    >  2:235334/1‑244 (MQ=255)
                     cgATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTACTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACAATCGCTGGCGATGCACCGAGCGCCAGCAAGGCATTGGCGGTAAAAG‑TTTGCACCCCCACATTGGTcctatcgtg                                                                                                                                                                                                                               >  1:141979/1‑243 (MQ=255)
                       aTAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCATTGGTCATGCAGTGcc                                                                                                                                                                                                                             >  1:189719/1‑250 (MQ=255)
                          aTCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCATTTGTCCTGCAGGGcaca                                                                                                                                                                                                                           >  1:364033/1‑250 (MQ=255)
                            cGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAag                                                                                                                                                                                                                         <  1:52188/250‑1 (MQ=255)
                                cgcACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAgggg                                                                                                                                                                                                                    <  1:411079/251‑1 (MQ=255)
                                   ccccgccGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAt                                                                                                                                                                                                                 <  2:117869/246‑1 (MQ=255)
                                           cTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCCTTTGGCATGCAGTGCACAAGAGGGGGATGTTgggg                                                                                                                                                                                                          >  1:340083/1‑248 (MQ=255)
                                                ggATCAAGGGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCGGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑aaaaa                                                                                                                                                                                                    <  2:250446/250‑1 (MQ=255)
                                                       gCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAAc                                                                                                                                                                                            <  2:217519/251‑1 (MQ=255)
                                                         gTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCAGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTTAAGG‑GTTGCACCACATCATTGGTCTAGAAGAGCAAAAGAGGGGGAGGTTGGTG‑AAAAAAGGGTAAcgc                                                                                                                                                                                          >  2:109890/1‑251 (MQ=255)
                                                           ccAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGt                                                                                                                                                                                        <  1:44747/251‑1 (MQ=255)
                                                               ggTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTgcgc                                                                                                                                                                                    <  1:12802/251‑1 (MQ=255)
                                                                  gtTTGAGAGCTTTTTGCTTGCTCAACGGCAGCCCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGa                                                                                                                                                                                 <  2:257006/251‑1 (MQ=255)
                                                                      gagaGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTgc                                                                                                                                                                             <  1:406447/251‑1 (MQ=255)
                                                                        gagCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAG‑GTTGCACCACATCATTGGGCATGCAGTGCACAAGAGGGGGAAGTTGGTGGAAAAAGG‑GTAAAGCGTGCGCAGATTgtgc                                                                                                                                                                           >  2:121134/1‑248 (MQ=255)
                                                                                    tgctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagca                                                                                                                                                                <  1:315872/250‑1 (MQ=255)
                                                                                    tgctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagc                                                                                                                                                                 <  2:40942/249‑1 (MQ=255)
                                                                                     gctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTG‑ACCACATCATTGGGCATGCCGGGCACAAGAGGGGGATGTTGGGG‑AAAAAGGTGTAAAGCGTGCGCAGATTGCGCCGAACCCCGCAgg                                                                                                                                                              >  1:241098/1‑251 (MQ=255)
                                                                                     gctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGGG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagcag                                                                                                                                                               <  2:258968/250‑1 (MQ=255)
                                                                                        cAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                                                                                                                            <  2:315630/250‑1 (MQ=255)
                                                                                        cAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCg                                                                                                                                                           <  2:241508/251‑1 (MQ=255)
                                                                                             gcagcaCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTg                                                                                                                                                      <  1:391426/251‑1 (MQ=255)
                                                                                               agcaCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTG‑GAAAAGGTGTAACGCGTGCGCAGATTGCGCTGGACCCAGCAGGGCGACTTGcc                                                                                                                                                    >  1:392117/1‑250 (MQ=255)
                                                                                                gcaCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAt                                                                                                                                                   <  1:135652/251‑1 (MQ=255)
                                                                                                    gCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCAAAACATTGGGCATGCAGTGCACAAGAGGGGGATGGTGGGG‑AAAAAGGGGGAACGCGGGCGCAGATTGGGCTGAACCCAGCCGGCCGACTTggaaaggt                                                                                                                                               >  1:32793/1‑244 (MQ=255)
                                                                                                     cATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGttt                                                                                                                                             <  1:306604/251‑1 (MQ=255)
                                                                                                         gCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTc                                                                                                                                         <  2:340083/251‑1 (MQ=255)
                                                                                                         gCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCATTGGTCATTCAGTGCACAAGAGGGGGATGTTGGGG‑AAAAAGGGGTAAAGCGTGCGCCGATTGCGCCGAACCCAGCAGGGCGACTTGCAT‑AGTTTtctc                                                                                                                                         >  1:97662/1‑251 (MQ=255)
                                                                                                             gAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGc                                                                                                                                     <  1:242943/251‑1 (MQ=255)
                                                                                                               gcgcGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCCCATCAATGGTCATGAAATGCACAAGAGGGGGATGTTGGGG‑GAAAAAGTGTAACGCGGGGGCAGATTGCGCTGAACCCCGCAGGTCGGCTTGCCT‑AGTTTTCTCCCGccc                                                                                                                                   >  1:400970/1‑250 (MQ=255)
                                                                                                                  cgTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATaa                                                                                                                                <  1:98348/251‑1 (MQ=255)
                                                                                                                       cTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGGATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAAACCAGCAGGGCGGCTTGCAT‑AGTTTGCTCCTGCCATAACgcg                                                                                                                            >  1:247315/1‑248 (MQ=255)
                                                                                                                         gCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTgaa                                                                                                                          <  1:418668/250‑1 (MQ=255)
                                                                                                                             cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TATGCACCACATCATTGGTCATGCAGTGCACAAAAGGGGGATGGTGGTG‑AAAAAGGTGTTCAGCGTGTGCCGAGTGCGCTGGACCCAGCAGGGCGACTTG‑ATGAGTTTGCTCCTCCCCTAACGTgaaaaaac                                                                                                                     >  2:8681/1‑246 (MQ=255)
                                                                                                                               gtgtGCCAACGTTAATCAACAAGGCCCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCaa                                                                                                                   <  2:91597/251‑1 (MQ=255)
                                                                                                                                 gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACCTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATg                                                                                                                 <  1:43946/251‑1 (MQ=255)
                                                                                                                                    ccAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTTCACCACATCATTTGGCATGCAGTGCACAAGAGGGGGATGTTGGGG‑AAAAAAGTGTAACGCGTGCGCCGATTGCGCTGAACCCAGCAGGGCGACTTGCAT‑AGTTTGCTCCCGCCCTAACCTGAAAAAGCAATGAcc                                                                                                              >  1:303/1‑251 (MQ=255)
                                                                                                                                       aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtgg                                                                                                           <  1:50870/251‑1 (MQ=255)
                                                                                                                                       aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtgg                                                                                                           <  2:424674/251‑1 (MQ=255)
                                                                                                                                            aaTCAACAAGGCACTGGCGATAGCCGCAAACTGCCTGGCCTCTGCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTc                                                                                                      <  2:187730/251‑1 (MQ=255)
                                                                                                                                            aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTggtggt                                                                                                       <  1:35243/250‑1 (MQ=255)
                                                                                                                                            aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTggtggt                                                                                                       <  1:180342/250‑1 (MQ=255)
                                                                                                                                              tCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACCTCAATGGGTATTCAGGGGACAAAAGGGGGATGTTGGGGGAAAAAGG‑GCAACGCGTGCGCAGATTGCTCGGGAGCCCCCAGGGCGGCTTGCAT‑AGTTTTCTCCTTCCATAACGTGacgaaggaactaccgggggggacg                                                                                                    >  2:417869/1‑233 (MQ=255)
                                                                                                                                               caacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCg                                                                                                    <  1:30902/250‑1 (MQ=255)
                                                                                                                                                aacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCACTGGGCATGCAGTGCACAAGGGGGGGATGTTGGGG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGGTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCAAAAAGGGAAGAAGGAATGAACCGGTGGGCCggg                                                                                                  >  2:133532/1‑249 (MQ=255)
                                                                                                                                                 acaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTGCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGa                                                                                                <  2:99856/251‑1 (MQ=255)
                                                                                                                                                   aaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCACATCATTGGGCCTGCAGTGCAAAAGAGGGGGATGTTGGTG‑TAAAAAGAGTAACGCGTGGGCAGCTGGCCCTGAACCCAGCAGGGGGACTTGCCT‑AATTTGCTCCCGCCATAACGAGAAGAAGCAATGACCTGGTGGGCCG‑Tgcct                                                                                              >  2:253552/1‑248 (MQ=255)
                                                                                                                                                        ccTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCt                                                                                         <  2:244044/250‑1 (MQ=255)
                                                                                                                                                             cGATAGCCGCAAACTGCCTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCt                                                                                    <  2:71049/251‑1 (MQ=255)
                                                                                                                                                             cGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTGCACCCCATCATTGGTCATGCCGTGCACAAGAGGGGGATGTTGGTG‑AAAAAGGGGTAAAGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTTCAT‑AGTTTTCTCCTGCCATAAAGTGAAGAAGCAATGACCTGGGGGTCCG‑TGACTTCCCTAcccc                                                                                    >  1:76032/1‑248 (MQ=255)
                                                                                                                                                                 ggCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCa                                                                                <  1:10139/250‑1 (MQ=255)
                                                                                                                                                                 aGCCGCAAACTGACTGGCTTCTTCGGTTTCGATATCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAGGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCa                                                                                <  2:195956/251‑1 (MQ=255)
                                                                                                                                                                         acggactgGCCGCTTCGGCTTTGAGAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGG‑TTTGAACCACAGCATGGGTCATGCAGTGCCCAAGAGGGGAATGTGGGGG‑AAAAAGGTGTAACGCGTGCGCCGATTGCGCGGAACCCAGCAGTTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGCGGTCCG‑TGACTTCCCGACGCTGGCATTGTCCAg                                                                        <  2:354642/248‑1 (MQ=255)
                                                                                                                                                                           tgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAG‑GTTGCACCACATCCATGGTCATGCAGTGCACAAGAGGGGGATGTTGGGG‑AAAAAGGGGTAACGCGGGCGCAGATTTCGCTGAACCCCGCCGGGCGACTTGCAT‑AGTTTGCTCCTGCCCTAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACCCCGGCATTATcccagt                                                                      >  2:102404/1‑247 (MQ=255)
                                                                                                                                                                             actgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGGACCACATCATTTGGCATGCAGGGCCCAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCCGATTGCGCTGAACCCAGCAGGGCGACCTGCCT‑AGTTTGCTCCCGCCATAAAGTGGAGAAACAATGACCTGGTGGTCCG‑TGACTTCCCTACGCCGGCATTATCCAGATcc                                                                    >  1:371070/1‑250 (MQ=255)
                                                                                                                                                                                 gCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑GTTTCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGGG‑AAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGAATTGCAT‑AGGTTGGTCCTGCCATAACGGTAAGAAGCAAAGACCTGGGGGGCCG‑GGACTTCCCTACCCTGGCCATATCCAAATCAAGTg                                                                >  1:55734/1‑251 (MQ=255)
                                                                                                                                                                                        cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACggg                                                         <  1:318183/251‑1 (MQ=255)
                                                                                                                                                                                          gTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAGGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTa                                                       <  1:145999/251‑1 (MQ=255)
                                                                                                                                                                                              cGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGTTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTTCGCTGGCATTATCCAGATCAGGTGATACGGGTATTtc                                                   <  2:128170/251‑1 (MQ=255)
                                                                                                                                                                                              cGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTAttt                                                    <  1:22664/250‑1 (MQ=255)
                                                                                                                                                                                                aTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAG‑GTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTG‑AAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCCGGGCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGGCCG‑GGACTTCCCTCCGCTGGCCTTATCCAGAACAGGGGATACGGGGATTtctc                                                 >  2:40228/1‑251 (MQ=255)
                                                                                                                                                                                                    ccATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGcc                                             <  2:334554/251‑1 (MQ=255)
                                                                                                                                                                                                                          cgcCAGCAAGGTATTGGCGGTAAAGG‑GTTTCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTTCG‑GGACTTCCCTACGCTGGCATTATCCAGATCAGGGGATACGGGTATTTTTCAGCCTTCCCGCAGAAGGGCACCCCgg                       >  1:69480/1‑250 (MQ=255)
                                                                                                                                                                                                                             cAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTc                    <  1:120778/251‑1 (MQ=255)
                                                                                                                                                                                                                             cAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTc                    <  1:264770/251‑1 (MQ=255)
                                                                                                                                                                                                                                cAAGGCATTGGCGGTAAAGG‑GTTTCACCAAAACTTTGGTCCTCCAGTGGACAAGAGGGGGATGTTGGTG‑AAAAAGGTGTAAATCGTGCGCAGATTGCCCTGAAACCAGCCGGGCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGCAGAAGGCATGACCTGGTGGCCCGGTG‑CTTCCCTACCCTGGCCTTATCCCGATCAGGGGGAACGGGTCTTTTTCCGCCTTCACTCAGACGGGcccccccggtccgt                 >  2:69911/1‑243 (MQ=255)
                                                                                                                                                                                                                                 aaGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGttt                <  1:235334/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aaGGTATTGGCGGTAAAGG‑GTTGCACCACATCATTTGTCATGCAGTGCACCAGAGGGGGATGTTGGGG‑GAAAAGGTGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAAAGGGAAGAAGCAAAGAACCGGGGGGCCG‑GGAATTCCCCCCGCCGGCCTTATCCAGATCAGGGGGTTCGGGGAATTCTCAGCCTTCACGCAGAAAGGCCCCCCgcggcggtt                >  2:240389/1‑243 (MQ=255)
                                                                                                                                                                                                                                          gCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgat       <  2:229395/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATaa  <  2:400970/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aaaGG‑TTTGCACCACATCATTGGTCATGCAGAGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑CGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGAATTTCTCAGCCTTCACGCAGACGGGCACCCCGAGTCGCTTGGTTGCGATGATa   >  1:32387/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
AGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGG‑TTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTG‑AAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAT‑AGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCG‑TGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAA  >  NC_000913/2184955‑2185444

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