Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,323,258 A→C 6.5% intergenic (‑28/+78) phnE ← / ← phnD defective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily/phosphonate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,2580AC 6.5% 217.7 / 3.0 76intergenic (‑28/+78)phnE/phnDdefective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily/phosphonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base C (4/3);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCGCCCAGC  >  NC_000913/4322962‑4323550
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                    
aaGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATGTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCCCAGAAGCGAGAACCAGCTGAGCGTGGGTGGGGCGAGGGGGATGGTTTGCATGTGTGGATCCGGTCTGGAGGCCGGATCAGGTg                                                                                                                                                                                                                                                                                                  >  1:449779/1‑301 (MQ=255)
  ggggATGGAGAGACACACCGCCAGCGCGTTGCCCCAGCCGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGACCGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCt                                                                                                                                                                                                                                                                                                <  2:450343/301‑1 (MQ=255)
  ggggATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCCCTCGGGTGGGGCGATGGTGATGGTTTGACAGTTTGACTCCGGTCAGTAAGCCGGATCAgggcct                                                                                                                                                                                                                                                                                                >  1:589819/1‑297 (MQ=255)
          agagaACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGc                                                                                                                                                                                                                                                                                        >  2:123959/1‑301 (MQ=255)
             gaACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGCAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCCGCGCTTGGGCGGGGCGCTGGTGATGGTTGGCATGGTTGGCTCCGGTCGGTCGGCCGGAACAGGGACTTGCAccacctc                                                                                                                                                                                                                                                                                     >  1:67895/1‑296 (MQ=255)
             gaACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGCCGACCAACACGGCGAGTACCACGGCCCAGCCCAGAAGCGAGCACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGGTTGCATGTTTGGCTCCGGTCCGTAGGCCGGAAAAGGATCTTGCCCCCCacc                                                                                                                                                                                                                                                                                     >  1:36233/1‑299 (MQ=255)
                ccaccGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCgg                                                                                                                                                                                                                                                                                  >  2:8242/1‑301 (MQ=255)
                  accGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACCCGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCCTGGGTGGGGCAATGGTGAAGGTTTGCATGATTGGCCCCGGTCTATAGGCCGGAAACGGACCTTGCACCGCATCCGGcc                                                                                                                                                                                                                                                                                >  1:553239/1‑300 (MQ=255)
                        gtggggggtccccGGACGGGGTTTGGGAGGGTGACGGCCATTGCGGGGGGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCAGGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATGCGGTGCAAGCACCTTATCCGGcctacaga                                                                                                                                                                                                                                                                          <  1:312833/292‑8 (MQ=255)
                            cgttGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATGCGGTGCAAGCACCTTATCCGGcctacagaccgg                                                                                                                                                                                                                                                                      <  2:243023/298‑12 (MQ=255)
                                 gtcgcgtttccggtctGCAGCGTGACGGTCAGTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACAt                                                                                                                                                                                                                                                                 <  1:456966/286‑1 (MQ=255)
                                       gatggtcttccggttttGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCgg                                                                                                                                                                                                                                                           <  1:279495/285‑1 (MQ=255)
                                        cGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGa                                                                                                                                                                                                                                                          >  2:600532/1‑301 (MQ=255)
                                           cgttttGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGc                                                                                                                                                                                                                                                       <  2:36284/298‑1 (MQ=255)
                                           cGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGc                                                                                                                                                                                                                                                       >  2:120837/1‑301 (MQ=255)
                                              tttGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGAc                                                                                                                                                                                                                                                    <  2:67985/301‑1 (MQ=255)
                                                 gCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACAtagatcggaagagcaca                                                                                                                                                                                                                                                 >  2:457530/1‑285 (MQ=255)
                                                      gTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGc                                                                                                                                                                                                                                            <  2:590580/301‑1 (MQ=255)
                                                       tGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGagatcggaagagcaca                                                                                                                                                                                                                                           >  2:279798/1‑286 (MQ=255)
                                                          cGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGAAGCGAAGCCGGGGCGTCt                                                                                                                                                                                                                                        >  2:266756/1‑301 (MQ=255)
                                                           ggacttttCGGTGGGGTAATCTTGCCACGGGCTGAAATGGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCtt                                                                                                                                                                                                                                       <  1:599771/296‑1 (MQ=255)
                                                               gtgtCGGGGGGGTATTCCGGGCCCGGGCTGAAGGCGGGAGGGAAGAAGTGGGCGGGGGACGTCGCCATGTTGCCGCCGGCTTTGATCAGCGTAAGCGGGGCCATGTCCGGGCCCGGCCACGGGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATcc                                                                                                                                                                                                                                   <  1:273797/298‑1 (MQ=255)
                                                               aTTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATcc                                                                                                                                                                                                                                   >  2:76905/1‑301 (MQ=255)
                                                                tttCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCg                                                                                                                                                                                                                                  <  2:610591/301‑1 (MQ=255)
                                                                 ggtgggggggTAATCTTGCCATTGGTTGAATTCGGGCGGGAAGAAGTGGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGGAAGCGGGGCCATGTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATAAGACGCCCCGGCGTCGCATCCGa                                                                                                                                                                                                                                 <  1:610100/293‑1 (MQ=255)
                                                                   cGGTGGGGTACTCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGGTGGTGATGGTTTGCCTGTTTGGCTCCGGTCTGTCGGCCGGATACGGGGCTTGCCCCACCACCGCCCTCCACGCCTGCCCATTGTATGTCGCGACGCCcggcagtcctcaccgcac                                                                                                                                                                                                                               >  1:105433/1‑286 (MQ=255)
                                                                         ggTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGaa                                                                                                                                                                                                                         >  2:70634/1‑300 (MQ=255)
                                                                            aaTCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCCACAACGCCTGCACATTGTCGGATGCGACGCCCGGGCGTCCTATCCGACCTACGATTTa                                                                                                                                                                                                                      >  1:302579/1‑301 (MQ=255)
                                                                              ctgggcgagtGGCTGAATGCGGGGGGGAAGAAGGCGGGGGCGAACTTGGCCAGGTTGCCGCCGTCTTTGAGCAGCGTAAGCGGGGCCATTTCCGCGCCCGGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCGCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACt                                                                                                                                                                                                                    <  1:76815/292‑1 (MQ=255)
                                                                              cgctcttccgggttgacagcgggcgggAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATGTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACt                                                                                                                                                                                                                    <  1:401547/282‑1 (MQ=255)
                                                                                ttGCCAGTGGCTGAAATCGGGCGTGAAGAAGTCGGCGTCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGc                                                                                                                                                                                                                  <  2:434983/300‑1 (MQ=255)
                                                                                       gggCTGAGATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCGGCGTAAGCGGGGCCATGTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCtt                                                                                                                                                                                                           <  1:526251/300‑1 (MQ=255)
                                                                                        ggCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCCGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTCGGCCGGATACGGTGCTTGCACCGCATCCGGCATCAACGCCCGAACATTCTCGGATGCGACGCCGGGGCGGCATATCCGACCCCACATTCACTGCACCCCtct                                                                                                                                                                                                          >  1:377616/1‑299 (MQ=255)
                                                                                        ggCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGCTGGTGATGGGTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCAACGCATCCGGCATCAACGCCTGCACATTGTCCGACGCGACGCCGGGGCGCCCTCTACGACCCACGATTTACTGCACCCCttt                                                                                                                                                                                                          >  1:64320/1‑301 (MQ=255)
                                                                                        ggCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGCTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACTCCTGCACCTTGTCGGACGCGACACCCGGGCCTCCTATCCGACCCACGCTTTACTACACCGCctt                                                                                                                                                                                                          >  1:76875/1‑301 (MQ=255)
                                                                                          ttGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTAc                                                                                                                                                                                                        <  2:302983/300‑1 (MQ=255)
                                                                                          cTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTAACCGACCTACGATTTACTGCACCGCTTTaa                                                                                                                                                                                                        >  2:179/1‑300 (MQ=255)
                                                                                           ggAAGTCGGGCGGGAGGAGGTCGGCGGCGAACGTCTCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACt                                                                                                                                                                                                       <  1:266488/300‑1 (MQ=255)
                                                                                            gAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACtagatcggaagagca                                                                                                                                                                                                      >  2:402061/1‑287 (MQ=255)
                                                                                             aagtCGGGGGGGATGAATGCGGCGGCGAACGTGGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCa                                                                                                                                                                                                     <  1:70544/298‑1 (MQ=255)
                                                                                             aaaTCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCa                                                                                                                                                                                                     >  2:365287/1‑301 (MQ=255)
                                                                                                              gTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCta                                                                                                                                                                                    >  2:441923/1‑300 (MQ=255)
                                                                                                                       gAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgttt                                                                                                                                                                           >  2:69098/1‑300 (MQ=255)
                                                                                                                             cGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                                                                                                                     <  2:105557/301‑1 (MQ=255)
                                                                                                                                      gccgccGTCTTTGATCATCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAgg                                                                                                                                                            <  2:64388/301‑1 (MQ=255)
                                                                                                                                                 tGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTg                                                                                                                                                 >  2:567808/1‑301 (MQ=255)
                                                                                                                                                 ctttcccgacgtgcggcgcttccgatctgcCCTGCCACGAGACGACCAACGCGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTg                                                                                                                                                 <  1:640896/273‑1 (MQ=255)
                                                                                                                                                       gCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCCGCTGCGCCTg                                                                                                                                           >  2:52924/1‑301 (MQ=255)
                                                                                                                                                          tAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCgcagctggatgacctggaccgcctgaac                                                                                                                                        <  2:378112/301‑28 (MQ=255)
                                                                                                                                                                  gCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCCTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCCGGCGGTCCAGGTCATCCAGCATCGCCTGAATCCCggtgg                                                                                                                                >  1:318341/1‑301 (MQ=255)
                                                                                                                                                                       ttCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGttttt                                                                                                                           <  2:553912/301‑1 (MQ=255)
                                                                                                                                                                         ccGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGCTAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCCACGAACCCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGATCTTTGc                                                                                                                         >  1:237467/1‑301 (MQ=255)
                                                                                                                                                                             gcCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTgagatcggaagagcacacgtctgaactcc                                                                                                                     >  2:641773/1‑273 (MQ=255)
                                                                                                                                                                              ccgggcagcGAGACGGCCAACAGGGCGAGTAAAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCt                                                                                                                    <  1:276422/293‑1 (MQ=255)
                                                                                                                                                                                   ccACGGGACGGACAGCACGGGTGGTACGAGGGCCCAGCTCGGAAGCGGGAAGCGGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTcc                                                                                                               <  1:123835/301‑1 (MQ=255)
                                                                                                                                                                                              ccAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCaa                                                                                                    >  2:632386/1‑300 (MQ=255)
                                                                                                                                                                                                       cGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTAtt                                                                                           >  2:135579/1‑301 (MQ=255)
                                                                                                                                                                                                                        ctttttccctacactgggttcttcggttctgggGCGATGGTGCTGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCa                                                                          <  1:608367/272‑1 (MQ=255)
                                                                                                                                                                                                                              aGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCCTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCtt                                                                    <  2:318745/301‑1 (MQ=255)
                                                                                                                                                                                                                                 gggagcgcgcgggggttGGGGGGGGGGGTGGTGGTGGTTTGCATGTTTGGGTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCAGCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGGCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGGTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTaa                                                                 <  1:8242/286‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGcgcg                                                         >  2:466639/1‑301 (MQ=255)
                                                                                                                                                                                                                                                ttGGTTTGGGTGTTGGTGTTGGTTTGCCTGTTTGGCTCCGGTCTGTAGGCCGGATAAGTTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGc                                                  <  1:135454/301‑1 (MQ=255)
                                                                                                                                                                                                                                                   cgtgcttttcgtgtgttggtTTTCTTTTTTTGCTCCGGTCTGTCGGCCTGTGTCTGTGCTTGCACCGCTTCCGGCATCAACGCCTGCACATTGTCTGATGCGACGCCGGGGCGTCTTATCCGACCTACGTGTTACTGCCCCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGaa                                               <  1:114467/289‑1 (MQ=255)
                                                                                                                                                                                                                                                    gtgggGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCCGGCGGTCCAGGTCATCCAGCTGCGCCTGCATCTCGGTGGTTTTCGCCCGCTTGCCCTTCCCATTCAGTCCTTTATTACTTTTCACTCCCTGCATCTCTTTAAAcctcaccacctggctacc                                              >  1:633109/1‑284 (MQ=255)
                                                                                                                                                                                                                                                      ggggCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCaagatcggaagagcacaagtctgaactccaa                                            >  2:609225/1‑271 (MQ=255)
                                                                                                                                                                                                                                                        ggCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGc                                          >  2:619113/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 gttggtTTGCATTTTGGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGAGGCGACGCCGGGGCGTCTGATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGc                                 <  1:631509/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAg                               >  2:14963/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                   tggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAACTGGCGaaagaacggaagagcaca                               >  2:114591/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                    ggtttgggttgttGGCTCCTGTCTGTATGCCTGATATGGTTCTTGCACCGCATCCGGCATCAACGTCTGCACATTGTCGGATGCGACGCCGGGGCGTCTGATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTGATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAgg                              <  1:120713/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gtgtttttGGCTCCGGTCTGTGTGTCGGTTATGTTGCTTGCACCGCATCCGGCATCAACGCCTTCCCTTTTTCGGATGCGACGCCGGGGCGTCTTATCCGACTTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCt                         <  1:14963/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGa              >  2:61719/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     gtggggtGTAGGCTGGATATGGTGCTGTCACTGCAGCCGGCATCAAGGCTTGCACGTGGGCGGGTGCGACGCCGGGGCGTCTTATCCGACCTACGATGTACTGCACCGCTTTACTCACCTAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGaa             <  1:179/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             gaGGCCGGTTAAGTTGCTTGGCCGGCATCCGGCATCAACGGCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCgccc     <  1:69008/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              aGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCgccca    >  2:413499/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                gCCGGATGCGGTGCAAGCACCTTATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCgcccagc  >  2:499470/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                    
AAGGGGATGGAGAGAACCACCGCCAGCGCGGTGCCCCAGACGGCGATTTGCAGCGTGACGGCCATTTCGGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCGCCCAGC  >  NC_000913/4322962‑4323550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: