Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,745,0890AC34.7% 138.7 / 18.3 77V84G (GTG→GGG) trmDtRNA m(1)G37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/16);  new base C (1/25);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 6.65e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCAC  >  NC_000913/2744843‑2745337
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
gCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAGTCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCt                                                                                                                                                                                                                                                       <  2:556299/250‑1 (MQ=255)
gCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCtt                                                                                                                                                                                                                                                      <  2:151762/251‑1 (MQ=255)
 ccGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCttt                                                                                                                                                                                                                                                     <  1:17984/251‑1 (MQ=255)
 ccGAGGCATCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCttt                                                                                                                                                                                                                                                     <  1:43251/251‑1 (MQ=255)
  cGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTg                                                                                                                                                                                                                                                    <  2:227402/251‑1 (MQ=255)
    agtcTTCATGTCCCCGTACTCCCGGAATAAACCGGGATACGGAGTCAATCAGCGTCATTGCTGGTAACTAACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGCCGCCCGCTTGATCCAGCGTGCGTCCCTGTGCGGACAGATAAATCACCTTTgcg                                                                                                                                                                                                                                                  >  2:357039/4‑251 (MQ=255)
      gCTTCATGTCCCAGTACTCCCGGAATAAACTGGGAAACGGAGTCAATCAGCGTCATTGCTGCTAACTCACCACCACTGAGAACGTAATCGCCGAATGACCATTATTCGTCAATTTCGGATTGGATCACGCGCTCATCTATACCTTCGTAGCCACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTAGCGGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGAGACAGATAAATCCCCTTTGCGcc                                                                                                                                                                                                                                                >  2:458326/1‑251 (MQ=255)
          cATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTc                                                                                                                                                                                                                                             >  2:78499/1‑250 (MQ=255)
            tGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCAcc                                                                                                                                                                                                                                          >  2:106627/1‑251 (MQ=255)
             gTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCAccc                                                                                                                                                                                                                                         <  1:509466/251‑1 (MQ=255)
             gTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCAcc                                                                                                                                                                                                                                          <  1:185220/250‑1 (MQ=255)
                                     gggAAACGTAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCCCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGTTTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTTACAAATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATg                                                                                                                                                                                                                 <  1:447096/251‑1 (MQ=255)
                                         aaCGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCACGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGt                                                                                                                                                                                                             >  2:442248/1‑251 (MQ=255)
                                          aCGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTc                                                                                                                                                                                                            >  1:94899/1‑251 (MQ=255)
                                                     cagcgcattggCTGGTAACCCCCCCCCCCTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGAATCACGCGCTAATCTACACCTTCGTAGCGACCGCCCACCAAATTCAATTTTGGATTCGTGGCCAGTTCGCTGACGCCCGCTTGATAAAGCTGGCGTCCCTGGGGTGACAGAAAAACCCCCTTGGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                                                 <  2:253296/241‑1 (MQ=255)
                                                     cAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTTTTCTTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTTCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGtt                                                                                                                                                                                                 <  1:285042/251‑1 (MQ=255)
                                                      aGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTg                                                                                                                                                                                                >  1:115482/1‑251 (MQ=255)
                                                       gCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCACTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTCCATGAATGCCGTCCCGCAAGGGttcc                                                                                                                                                                                               >  2:37468/1‑249 (MQ=255)
                                                        cGTCATTGCTGGTAACTCCCCACCCCTGAGAACGTAATCCCCGATTTCCCATTCTGCGTCAATTTCGGTTTGGACCACGCGCTCTTCTATACCTTCGTAGCTCCCGCCCCCCAAAATCAATTTTTGATCTGTTGCGATTTCGCTGACGCCCGCTTTATCAAGCTTGCGTCCCTGTGGTGACATATAAATCGCCTTTGCCCCTTCTCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCa                                                                                                                                                                                              <  2:124768/251‑1 (MQ=255)
                                                         gTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCCAGGGTTGCAc                                                                                                                                                                                             >  2:55723/1‑251 (MQ=255)
                                                            aTTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcat                                                                                                                                                                                          >  1:6963/1‑251 (MQ=255)
                                                              tGCTGGTAACCCACCACCACTGAGAACGTAACCGCCGATTGACCTTTCTTCGTCAATTTCGGTTTGGATCCCGCGCTCATCTATCCCTTCTTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTGGCCGGTTCGCTGACGCCCGCTTGATCAAGCTGGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCTCcatc                                                                                                                                                                                         <  2:440183/250‑1 (MQ=255)
                                                                      ccgcccccccacggaaaactgaaccGCCGATGGACCATTCTTCGTCAATGTCGGTTTGGATCCCGCGCTCATCTATCCCTTCGTAGCGCCCGCACCCCAGAAACAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATc                                                                                                                                                                                <  2:357339/227‑1 (MQ=255)
                                                                        ccaccaccacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCTCCCGCACACCAGAATCAATTTTTGATTCGTTGCCGGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATccc                                                                                                                                                                              <  2:487018/250‑1 (MQ=255)
                                                                         cccccccacTGAGACCGTAACCGCCGATTGCCCATTCTTCGTCAATTTCGGTTTGGCTCACGCGCTCATCTATCCCTGCGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGTTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATccc                                                                                                                                                                              <  2:339328/245‑1 (MQ=255)
                                                                              ccacTGAGAACTTAATCCCCAATTGACCATTCTTCTTCAATTTCGGTTTGAATCACGCCCTCATCTATACCTTCGTAGCGCCCCCACACCAAAATCAATTTTTTATTCTTTGCCAGTTCGCTTACCCCCGCTTTATCAAGCTTTCGTCCCTTTGGTGACAAATAAATCCCCTTTTCCCCTTCACCCGCCGCCGCTTTTGCTTCATGAATGGCTTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTcc                                                                                                                                                                        <  1:364897/251‑1 (MQ=255)
                                                                              ccacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTcc                                                                                                                                                                        >  1:188044/1‑251 (MQ=255)
                                                                              ccacTGAGAAAGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTcc                                                                                                                                                                        >  1:544023/1‑251 (MQ=255)
                                                                               cacTGAGAACGTAATCGCCTATTGCCCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccg                                                                                                                                                                       <  2:68025/251‑1 (MQ=255)
                                                                               cacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccg                                                                                                                                                                       >  1:482153/1‑251 (MQ=255)
                                                                                          taaccGCCAATTGCCCATTCTCGGTCATTTTCGTTTTGGATCCCGCCCTCATCTATCCCTTCGTAGCGCCCGCCCACCAGAATAAATTTTTGATTCGTTCCCAGTTCGCTGCCGCCCGCTTGATCAAGCTGGCGCCCCTGTGGGGACAGATAATCCCCCTTGGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTaa                                                                                                                                                              <  2:479719/245‑1 (MQ=255)
                                                                                                       cccATTCTCCGTAAATTTCGGTTTGGAACCCGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAACCAATTTTGGATTCGTTGCCGGTTCGCTGACGCCCGTTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCCCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacg                                                                                                                                               <  2:536065/250‑1 (MQ=255)
                                                                                                               tcGTCAATTTCGGTTTGGATCACGCGCTCATCTAGACCTTCGTAGCGACCGCACACCAGAATAAATTTTGGATTCGTTGCCAGTTCGCTGACGCCCGCTTGAACAAGCTTGCGTCCTGGGGGGGACAGAAAAATCACCTTGGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggtacg                                                                                                                                          <  1:281613/248‑1 (MQ=255)
                                                                                                                cGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGc                                                                                                                                      >  2:328029/1‑251 (MQ=255)
                                                                                                                cGTCAATTTCGGTTTGGATCACCCCCTCACCTATCCCTTCGTAGCGACCCCACACCAGAATCAATTTTTGATGCGTTGCCAGTTCGCTACCGCCCGCTTGATCAATCTTGCGTCCCTGTGCGGACAGAAAAATCTCCTTGGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGATTGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACTGTACGGTg                                                                                                                                       <  1:276627/250‑1 (MQ=255)
                                                                                                                    aaTTGCGGTTTGGATCACCCGCTCAACTATACCTTCTTAGCGACCCCCCACAAAAACCAATTTTTGATTCTTTGCCAGTCCGCTGACCCCCGTTTGATAAAGCTTGCGTCCCTGGGGTAAAAAAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGt                                                                                                                                  <  2:243493/251‑1 (MQ=255)
                                                                                                                              tGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGa                                                                                                                         >  2:541697/1‑250 (MQ=255)
                                                                                                                                 aTCACGCGCTCTCCTATACCTTCGTAGCGCCCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGCCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTc                                                                                                                     <  1:46907/251‑1 (MQ=255)
                                                                                                                                   cACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgc                                                                                                                   >  1:65767/1‑251 (MQ=255)
                                                                                                                                      gcgcTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTcgcg                                                                                                                  >  2:269925/1‑249 (MQ=255)
                                                                                                                                            aTCTATCCTTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCATTTCGCTGACGCCCGCTTGATCAAGCTTGCTTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCCCCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTc                                                                                                           <  1:541902/250‑1 (MQ=255)
                                                                                                                                             tctatcccttggtcgcacccGCCCACCATAATCAATTTTTGTTTCGTTGACGGTTCGCTGACGCCCGCTGGATAAAGTTGGCGCCCTGGGGGTGACAGAAAAATCCCCTGTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACAGCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTcc                                                                                                          <  2:295827/232‑1 (MQ=255)
                                                                                                                                               tataCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGc                                                                                                       >  2:155703/1‑251 (MQ=255)
                                                                                                                                                   ccTTCGTAGCGCCCGCCCACCAGAATCAATTTTGGATTCGTTGCCAGTCCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTg                                                                                                   <  1:418190/251‑1 (MQ=255)
                                                                                                                                                        gTAGCGACCGCACACCAGAATCAATTTTTGATTCGTGGCCAGTTCGCTGCCGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                                              <  1:55723/251‑1 (MQ=255)
                                                                                                                                                          aGCGACCGCACCCCAGAATCAATTTTTGATTCGTTCCCAGTCCGTTGACCCCCGCTTGACCAAGCTTGCTCCCCTGTGGTGACAGATAAACCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCt                                                                                             <  1:456404/250‑1 (MQ=255)
                                                                                                                                                             gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                                              >  2:175312/1‑246 (MQ=255)
                                                                                                                                                             gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGGCGCGAGGACTCCAGCTCTGGATGCTcagc                                                                                         >  1:90377/1‑251 (MQ=255)
                                                                                                                                                             gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGCCCCGGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                                              <  1:175312/246‑1 (MQ=255)
                                                                                                                                                              cccGCACACCAGAATCAATTTTTGATTCGTTGCCATTTCGCTGCCGCCCGCTTGACCAAGCTTGCGCCCCTGTGGGGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTcagca                                                                                        <  1:562911/250‑1 (MQ=255)
                                                                                                                                                                   acacCAGAATCAATTTTTGATTCGTTGCCAGCTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGccc                                                                                   >  1:144123/1‑250 (MQ=255)
                                                                                                                                                                          aaTCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTa                                                                             >  1:276903/1‑250 (MQ=255)
                                                                                                                                                                          aaTCAATTTTTGATTCGTTGCCAGTTCGCTGACCCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTaa                                                                            <  1:346331/251‑1 (MQ=255)
                                                                                                                                                                                     gATTCCTTGCCATTTCGCTGACGCCCGTTTAATCAGTCTTGGTCCCTGTGGATAACAGATAAATCACCTTTGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGTTCCGCCGCCGTAAGGACGATCGTCCTCGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGcc                                                                  <  1:574500/250‑1 (MQ=255)
                                                                                                                                                                                          gTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTAc                                                            >  1:98075/1‑251 (MQ=255)
                                                                                                                                                                                          gTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTAc                                                            >  2:278171/1‑251 (MQ=255)
                                                                                                                                                                                           ttGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTAcc                                                           >  1:495576/1‑251 (MQ=255)
                                                                                                                                                                                                aGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGt                                                       >  1:298084/1‑250 (MQ=255)
                                                                                                                                                                                                          cGCCCGCTTGATCAAGCTGGCGTCCTTGTGGTGACAGATAAACCCCTTTTGCGCCTGCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAAtt                                            <  1:80399/251‑1 (MQ=255)
                                                                                                                                                                                                           gCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTg                                           >  1:92115/1‑251 (MQ=255)
                                                                                                                                                                                                              ggtttGATCAAGCTTGTGTCCTTGGGGTGACAGAAAAATCCCCTTCGCCCCTCCCCCCGCCGCGGCTTTTGCTGCTTGAATGGCGTCCCGCAAGGGTTGCCCCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTGAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcgc                                        <  2:417590/248‑1 (MQ=255)
                                                                                                                                                                                                              cGCTTGATCAGCCTTGCTTCCCTGGGTTGACAGATAAATCCCCTTGGCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcgc                                        <  1:328029/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAAc                                   >  1:542693/1‑250 (MQ=255)
                                                                                                                                                                                                                    aTCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAAc                                   >  2:406094/1‑250 (MQ=255)
                                                                                                                                                                                                                               gTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                       <  1:106627/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGACCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACa                      <  1:37996/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cccTTTGGTGACAAATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACGGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAg                     <  1:458326/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                       >  2:369190/1‑249 (MQ=255)
                                                                                                                                                                                                                                 cccTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc                       >  2:574659/1‑249 (MQ=255)
                                                                                                                                                                                                                                 catgggggtgACAGAAAATCCCCCTTGGCGCCTCCTCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATGAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAg                     <  2:382900/245‑1 (MQ=255)
                                                                                                                                                                                                                                    ggtggtgACAGATAAACCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCt                  <  1:355065/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gtgACAGATAAATCCCCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAAtt              <  1:159527/251‑1 (MQ=255)
                                                                                                                                                                                                                                              agttAAATCCCTTTTGCGCCTACCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCTCCATTAACACCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCAGGCGTGAAGTCGCGAGGACTCCAGCTCGGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGcc         <  1:450640/248‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaaaccACCTTTGCGCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAAt      <  1:291302/246‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaaTCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATc     >  2:340389/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   aaTCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATcc    >  1:502703/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      cACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCcac  >  1:98272/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         
GCCGAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCAC  >  NC_000913/2744843‑2745337

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: