Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,901,8170C.35.6% 138.1 / 103.0 75coding (404/1617 nt)bglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (20/27);  new base . (13/13);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-02
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTTTGGTGGTAACGTACATATCGGAGAACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCGTCGCTGGCGTATACTTTTTCTTCTCTTCATCT  >  NC_000913/3901567‑3902065
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
aCCTTTGGTGGCCACGTACATATCAGAAACCTGCATATAGATCTCGCCCCCGACGTTATCGCCAAACCAGCAGGTAGAGTACTGATGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCCAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGACCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTccccc                                                                                                                                                                                                                                                          >  2:308389/1‑251 (MQ=255)
 ccTTTGGTGGTAACGCACATATCGGAGAACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAAATGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcccccc                                                                                                                                                                                                                                                         >  2:220309/1‑250 (MQ=255)
     gggtggtAACGTACATATCGGAGAACGGCATATAGTTCCCGCCCCCGGCGTTAGCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTGCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAg                                                                                                                                                                                                                                                    <  2:421019/250‑1 (MQ=255)
          gtAACGTACATATCGGAGAACTGAATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGCAGAGTACTGCTGACCAACGTTAACATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATAGCCCATGATTTAGGAGGACCATGAGATGCGGTCCCCAAACCATAACg                                                                                                                                                                                                                                                >  2:507167/1‑251 (MQ=255)
                     atCGGAGAACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGG‑GCCCCAACCAGAACGGTAATAGCCGtt                                                                                                                                                                                                                                    >  2:442994/1‑251 (MQ=255)
                       cGGATAACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTGTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGATCCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTg                                                                                                                                                                                                                                   <  1:250017/250‑1 (MQ=255)
                        ggAGAACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGt                                                                                                                                                                                                                                  >  1:174227/1‑250 (MQ=255)
                         gagaACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGgtgtt                                                                                                                                                                                                                                 >  2:9053/1‑250 (MQ=255)
                         aagaACTGCTTATATTGCGCGCCACCGGCGTTAGCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGGAATAGCCGTTGTa                                                                                                                                                                                                                                 <  2:236240/250‑1 (MQ=255)
                          agaACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGCACGGTAATAGCCGTTGTa                                                                                                                                                                                                                                 >  2:232897/1‑250 (MQ=255)
                          agaACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTaa                                                                                                                                                                                                                                <  1:122924/251‑1 (MQ=255)
                             aCTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTa                                                                                                                                                                                                                             <  1:566318/251‑1 (MQ=255)
                             aCTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTaa                                                                                                                                                                                                                            <  1:249567/251‑1 (MQ=255)
                                 cATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATCCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCCGAACGGTAATAGCCGTTGTAACTAAAGc                                                                                                                                                                                                                         >  2:479079/1‑251 (MQ=255)
                                    ataGTTCTCGCCCCCGGCTTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCTTGGTAGACACGTTGTTTTAACTGCAAATCAAACCGGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAAGCACCCATGATTGAGGTTAACCATGAGAGGCGGTCCCCCAACCATATCGGTAATAGCCGTTGTAACTAAAGCCa                                                                                                                                                                                                                       <  2:602953/250‑1 (MQ=255)
                                      aGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAt                                                                                                                                                                                                                     <  1:120458/250‑1 (MQ=255)
                                            cGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCt                                                                                                                                                                                                              <  1:181155/251‑1 (MQ=255)
                                             gCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTCCTGACCAACGATACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCtt                                                                                                                                                                                                             >  1:102407/1‑251 (MQ=255)
                                                      gTTATCGCCAACCCAGCCGGTAGAGTACTGCTGACCAAAGTTACCATCCATCATCACAACGGCATCAACCCGTTTGACGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGACGGAGTATTCGTTACCAGAGCGTCCCAGAGACCCACTCGCCCATGATTTAGGGGAAACATGAGCGGCGGTCCCACCACCAGAACGGCAATAGCCGTTGTAACCAAAGCCAATTTCCTCTTACACaaaat                                                                                                                                                                                                    >  2:56428/1‑249 (MQ=255)
                                                        tATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAttt                                                                                                                                                                                                   <  2:242623/250‑1 (MQ=255)
                                                         aTCGCCAAACCAGCCGGTAGAGGACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACTTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGAGTTAGGTGAACCATGATAGTCGG‑GCCCCAACCAGAACTGTAATAGTCTGTGTAACTAAAGCCAATTTCATCTTTCACAAATTTAc                                                                                                                                                                                                 <  1:354773/250‑1 (MQ=255)
                                                            gCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGc                                                                                                                                                                                                                                        <  1:386166/208‑1 (MQ=255)
                                                            gCCAAACCAGCCGGTAGAGAACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGc                                                                                                                                                                                                                                        >  2:386166/1‑208 (MQ=255)
                                                             ccAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaa                                                                                                                                                                                             <  1:245176/251‑1 (MQ=255)
                                                              cAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAAt                                                                                                                                                                                           <  2:333210/251‑1 (MQ=255)
                                                                    aGCCGGTTGAGTACTGCTGACCAACGTTTCCATCCTTCATCACAACGGCATCAACCCGTTTGCTGTTTTCGTTGTTGACACGTTTTTTTAACTGCAAATCAAACCAGCTGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGACCCATTAGAGGCGGTCCCCCAACCATAACGGTAATAGCCGTTGTATCTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTc                                                                                                                                                                                      <  2:528972/251‑1 (MQ=255)
                                                                    aGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTc                                                                                                                                                                                      >  1:110563/1‑251 (MQ=255)
                                                                           agattACTGCTTACCAACGTTCCCATCCATCATCACAACGGCATCAACCCTTTTGCCGGTTTCGTTGTAGCCCCGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAGTTCATCTTTCACAAATTTACTGAAATCTTTCAGTGTc                                                                                                                                                                                <  2:179762/246‑1 (MQ=255)
                                                                           agagTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGTTCCCCCAACCAGAACGGTAATAGCCGGTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTc                                                                                                                                                                                <  1:271817/250‑1 (MQ=255)
                                                                           agagTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTc                                                                                                                                                                                <  1:9053/250‑1 (MQ=255)
                                                                           agagTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGGAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATTTTTCAGAGTcc                                                                                                                                                                               >  2:169933/1‑250 (MQ=255)
                                                                             agTACTGCTGACCAACGTTCCCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGTCCCAGAGAACCAATCGCCCATTATTTAGGTGAACCATGATAGGCGGTCCCCCATCCAGAACGGTAATAGCCTTTGTAACTAAAGCCAAGTTCATCTTTCACAAATTTACTGGAATCTTTCAGAGTCATg                                                                                                                                                                             <  2:597008/251‑1 (MQ=255)
                                                                              gTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTTAACCAGGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATgg                                                                                                                                                                            <  1:584712/251‑1 (MQ=255)
                                                                               tACTGCCGACCAACGTTACCATCCATAATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGCCACGTTGTTTTAAATGCAAATCACAACAGCCGGAGTATTAGTTCCCAAAGCGGACAAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGG‑GACCAAACCAGAACGGTAATAGCCGATGGAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCagaggactgga                                                                                                                                                                           >  2:63086/1‑243 (MQ=255)
                                                                                         cAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATc                                                                                                                                                                 <  1:506940/251‑1 (MQ=255)
                                                                                         cAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCg                                                                                                                                                                >  2:121160/1‑251 (MQ=255)
                                                                                                aCCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCCGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCCTCTTTCACAAATTTACTGAAATTTTTCAGAGTCATGGAAGAATAAATCGAGCCtt                                                                                                                                                          >  2:77897/1‑250 (MQ=255)
                                                                                                  catccatcATCACAACGGCATCAACCCGTTTGCCGTGTTCGTTGTATACACGTTGTTTTAACTGCAAATCAATCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGAGTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTTCTGAAATCTTTCAGAGTCATGGAAGAAGAAATCTAGCCTGtc                                                                                                                                                        <  2:213612/251‑1 (MQ=255)
                                                                                                      catcatcttagcGGCATCCACCCGTTTGCCGGGTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAACCCAGCCGGAATATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGGGGTCCCCCAACCAGAACGCTAAAAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAACTGACTGAAATCCATAAGAGTCATGGAAGAATAAATCGAGCCTGTCTCAc                                                                                                                                                    <  1:56428/240‑1 (MQ=255)
                                                                                                       atcatcACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AAACAGAACGGTAATAGCCGTTGTAACTAAAGCCACTTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTg                                                                                                                                                  >  1:512935/1‑251 (MQ=255)
                                                                                                            cacaAAGGCATCAACCCGGTTGACGTTTTCGTTGTAGACACGTTGTTTTAACTGAAAATCAAACCCGCCGGTGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGG‑GACCCAACCAGAACGGTAATAGCCGATGTAACTAAAGCCAATTTCATCTTTAACAAATTTACTGAAATATTTCAGAGTCATGGACGAATAAATCGAGCCTGTCTCACTGGcctt                                                                                                                                             >  2:221227/1‑248 (MQ=255)
                                                                                                              caACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCAtttt                                                                                                                                           >  1:345366/1‑251 (MQ=255)
                                                                                                               aaCGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAACCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCCGAACGGTAATAGCCGTTGTAACTAAAGACAATTTCATCTTTCACAAATTTACTGCAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCCCTGGCAttttt                                                                                                                                          >  2:240286/1‑251 (MQ=255)
                                                                                                                           ccGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCAtttt                                                                                                                               <  1:77897/251‑1 (MQ=255)
                                                                                                                           ccGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAAAATAAATCGACCCTGTCTCACTGGCATTTTTCTCTTCAttttt                                                                                                                              >  1:107831/1‑251 (MQ=255)
                                                                                                                                 gCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCATAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACGGGCATTTTTCTCTTCATTTTTGAcc                                                                                                                          <  2:54763/250‑1 (MQ=255)
                                                                                                                                  ccGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGt                                                                                                                       <  1:240286/251‑1 (MQ=255)
                                                                                                                                      tttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCTTTTTTCTCTTCATTTTTGACCAGTAcc                                                                                                                    >  2:3954/1‑250 (MQ=255)
                                                                                                                                      tttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACa                                                                                                                    <  2:174227/250‑1 (MQ=255)
                                                                                                                                               aGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAgg                                                                                                            <  1:111785/250‑1 (MQ=255)
                                                                                                                                                     gttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTCCCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTTTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGct                                                                                                     <  2:341628/251‑2 (MQ=255)
                                                                                                                                                            tAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGGAAAAGCCGTTGTAACTAAAGCCAATTTCCTCTTTCACAAATTTACTGAAATCTTTCCGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTTACCAGTACAGCATCAGGTTTTGCGGCactac                                                                                              >  2:428841/1‑251 (MQ=255)
                                                                                                                                                            aaaCTCCAAATCAATCCAGCCGAGGTATTCGTTACCAAACCGGCCCAGAGACCCAATCGCCCAAGATTTAGGTGAACAATGAGGGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACGAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCactac                                                                                              <  2:86898/249‑1 (MQ=255)
                                                                                                                                                              aCTGCAAATCAACCCAGCCGGAGTATTCGGTCCCAAAGCGGCCCAGAGAACCAATCGCCAATGATTTAGGTGAACCATGAGTGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATGTCATCTGTCACAAATGTACTGAAATCTTGCAGAGTCATGGAAGAATAAATCTAGCCTGTCTCACTGGCATGGTTCTCTTCATTTTTGACCAGGACAGCATCAGGGTTTGAGTCACTACTgg                                                                                           <  2:576164/251‑1 (MQ=255)
                                                                                                                                                               cTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGAGGCCCAGAGAACCAATCGCCCATGATTTAGGCGAACCATGAGAGGCGGTCCCC‑AACCCGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGGGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTAAAGCATCAGGTTTTGCGGCACTACTGGt                                                                                          >  1:397091/1‑251 (MQ=255)
                                                                                                                                                                tGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTTTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTc                                                                                         >  1:217253/1‑251 (MQ=255)
                                                                                                                                                                 gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGAGTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTGTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTc                                                                                         <  2:55508/251‑1 (MQ=255)
                                                                                                                                                                   aaaTCAAACCAGCCGGAGTATTCGTTACCCAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCCTGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCCGTACAGCATCAGGTTTTGCGGCCCTCCTGGTCgg                                                                                       >  2:265970/1‑251 (MQ=255)
                                                                                                                                                                       cAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAAACCATAGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCATGCCTGTGTCACTGGCATTTTTCTCTTCCTTTTTGACCAGTAACGCCTCAGGTTTTGCGGCACTACTGGTCGGGAAc                                                                                   >  1:546294/1‑251 (MQ=255)
                                                                                                                                                                               ccGGGGTATCCGTTACCAAAGCGGCCCAGAGAACCAATCTCCCACGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAAAACGGTAAGAGCCGTTGTAACTAAAGCCAATTTCATGTTTCACAATTTTACTTAAATCGTTCAGTGTCATGGAAGAATAAAGCGAGGCTGTCTCACTGGCATTTTGCTCTTCATTTTGGACCAGTACAGCATCAGGTTTTGCTGCGCTTTTGGTCGGCAACGGATTggc                                                                           <  2:62483/251‑1 (MQ=255)
                                                                                                                                                                                      aTTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCAt                                                                                                                                  >  2:317837/1‑188 (MQ=255)
                                                                                                                                                                                      aTTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCAt                                                                                                                                  <  1:317837/188‑1 (MQ=255)
                                                                                                                                                                                         cGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGAGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAAATAAAGCCAATTTCATCTTTCCCAAATTTACTGAAATCTTTCAGAGTCATGGGAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGGTTTGCGGCACTACTGGGCGGGAACGGATTGGCGGCATACCCt                                                                 >  2:51572/1‑251 (MQ=255)
                                                                                                                                                                                                    cGGCCCATAGACCCACTCACCCATGATTGAGGTGAACCATGTGAGCCGGTCACCCATCCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGATTCATGCAAGAATAAATCGAGCCTGTCTCACTGGCATGTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTGGATCATTCGt                                                      <  1:399959/251‑1 (MQ=255)
                                                                                                                                                                                                       cccAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATATTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCCTTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACtt                                                  >  1:195722/1‑251 (MQ=255)
                                                                                                                                                                                                          agagAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGTTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCGTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCAGCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTAc                                                <  1:220309/251‑1 (MQ=255)
                                                                                                                                                                                                          agagAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAAATAAAGCCAATTTCATCTTTCACAAATTTCCTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCATGTCTCACTGGCATTTTTCTCTTCATTTTTGGCAAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAAAGGATTGGCGGCATACCCTTGATCATTCGGACTTAc                                                >  1:74663/1‑251 (MQ=255)
                                                                                                                                                                                                           gagaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGc                                              >  1:436413/1‑251 (MQ=255)
                                                                                                                                                                                                            agaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGc                                              >  2:433564/1‑251 (MQ=255)
                                                                                                                                                                                                             gaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCAt                                                          >  2:429712/1‑237 (MQ=255)
                                                                                                                                                                                                             gaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAAGCGAGCCTGTCTCACTGGCATTTTTCTCGTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCAt                                                          <  1:429712/237‑1 (MQ=255)
                                                                                                                                                                                                              aaCCAATCGCCCATGATTTAGGTGAACCATGGGAGGCGGTCCCC‑AACCAGAACGGTAATAGCCGTTTTAACTAAAGCCAATTTCAGCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGGCCAGTACAGCATCAGGTATTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTAc                                           <  1:598571/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGGAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCg                                  >  1:234149/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCGTc                                <  1:232897/251‑1 (MQ=255)
                                                                                                                                                                                                                             aTTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAAAGGATTGGCGGCATAACCTTGATCATTCGTCCTTCCGCCACGATTCACCggccct                             >  2:346088/1‑246 (MQ=255)
                                                                                                                                                                                                                                                 gCGGTCCCC‑AACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCGTCGCTGGCGTATACTTTTtcttctct        <  1:90699/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ggTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCGTCGCTGGCGTATACTTTTtcttctctt       >  1:521878/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       ccccAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCGTCGCTGGCGTATACTTTTTCTTCTCTTCATc   <  1:265970/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        cccAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCGTCGCTGGCGTATACTTTTTCTTCTCTTCATCt  >  1:456897/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
ACCTTTGGTGGTAACGTACATATCGGAGAACTGCATATAGTTCTCGCCACCGGCGTTATCGCCAAACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTACAGCATCAGGTTTTGCGGCACTACTGGTCGGGAACGGATTGGCGGCATACCCTTGATCATTCGTACTTACGCTACGATTCACCGTCGCTGGCGTATACTTTTTCTTCTCTTCATCT  >  NC_000913/3901567‑3902065

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