Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,585,673 C→T W363* (TGG→TGA hsdR ← type I restriction enzyme EcoKI endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,585,6730CT98.7% 246.6 / NA 78W363* (TGG→TGAhsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (37/39);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATT  >  NC_000913/4585426‑4585914
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
gCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGGTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCTCGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAg                                                                                                                                                                                                                                                <  1:307294/251‑1 (MQ=255)
  aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCt                                                                                                                                                                                                                                              >  1:413309/1‑251 (MQ=255)
  aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCt                                                                                                                                                                                                                                              <  1:304903/251‑1 (MQ=255)
     cccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGt                                                                                                                                                                                                                                           <  1:222069/251‑1 (MQ=255)
      ccAGCTCGCTCATGCCAGGGTTATCGGCAACCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAATGCTTTCGACATATTGCGGGTATCACGAACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGtt                                                                                                                                                                                                                                          >  2:458515/1‑251 (MQ=255)
                       gggTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGc                                                                                                                                                                                                                         <  1:388406/251‑1 (MQ=255)
                        ggatATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATAACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCACGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCt                                                                                                                                                                                                                        >  2:418554/4‑251 (MQ=255)
                           tagcGGCAAACCACTGATTCTGGTTTTGCGGGGCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTGGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCAGACTCTGGCACTGCTTc                                                                                                                                                                                                                      <  2:31913/247‑1 (MQ=255)
                            aTCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACACGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCAt                                                                                                                                                                                                                    >  2:39414/1‑251 (MQ=255)
                                  aaaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACtt                                                                                                                                                                                                              <  1:411749/251‑1 (MQ=255)
                                     ccACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTTCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAt                                                                                                                                                                                                           <  2:431966/251‑1 (MQ=255)
                                     ccACTGAGTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCGTCCGGGCGGGGCCACTCGGGTAAGGCTTTCGACATAGTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTGGGTCTTCATTGTTGCGCGGTATTCTCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTTTTTGCCGCTGTTGTCCTGTCAGCTGGTTTCAGACTCTGGCACTGCGTCATGCACTTCag                                                                                                                                                                                                           <  2:565763/251‑2 (MQ=255)
                                       aCTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATcc                                                                                                                                                                                                         >  1:348408/1‑251 (MQ=255)
                                       aCTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCGTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCGGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATcc                                                                                                                                                                                                         <  1:562404/251‑1 (MQ=255)
                                             tCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGa                                                                                                                                                                                                    <  2:510122/250‑1 (MQ=255)
                                               tGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGCATTCGCGCCCGTTGGTCGAGTAGCAGACGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagt                                                                                                                                                                                                  >  1:252426/1‑250 (MQ=255)
                                                gTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGTGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagaat                                                                                                                                                                                                >  2:4783/1‑248 (MQ=255)
                                                ctttttGCGGTTCGCTGCCGAGCATTTCCAGAAGCTCTTCAGGGCGGTGCCACTCGGGAACGTCTTTCGACATATCGCGGGTATCATGCACGCCGCGATACCAGATGACGCTTTTGGTCTTCATTGTTGCGCGGTATTCACGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGGGTCCTGCCAACTGGTTTCATACGATGGAACTACTTCATGCACTTCATCCGGTGagt                                                                                                                                                                                                  >  2:359486/2‑249 (MQ=255)
                                                  ttttGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagtagt                                                                                                                                                                                               <  1:577830/250‑1 (MQ=255)
                                                  ttttGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGGTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTg                                                                                                                                                                                              >  2:1891/1‑251 (MQ=255)
                                                      gccgTTCGCTGCCGAGCATTTCCCGCAGCTTTTTCGGGCGGTGTGACGCGGGTCACGATTTCGACTTATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCTTCCGTTGGGCGAGTAGCAGCAGGGGAATTTAAACCGTTGTTTGCCGCTAGTGTCCTGGCAGCGGGTTTCATACTCTGTCACTGCTTCATCCACTTCATCCGTTGAGTATTGCTc                                                                                                                                                                                           <  2:317880/247‑1 (MQ=255)
                                                           tCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaagca                                                                                                                                                                                      >  2:167110/1‑250 (MQ=255)
                                                            cGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaagcaag                                                                                                                                                                                    <  1:271527/251‑1 (MQ=255)
                                                                gCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCTCTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGGCTTCATTGTTGCGCGGTATTCGCTCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGGTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATCCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGtt                                                                                                                                                                                 <  2:411933/250‑1 (MQ=255)
                                                                    ggCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATCCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCg                                                                                                                                                                            <  2:95382/250‑1 (MQ=255)
                                                                    aGCATTTACAGCCGCTCTTCCGGCCGGCGCCACCCGGCTAAGGCTCTCGACATATTGCGGGTATCACGCACGTAGCGATACCAGATGCCGCTTTTGGCCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGCAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCACGCAAGGTTTCCCg                                                                                                                                                                            >  2:342684/1‑251 (MQ=255)
                                                                          tCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGa                                                                                                                                                                       >  1:476329/1‑250 (MQ=255)
                                                                           ccAGCAGCTCGTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGaa                                                                                                                                                                      <  1:284712/250‑1 (MQ=255)
                                                                                agcTCTTCCGGGCGGTGCCACTCGGGTACGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGCTACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAtt                                                                                                                                                                >  2:413797/1‑251 (MQ=255)
                                                                                      cccGGGCGGTGCCACGCGAGTAAGGCTTTCGACATATTGCGGGTATCACGCACTTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGGCGGTTGGTCGAGTAGCATAAGGGGATTTTAAACCGTTGGTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATGCTCTGGCAGTGCTTCATGCACTGCATCCTGTGAGAAGTGCGCAAGCAAGGTTGCCTGCAGGAAGTCATCATCGaa                                                                                                                                                          <  2:67796/250‑1 (MQ=255)
                                                                                       ccGGGCGGTGCCACTCGGGTAAGGCTGTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCTTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTACATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGaa                                                                                                                                                          <  2:348408/250‑1 (MQ=255)
                                                                                                gCCACTCGGGTAAGACTTTCGACATATTGCGGGTATCACGCACGTCGCGCTACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGAATTCGCGCCCGTTGGTCGAGTCGCAGAAGGGGATTTTAAACAGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCTTGCACTTCATCCGGTGGGTAGTGCTCAAGCAAGGTTTCCCGCCGGAAGCCATTATCGACACATTTACTg                                                                                                                                                >  2:562201/1‑251 (MQ=255)
                                                                                                  cACTCGTGTAAGGCTTTCGACAGATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAGTTACTGTa                                                                                                                                              <  2:472287/251‑1 (MQ=255)
                                                                                                   aCTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATAACAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGt                                                                                                                                               >  2:296132/1‑249 (MQ=255)
                                                                                                         gTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCATGAAGCCATTATCGAAACATTTACTGTAGCGATAc                                                                                                                                       <  2:413546/251‑1 (MQ=255)
                                                                                                         gTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTATCGATAc                                                                                                                                       <  2:513706/251‑1 (MQ=255)
                                                                                                          tAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACg                                                                                                                                      >  1:555998/1‑251 (MQ=255)
                                                                                                             ggCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACt                                                                                                                                   <  1:20256/251‑1 (MQ=255)
                                                                                                                   cGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCTGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCTTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGa                                                                                                                             <  1:585374/251‑1 (MQ=255)
                                                                                                                         atTGCGGGTATCACGCACGTCGCGATTCCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGGCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCGCTGCTTCATGCACTTCATCCGGTGAGTAGTGTTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATGTTCTTTAGCGATACGACTCATTGAGCCTg                                                                                                                        <  2:297155/250‑1 (MQ=255)
                                                                                                                         atTGCGGGTATCACGCACGTCGCGATGCCAGATGGCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTGCTGTAGCGATACGACTCATTGATCCTgg                                                                                                                       <  2:547674/251‑1 (MQ=255)
                                                                                                                         atTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTgg                                                                                                                       >  2:378551/1‑251 (MQ=255)
                                                                                                                           tGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCg                                                                                                                     >  1:249427/1‑251 (MQ=255)
                                                                                                                              gggTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCggg                                                                                                                   >  2:456740/1‑250 (MQ=255)
                                                                                                                               ggTATCACGAACGTCGCGATACCAGATGCCGCTTTTGGTCTTAATTGATGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGac                                                                                                                 >  2:588814/1‑250 (MQ=255)
                                                                                                                                     acgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCg                                                                                                            >  1:180001/1‑250 (MQ=255)
                                                                                                                                      cgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGat                                                                                                          <  1:447904/251‑1 (MQ=255)
                                                                                                                                       gcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGat                                                                                                          >  2:40910/1‑250 (MQ=255)
                                                                                                                                          cgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTg                                                                                                      >  1:215885/1‑251 (MQ=255)
                                                                                                                                           gtcgagatgacAGATGCCGCTTTTGGTCTTCATTCTTGCTCAGTATTCGCGCCAATGGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGt                                                                                                     <  1:149255/241‑1 (MQ=255)
                                                                                                                                               cgATACCGGATGCCGCTTTTGGTCGTCATGCTTGCGCGGTATTCGCGACCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGTTGTCCTGTCAGCTGGTTGCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCTATATTGTAACg                                                                                                 <  2:166100/251‑1 (MQ=255)
                                                                                                                                                   aCCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTc                                                                                             <  1:34273/251‑1 (MQ=255)
                                                                                                                                                            cGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTacc                                                                                    <  1:133085/251‑1 (MQ=255)
                                                                                                                                                               ttttGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccacc                                                                                 >  2:201627/1‑251 (MQ=255)
                                                                                                                                                                tttGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGCTATTGTTACGTTTCGCATCTACCACcg                                                                                >  1:368349/1‑251 (MQ=255)
                                                                                                                                                                    gTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGAt                                                                            <  1:373249/251‑1 (MQ=255)
                                                                                                                                                                     tCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGTTATTGTTACGTTTCGCCTCTACCACCGCGATg                                                                           >  2:147065/1‑251 (MQ=255)
                                                                                                                                                                     tCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCg                                                                                            <  1:285565/234‑1 (MQ=255)
                                                                                                                                                                     tCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCg                                                                                            >  2:285565/1‑234 (MQ=255)
                                                                                                                                                                            tgttgCGCGGTATTCGCGCCCGTTGCTCGAGCACCAGAAGGGGATTTTAAACCTTTGTTTCACGCTGTTGTCCTGTCCGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGCGTATTCCTCAAGCAAGGTTTCCCGCAGGAAGCCATTCTCGAACCATTTACTGTAGCGATACGCCTCATTGAGCCTGGCGGGCACGTCGATCTTTTTCCGTTTCGCCTCTACCCCCGCGATGGGTTTgc                                                                    >  2:588627/1‑250 (MQ=255)
                                                                                                                                                                             gttgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCCCCGCGATGGGTTTg                                                                     >  1:552174/1‑249 (MQ=255)
                                                                                                                                                                                                    ggTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATCCGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGACGACCAACAGCACATAAtacg                                            >  1:512791/1‑248 (MQ=255)
                                                                                                                                                                                                        gAGTAGCAGATGGGCCTTTTAAAACGTTGTTTGCCGCTGGTGTCCCCTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGATCAAGCAAGGTTTCCCGCTGGAAGCCATTATCGCAACATTTACTGTAGAGATACGACTCATTGAGCCTGGCGGGAACGTCGCTATTGTTACGTTTCGCCTCTACCACAGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCaa                                         >  2:89538/1‑250 (MQ=255)
                                                                                                                                                                                                          gTAGCAGAAGGGGATTTTAATCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGTTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTGCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACGCATTGAGCCTGGCGGGAACGTCGATATTGTTACGGTTCGCCTCTACCTCCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGc                                      <  2:555998/251‑1 (MQ=255)
                                                                                                                                                                                                               aGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATACTCCGCAAAGCCCTg                                  >  1:210971/1‑250 (MQ=255)
                                                                                                                                                                                                                 aaGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGAtt                               >  1:43138/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGAtt                               <  1:40910/251‑1 (MQ=255)
                                                                                                                                                                                                                  aGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCCCATAATCCGCAAAGCCCTggtt                               >  2:143334/1‑247 (MQ=255)
                                                                                                                                                                                                                     ggATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCGCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTGGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTAccc                           <  1:90883/251‑1 (MQ=255)
                                                                                                                                                                                                                       aTTTTAAACCGTTGTGTGCCGCTTGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACGCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGt                         <  1:229806/251‑1 (MQ=255)
                                                                                                                                                                                                                        ttttAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGCTACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGGTGGGTTTGAGGCCGACAAAAAGAACATAATCATCAAAGCCCTGATTCCCCGtt                        >  2:48088/1‑251 (MQ=255)
                                                                                                                                                                                                                          ttAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGttt                       <  2:215885/250‑1 (MQ=255)
                                                                                                                                                                                                                          ttAAACCGTTGTTTGCCGCTGGGGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGGAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGttt                       <  2:6868/250‑1 (MQ=255)
                                                                                                                                                                                                                            aaaCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAATCCCTGATTACCCGTTTCa                     <  1:163081/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCTGTTTCATCTTTTCCg            >  1:102464/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCg        <  1:167110/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ccGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCgg       >  2:550719/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ggTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCAtt  >  2:115664/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
GCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATT  >  NC_000913/4585426‑4585914

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