Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,874,4040AC29.1% 188.3 / 10.6 86V339G (GTG→GGG) cysNsulfate adenylyltransferase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (50/11);  new base C (2/23);  total (52/34)
Fisher's exact test for biased strand distribution p-value = 1.26e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAG  >  NC_000913/2874156‑2874648
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
cgatAAAAATCAGTCCACCCGTCACCGGATTTTGTTGATAACGATCTAACTCCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGTCAGGTTTTCAACTTAACGCTGGGTAAGGTTATTAATATCATCCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCAc                                                                                                                                                                                                                                                     <  1:515877/250‑1 (MQ=255)
      aaaTCAGCCCCCCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGAcgcg                                                                                                                                                                                                                                              <  2:536606/251‑1 (MQ=255)
              ccACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGGTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACCCTGCCCTGGAGAAAGCGGCCGTTCCGCCATCCATTCCACATCCACCCACgcgccccgcac                                                                                                                                                                                                                                      >  2:317979/1‑244 (MQ=255)
                aCCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCAccg                                                                                                                                                                                                                                    >  1:237446/1‑251 (MQ=255)
                  ccGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCAccgcc                                                                                                                                                                                                                                  <  1:29576/251‑1 (MQ=255)
                   cGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCAccgccg                                                                                                                                                                                                                                 >  1:708470/1‑251 (MQ=255)
                    gTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCCACCCGCTCTCCACCGCCgg                                                                                                                                                                                                                                >  1:441667/1‑251 (MQ=255)
                          ggATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGc                                                                                                                                                                                                                          >  1:712145/1‑251 (MQ=255)
                              tttGTTGATACCGATCTACCACCAGCGGCTGGTCAAAAGGGAGTCCCACGGGCCCAACCCCATTCAGTGGCAGGTTTTCACCTTCCCGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCACCACCACGCGCGCGCGTCTTCTTCCGGCCAATTTTGATGTCGTAACCCTCCCCTGGAGAAAGCGGCTGTCCCGCCACCCATCCCACACCCCCCGACGCCCTCTGCACCGCCGGTAACGCTTc                                                                                                                                                                                                                       <  2:601432/250‑1 (MQ=255)
                                       aaCGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTGAGTGGCAGGTTTTCAACTTCACGCTGGGGAAGGTTATTAATATCAACCTGATAGCGGATGCCCTCAACAACCGCGCGCGGCTTTTTACCGGCAATTTTTATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTTCGCCATCCCTACCACATCCCCCCACCCGGTCTGCCCCGCCGGGAACGCCTCGTCCgccgcc                                                                                                                                                                                                             >  2:460794/1‑251 (MQ=255)
                                         cGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGGTTTCAACTTCACGCTGGGTAAGGTTATTAATAACAACCTGATAGCGAATGCCACCAACACGCGCGCGGGTCTTCTTACCGGAAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTTCGCCAACCATACCACATCACCCGACGCCCTTCGCACCTCCGGGAACGCTTcacccgccgcccc                                                                                                                                                                                                           >  2:647948/1‑239 (MQ=255)
                                                  cccAGCGGCTCGTCAAAAGTGAGACCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCACCTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATc                                                                                                                                                                                                  <  1:649988/250‑1 (MQ=255)
                                                      gCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCa                                                                                                                                                                                              <  1:138744/251‑1 (MQ=255)
                                                        ggCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACg                                                                                                                                                                                            >  1:245548/1‑251 (MQ=255)
                                                              aaaaaagggggacccacgggcccGATCCCTTTAGGTGGCAGGTTTTCAACTTCAGGCTGGGTAAGGTTATTAATATCACCCTGATAGCGAATGCCCTCAACACGCGCGCGCGTCTTCTTCCCGGAAATTTTGATGTCGTACCTCGGCCGGGGAGAAAGGGGCTTTTCCGCCATCCATACCACCCCCCCCGCCGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGAt                                                                                                                                                                                      <  2:573729/230‑1 (MQ=255)
                                                               aaaaggTGAGACCCACGGTGCCGTTCCCATTGAGTGGCAGTTTTTCAATTTCACGCGGGGTAGGGTTTTTAATATCACCCTGATAGAGATGGCCATCAACACGCGCGCGCGTCTTCTTCCCGGAATTTTTGAGGTCGTAACTCTCCCCTGGAGAAAGCGGCTGTCCCCCCACCCATCCCACATCCCCCGACGCGCGCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATg                                                                                                                                                                                     <  2:256825/246‑1 (MQ=255)
                                                                aaaaGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCCCCAGCAGCAGAGCGCCACGGCTGATGt                                                                                                                                                                                    >  2:608362/1‑251 (MQ=255)
                                                                   aGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCg                                                                                                                                                                                  >  1:46768/1‑250 (MQ=255)
                                                                    gTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGAt                                                                                                                                                                                >  1:102780/1‑251 (MQ=255)
                                                                        gaTCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAAATTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGAtctcc                                                                                                                                                                            >  1:637078/1‑250 (MQ=255)
                                                                        aaTCCACGGGGCCGACCCCATTCAGTGGCAGTTTTTCACCTTCACGCTGGGTAATGTTATTAATATCACCCTGATAGCGAAGGCCATCAACGCGCGCGCGCGTCTTCTTGCCGGCATTTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCCCCACCCATACCCCACCCACCGCCGGGCTCTCCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCa                                                                                                                                                                            <  1:294686/250‑1 (MQ=255)
                                                                          cccACGAGGCCTACCCCTTTCAGTGGCAGTTTTTCAACTTCACGCGGGGTAGGGTTATTAATATCAACCTGATAGCGAATGCCAACAACCCGCGCGGGCGTCTTTTTCCCGGAAATTTTGATGTCGTAACTCTGCCCGGGAAAAGGCGGCTGTCCCGCCATCCACCCCACACCCACCGACGCGCTCGGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATc                                                                                                                                                                          <  2:205857/250‑1 (MQ=255)
                                                                           ccACGGGGCCTATCCCTTTTAGTGGCAGTTTTTAATCTTCCCGTGGGGTAAGGTTATTAATACCAACCTAATGGCGAATGCCATCAACACGCGCGCGCGTCTTCTTCCGGGCAGTTTTGTTGTCGTACCTCTGCCCTGGAAAAGGCGGCTGTCCCGCCCTCAATCCCACACCCCCCGACGCGTTCGGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCAt                                                                                                                                                                           <  2:584650/249‑1 (MQ=255)
                                                                           ccACGAGGCCGATCCCATTCAGGGGCAGGTTTTAAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCAGCTGTTCCGCCATCCATACCACATCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATcc                                                                                                                                                                         <  1:28301/251‑1 (MQ=255)
                                                                             aCGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCCCGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCCACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGGAACGCTTCGTCTGCCCCCCGCAGCAGATCGCCCCGGCTGATGTCGATCTCATCCgc                                                                                                                                                                       >  2:640827/1‑250 (MQ=255)
                                                                              ggAGGCCGTTCCCTTTCAGTGGCAGTTTTTCAACTTCACGTGGGGTAAGTTTATTAATATCAACCGGATAGCAAATGCCACCACCCCGCGCGCGCGTTTTCTTACCGGCAATTTTGATGTCGTACCTCTCCCTTGGAGAAAGCGGTTGTCCCGCCATCCATCCCCCACCCCCCGCCGCGCTCCGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTc                                                                                                                                                                      <  2:618846/250‑1 (MQ=255)
                                                                               gcggccgtgcccATTCAGTGGCAGGTTTTCAACGTCACGCTGGGTAAGGTTATTAATACCACCCTGATAGCGAATCCCATCAACACGCGCGCGCGTCTTCTTACGGGAAATTTTGATGTCGTAACTCTGCCCGGGAGAAAGCGGCTGTCCCGCCACCCATACCCCATCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCa                                                                                                                                                                     <  1:85389/242‑1 (MQ=255)
                                                                                  gCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCCGCAGCAGAt                                                                                                                                                                                                   >  1:229635/1‑218 (MQ=255)
                                                                                         ccATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTTCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCCgggg                                                                                                                                                           >  2:18908/1‑250 (MQ=255)
                                                                                                  ggCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGt                                                                                                                                                                                                                               >  1:504406/1‑174 (MQ=255)
                                                                                                  ggCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTtc                                                                                                                                                  >  1:437719/1‑251 (MQ=255)
                                                                                                  ggCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTtc                                                                                                                                                  >  1:27149/1‑251 (MQ=255)
                                                                                                                    cGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCttcttc                                                                                                                                >  1:576694/1‑251 (MQ=255)
                                                                                                                    cGCTGGGTAAGGTTATTAATACCAACCTGATAGCGAATGCCACCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCGGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATCCCACACCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCttcttc                                                                                                                                <  1:640827/251‑1 (MQ=255)
                                                                                                                        gggTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgcg                                                                                                                             >  1:684323/1‑250 (MQ=255)
                                                                                                                         ggTAAGGTTATTAATAACAACCTGATAACGAAAGCCAACACCCGGCGCGCGCGTCTTCTAACCGGCAATTTCGTTGTCGTAACTCTCCCCTGGAGAACGGGCCTGTCCCCCCACCCACACCACATCCACCGCCGCCTTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGAt                                                                                                                           <  1:269964/251‑1 (MQ=255)
                                                                                                                            aagtttattaATACCAACCGGAAAGAGAAGGCCACCAACCGGGGCGCGCTCCTTCTTACCGACAATTTTGATTGCGTACCTCGCCCCGGGAGAAGGCGGTTGTCCCCCCATCAATACCACACCCCCCGACGCGCTCGGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCAc                                                                                                                        <  2:88816/247‑1 (MQ=255)
                                                                                                                            aaGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCAc                                                                                                                        >  1:699154/1‑251 (MQ=255)
                                                                                                                            aaGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCAc                                                                                                                        >  1:656466/1‑251 (MQ=255)
                                                                                                                                ttattaATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATc                                                                                                                    >  2:49744/1‑251 (MQ=255)
                                                                                                                                       tatCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttct                                                                                                                                  >  1:649715/1‑230 (MQ=255)
                                                                                                                                        atccaCCTGATAGCGATGGCCATCACCCCGCGCCCGCGTCTTCTTCCCGGTATTTTTTATATATTCACTCGCCCCTGGAAAACGCTGGTGTTCCCCCCCCCATCCCACACCCCCCGCCGCGCTCTCCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGAGGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTc                                                                                                             <  2:412297/246‑1 (MQ=255)
                                                                                                                                        atCACCCTGATAGCGAATGCCATCAACCCGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATCCCACACCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCa                                                                                                            <  1:137180/251‑1 (MQ=255)
                                                                                                                                        atCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGc                                                                                                                                                                   >  1:469591/1‑196 (MQ=255)
                                                                                                                                           aaCCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACg                                                                                                          >  1:167658/1‑250 (MQ=255)
                                                                                                                                             ccggATAGCGAAGGCCATCACCACGCGCGGGCGTCTTCTTACCGGAAATTTTGATGTCTTACCGCTCCCCGGGAGAAGGCGGTTTTTCCCCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATc                                                                                                       <  1:679707/248‑1 (MQ=255)
                                                                                                                                                gATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGCCACGATctgc                                                                                                    >  2:654100/1‑248 (MQ=255)
                                                                                                                                                gATAGCGAATCCCACCACCCCGCGCGCGCGTCTTCTTACCGGAATTTTTGAGGCCGTAACCCGGCCCGGGAAAACGCGGCTGTCCCGCCCTCCATCCCACACCCCCCGCCGCGCTCCGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATcc                                                                                                      <  2:46768/249‑1 (MQ=255)
                                                                                                                                                 aTAGCGAATGCCACCAACACGCCCGCGCGCCTTCTTACCGGCAATTTTGATGTCGTACCTCTGCCCTGGAGAAAGCGGCTGTTCCCCCACCCATACCACACCCCCCGACGCGCTCGGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcg                                                                                                   <  1:680856/251‑1 (MQ=255)
                                                                                                                                                  tAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgc                                                                                                  >  1:591778/1‑251 (MQ=255)
                                                                                                                                                    gCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCCGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgcg                                                                                                 >  1:710950/1‑250 (MQ=255)
                                                                                                                                                      gAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGCTCTCATCCGTCAGCACCAGGGCGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACg                                                                                              >  1:117254/1‑251 (MQ=255)
                                                                                                                                                         tGCCATCAACACGCGCGCGAGTCTTCTTACCGGCACTTTTGATGCCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGACATCCATACAACATCCACCGACGCGCTCTGCACCCCCGGTAACGATTCGTCTGCCGCCAGCAGCAGATCGCCACGCCTGAt                                                                                                                                                                                      >  1:424988/1‑160 (MQ=255)
                                                                                                                                                               cAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGGGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGGTTGATTcc                                                                                     >  2:237471/1‑251 (MQ=255)
                                                                                                                                                               cAACACGCGCGCGCGTCTTATTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTcc                                                                                     >  2:307950/1‑251 (MQ=255)
                                                                                                                                                                 acacGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCcac                                                                                   >  1:187484/1‑251 (MQ=255)
                                                                                                                                                                            cgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGgc                                                                         >  1:370364/1‑250 (MQ=255)
                                                                                                                                                                            cgTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGACTCCACACCAGAGGgc                                                                         >  1:676609/1‑250 (MQ=255)
                                                                                                                                                                               cttcttACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGCAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTCACGCTTCGTCTGCCGCCAGCACCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGgcagc                                                                      >  2:412980/1‑250 (MQ=255)
                                                                                                                                                                                 tcttCCCGACAATTTTGATGTCGTAACTCTCCCCTGGAGAAAGCGGCTGTCCCGCCATCCATCCCACACCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCAc                                                                    <  1:339485/250‑1 (MQ=255)
                                                                                                                                                                                    tcccGGCAATTTTGATGTCGTAACTCTGCCCGGGAGAAAGCGGCTGTTCCGCCATCCATCCCCCACCCCCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACt                                                                   <  1:460794/246‑1 (MQ=255)
                                                                                                                                                                                     tacaggaaaTTTTGAGGCCGTACCTCTGCCCTGGAGAAAGCGGCTGTTCCGCCACCCATCCCCCATCCACCGCCGCGCCCCGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACttt                                                                 <  2:710950/243‑1 (MQ=255)
                                                                                                                                                                                        ggCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGacac                                                            >  2:21422/1‑251 (MQ=255)
                                                                                                                                                                                                ttgttgacGTAACTCCGCCCGGGAAAAAGGGCTTGCTCCCCCATCCATACCACATCCCCGGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGcc                                                      <  1:719749/243‑1 (MQ=255)
                                                                                                                                                                                                ttGATGTCGTACCTCTGCCCTGGAAAAAGCGGCTGTTCCGCCATCCATCCCACACCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgc                                                                                                    <  1:546381/204‑1 (MQ=255)
                                                                                                                                                                                                ttGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgc                                                                                                    >  2:546381/1‑204 (MQ=255)
                                                                                                                                                                                                 tGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCg                                                    >  1:547894/1‑251 (MQ=255)
                                                                                                                                                                                                  gATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCg                                                    >  1:711219/1‑250 (MQ=255)
                                                                                                                                                                                                    tGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACt                                                 >  1:98467/1‑251 (MQ=255)
                                                                                                                                                                                                           aCTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGt                                                          >  1:660118/1‑235 (MQ=255)
                                                                                                                                                                                                             tctGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGCTGCCCGACTTCCAcgcg                                         >  2:493300/1‑250 (MQ=255)
                                                                                                                                                                                                              ctGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGAc                                       >  2:618464/1‑251 (MQ=255)
                                                                                                                                                                                                               cgCCCGGGAGAAGGCGCCTGTCCCCCCCTCCATACCCCACCCCCCGCCGCGCCCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGAcc                                      <  2:148433/250‑1 (MQ=255)
                                                                                                                                                                                                                ccccTGGAGAAAGCGGCTGTCCCGCCATCCATACCCCATCCCCCGCCGCCCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGAc                                       <  1:49744/248‑1 (MQ=255)
                                                                                                                                                                                                                   cTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGAt                                  >  1:500325/1‑251 (MQ=255)
                                                                                                                                                                                                                      gagaAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGcc                               >  2:481961/1‑251 (MQ=255)
                                                                                                                                                                                                                      gacaatgCGGCTCTCCCGCCACCCATCCCACACCCCCCCACGCGCTCACCCCCGCCGGTAACGCTTCGTCTTCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGcc                               <  2:367029/245‑1 (MQ=255)
                                                                                                                                                                                                                       agaAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGcc                               >  2:295912/1‑250 (MQ=255)
                                                                                                                                                                                                                         acaccGGCTGTCCCGCCAGCCACCCCTCCTCCCCCGCCGCGCTCGCCACCGCCGGTAACTCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGGCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGc                            <  1:400680/247‑1 (MQ=255)
                                                                                                                                                                                                                            gctgCTGTCCCGCCATCCATGCCACACCCACCGACGCGCTCTGCGCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGtt                         <  1:307950/248‑1 (MQ=255)
                                                                                                                                                                                                                              ggCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTcc                       >  2:643480/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGc                    >  1:476647/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGAt                                                                                                        >  1:349511/1‑161 (MQ=255)
                                                                                                                                                                                                                                            tccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACg           >  1:310242/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 cccACATCCCCCGACGCGCTCTGCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATc     <  1:237471/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccacccACCGCCGCGCTCCGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCgag  <  1:608362/247‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
CGATAAAAATCAGCCCACCCGTCACCGGATTTTGTTGATAACGATCTAACACCAGCGGCTCGTCAAAAGTGAGATCCACGAGGCCGATCCCATTCAGTGGCAGGTTTTCAACTTCACGCTGGGTAAGGTTATTAATATCAACCTGATAGCGAATGCCATCAACACGCGCGCGCGTCTTCTTACCGGCAATTTTGATGTCGTAACTCTGCCCTGGAGAAAGCGGCTGTTCCGCCATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACCAGAGGGCAGCACTTTGACACGTTGCCCGACTTCCACGCGACCGGATGCCAGCGTTCCGGCGTAACCACGAAAATCGAG  >  NC_000913/2874156‑2874648

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