Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,678 G→A P18S (CCA→TCA)  dnaA ← chromosomal replication initiator protein DnaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,6780GA96.3% 336.1 / ‑0.9 112P18S (CCA→TCA) dnaAchromosomal replication initiator protein DnaA
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (43/62);  total (47/65)
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

Only 100 of 112 total aligned reads displayed.

GGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTC  >  NC_000913/3883437‑3883923
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
ggggCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATc                                                                                                                                                                                                                                                <  2:462826/251‑1 (MQ=255)
   gCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTg                                                                                                                                                                                                                                             <  2:326610/251‑1 (MQ=255)
    ccGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGc                                                                                                                                                                                                                                            >  1:651989/1‑251 (MQ=255)
        gACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATc                                                                                                                                                                                                                                        >  1:734447/1‑251 (MQ=255)
              gCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCa                                                                                                                                                                                                                                   <  1:342318/250‑1 (MQ=255)
                 cGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGac                                                                                                                                                                                                                               <  1:77170/251‑1 (MQ=255)
                  gTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAACTCAGGGGCATCAGCTCCACAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGCTAACTCCTCCTGCAATCGGGCAAGac                                                                                                                                                                                                                               >  1:89677/1‑250 (MQ=255)
                        gCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACActgctg                                                                                                                                                                                                                        >  1:142328/1‑251 (MQ=255)
                        gCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACACAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACActgct                                                                                                                                                                                                                         <  1:247987/250‑1 (MQ=255)
                              tgtgGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACTAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCaaa                                                                                                                                                                                                                   <  1:540572/250‑1 (MQ=255)
                                    gTTTGAGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGCTATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGAGCGTACAGGGCACGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTTGCTGCTAACTCATCCTGCACTCGGGCAAGACACTGCTGCCAACGCGaacg                                                                                                                                                                                                            >  2:54270/1‑249 (MQ=255)
                                            cACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGTCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCTAAAGTGACAcg                                                                                                                                                                                                     <  2:301470/250‑1 (MQ=255)
                                            cACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACGGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACAcgg                                                                                                                                                                                                    <  2:294847/251‑1 (MQ=255)
                                                ggTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGCATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACAcggcgg                                                                                                                                                                                                 >  1:682230/1‑250 (MQ=255)
                                                         ccGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATTCTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCg                                                                                                                                                                                       <  1:680688/251‑1 (MQ=255)
                                                            aCTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAAc                                                                                                                                                                                    <  1:555692/251‑1 (MQ=255)
                                                             cTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAAc                                                                                                                                                                                    >  2:436109/1‑250 (MQ=255)
                                                                cAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTGAGCAGTGCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACaaa                                                                                                                                                                                 <  2:262908/250‑1 (MQ=255)
                                                                   aaCGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTc                                                                                                                                                                             <  1:777864/251‑1 (MQ=255)
                                                                       cAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGAt                                                                                                                                                                          >  1:194688/1‑250 (MQ=255)
                                                                        aGCTGTTGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAAGCGATGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATaa                                                                                                                                                                        <  1:208719/251‑1 (MQ=255)
                                                                          gtgtgGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATg                                                                                                                                                                      <  2:254899/250‑1 (MQ=255)
                                                                            gtgGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTg                                                                                                                                                                                                          >  2:417079/1‑213 (MQ=255)
                                                                            gtgGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTg                                                                                                                                                                                                          <  1:417079/213‑1 (MQ=255)
                                                                              ggggCATCCGCTCCGCAGAAACTGGTTAGCAGTACATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCCCTCAGTTACGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGGTAACTCATCCTGCAATCGGGACAGACCCTGCTGACAAAGCGAACGTGACACGGCAGCCTCCACTCTAAAAAAAgacaaaaaagccg                                                                                                                                                                   >  2:555958/1‑238 (MQ=255)
                                                                                 gCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAgg                                                                                                                                                               <  1:597937/251‑1 (MQ=255)
                                                                                 gCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCGCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAgg                                                                                                                                                               <  2:236244/251‑1 (MQ=255)
                                                                                     ccgctccgcAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGa                                                                                                                                                           <  1:568647/251‑1 (MQ=255)
                                                                                          ccgcAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACAt                                                                                                                                                      <  1:585374/251‑1 (MQ=255)
                                                                                           cgcAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACAtt                                                                                                                                                     <  1:127073/251‑1 (MQ=255)
                                                                                           cgcAGAAACTGGGTAGCAGACCATTGATATTATTAAGGTACTTGTCCCGTACCCAATAGAGGACAAAACGGTTTGGCGCGTACAGGGACAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGATGATAACTCATCCTGCAATCGGGCAAGAAACTGCTGCCAAAGCGAAAGTGACACGGCGGAATCAACTGGAAAAAAAGTCGATCAAGACTAAAGCTGAAACAtt                                                                                                                                                     >  2:75390/1‑251 (MQ=255)
                                                                                             cAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCa                                                                                                                                                   <  1:784141/251‑1 (MQ=255)
                                                                                                aaaCTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATg                                                                                                                                                 >  1:104422/1‑250 (MQ=255)
                                                                                                 aaCTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATAACTCAGACTGAAACATTCAttat                                                                                                                                               >  2:57637/1‑248 (MQ=255)
                                                                                                     ggTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttgt                                                                                                                                            <  1:636195/250‑1 (MQ=255)
                                                                                                      gTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATAATGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGAttgttg                                                                                                                                          >  1:380875/1‑251 (MQ=255)
                                                                                                        tAGCAGTCCATTGATATTATTAAGGTTCTTTTCCCGTACCCAAACGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTAGCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACTGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGa                                                                                                                                         <  2:774627/250‑1 (MQ=255)
                                                                                                        tAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTc                                                                                                                                                                                        <  2:729258/203‑1 (MQ=255)
                                                                                                        tAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTc                                                                                                                                                                                        >  1:729258/1‑203 (MQ=255)
                                                                                                         aGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGGTTGTTGacg                                                                                                                                       >  2:353732/1‑251 (MQ=255)
                                                                                                           cAGTCCAATGGTATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGGGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACTGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACAGTCATGATTGTTGacgta                                                                                                                                     <  2:353326/251‑1 (MQ=255)
                                                                                                            agaacATTGATATTATTAAGGTCCTTGTCCCGTACCCAATCGCGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCAAAAGCGAAAGGGACACGGCGGACTCCAATCGAACAAAAGACGATAATGACAAAGGATGAAACATTCATGATTGTTGacgtac                                                                                                                                    >  2:528587/5‑251 (MQ=255)
                                                                                                             gTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGacgtacg                                                                                                                                   >  1:380173/1‑251 (MQ=255)
                                                                                                                cATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATAGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGAGACACGGCGGACTCCACTCGAACAAAAGTCGCTAATGACTAAGGCTGAAACATTCATGATTGTTGACGTCCGTCg                                                                                                                                >  1:111685/1‑251 (MQ=255)
                                                                                                                       attattAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGAc                                                                                                                         >  1:661106/1‑251 (MQ=255)
                                                                                                                       attattAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTAGCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGa                                                                                                                          <  2:194688/250‑1 (MQ=255)
                                                                                                                         tattaAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGAccc                                                                                                                       >  1:864673/1‑251 (MQ=255)
                                                                                                                              aGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCtt                                                                                                                  >  2:312372/1‑251 (MQ=255)
                                                                                                                              aGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCtt                                                                                                                  <  1:436109/251‑1 (MQ=255)
                                                                                                                                gTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTgg                                                                                                                >  1:390526/1‑250 (MQ=255)
                                                                                                                                   cTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATa                                                                                                                                                                         <  1:360907/191‑1 (MQ=255)
                                                                                                                                   cTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATa                                                                                                                                                                         >  2:360907/1‑191 (MQ=255)
                                                                                                                                       ccccGTACCCAATCGAGGTCAAAACGTTTTGGCGCGTCTAGGGCCATGTTGTTAGCGCCCAGTTCCGCCTGCAATGGGCGTATCCACAAACTGAATTC‑TGGGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAAAAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGGCGAAAAGACCCTGCTTGAGGGTGAc                                                                                                         <  2:402683/250‑1 (MQ=255)
                                                                                                                                            tACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGcccg                                                                                                    >  2:670013/1‑248 (MQ=255)
                                                                                                                                             aCCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGCATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGCAACATTCATGATTGTTGACGTACGTCGaaaaagatcggcagagcgccgtgtagggg                                                                                                   >  2:324523/1‑226 (MQ=255)
                                                                                                                                             aCCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGaaaa                                                                                                                            <  1:324523/226‑1 (MQ=255)
                                                                                                                                               ccAATCGAGGACAAAACGGGTTGGCGCGTACAGGGCCAGCGTGATATCACGCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGATGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGGGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGAACGACGAACAGACCCTGATTGAGGCTGACgaac‑gaac                                                                                                >  1:771969/1‑251 (MQ=255)
                                                                                                                                                     gAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTACCGTACGTCGAAAAGACCCTGCTTGAGGGGGACGCCC‑GAACCGcggt                                                                                          >  2:668609/1‑248 (MQ=255)
                                                                                                                                                       ggACAAAACGGTTTGGCGCGTACAGGGTCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCt                                                                                        <  1:179449/251‑1 (MQ=255)
                                                                                                                                                          aaaaaCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCg                                                                                     <  1:449270/250‑1 (MQ=255)
                                                                                                                                                           aaaaCGGTTTGGCGCGTACAGGGCCAGCGGGTTATCGATCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAACGCGAAAGTGACACGGCGGACTCCACTCGAACCAAAGACGATAATGAATAAGGCTGAAACATTCATGATTGTTGACGAACGTCGAAAAGACCCTGCATGAGGGTGAAGAAA‑GAACCGcaatcggggg                                                                                    >  2:629411/1‑242 (MQ=255)
                                                                                                                                                            aaaCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGt                                                                                   >  2:787828/1‑251 (MQ=255)
                                                                                                                                                              aaggTTTGGCGCGGAAAGGGCAAGCGTGTTATCGCTCAGTTCCGCCTGCAAAGGGCGTATCCACATACTGAATTC‑TGTGGCTGATACCTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTACGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGATTGAGGGTGACACAC‑GAACCGCTGTCAGCGGtt                                                                                  >  2:513871/3‑250 (MQ=255)
                                                                                                                                                              aCGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTta                                                                                 <  1:664979/251‑1 (MQ=255)
                                                                                                                                                                    tGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGc                                                                            <  1:851636/250‑1 (MQ=255)
                                                                                                                                                                    tGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGc                                                                            <  1:598815/250‑1 (MQ=255)
                                                                                                                                                                    tGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCAACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCACAGCGCCAGTGACACGGCGTACTCCACTCGAACACAAGTCGCTAATGCCTAAGGCTGAAACATTCATGATTGTTGACGTACGTCACAAAGACCCAGCTTGAGGGTGACGCCC‑GAACCGCTGTCTGCGGTTATATGcc                                                                           >  2:320097/1‑251 (MQ=255)
                                                                                                                                                                     ggCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGccc                                                                          <  1:60825/251‑1 (MQ=255)
                                                                                                                                                                        gcgTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGAt                                                                       <  2:199226/251‑1 (MQ=255)
                                                                                                                                                                          gTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCa                                                                     <  1:645371/251‑1 (MQ=255)
                                                                                                                                                                           tACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGCACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GCACCGCTGTCTGCGGTTATATGCCCGATCa                                                                     >  2:803280/1‑250 (MQ=255)
                                                                                                                                                                            ccAGGGCCAGTGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTTATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCa‑a                                                                   <  2:803720/249‑1 (MQ=255)
                                                                                                                                                                             cAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGa                                                                 <  2:864673/251‑1 (MQ=255)
                                                                                                                                                                                            aTCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATc                                                  <  1:730742/251‑1 (MQ=255)
                                                                                                                                                                                            aTCGCTCAGTTCCGCCTGCAATGGGCGTAGCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCGTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGAt                                                   <  1:142379/250‑1 (MQ=255)
                                                                                                                                                                                                cTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAA‑GAACCGCTGTCTGCGGGTATATTCCCGATCA‑AGATCCTGCAACACGATCggg                                               >  2:850800/1‑250 (MQ=255)
                                                                                                                                                                                                  cAGTTCCGCCTGCAATGGGCGTTTCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGAcc                                            <  2:529592/251‑1 (MQ=255)
                                                                                                                                                                                                   aGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACcg                                           <  1:850800/251‑1 (MQ=255)
                                                                                                                                                                                                      tCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACcgcg                                         <  1:788266/250‑1 (MQ=255)
                                                                                                                                                                                                      tCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACcgcg                                         <  1:320097/250‑1 (MQ=255)
                                                                                                                                                                                                        cGCCTGCAATGGGCGTATCCACATCCTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGAt                                      <  1:596952/251‑1 (MQ=255)
                                                                                                                                                                                                        cGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTtata                                                                               >  2:591466/1‑211 (MQ=255)
                                                                                                                                                                                                        cGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTtata                                                                               <  1:591466/211‑1 (MQ=255)
                                                                                                                                                                                                          ccTGCAATGGGCGTATCCACAGACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCa                                    <  1:811645/251‑1 (MQ=255)
                                                                                                                                                                                                           cggCAATGTGCGTATCCCTATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATTACTAAGGCTGAAACATTCATGATTGTTGACGTACGGCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCAt                                   <  2:310646/249‑1 (MQ=255)
                                                                                                                                                                                                             gCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATa                                  <  1:312372/250‑1 (MQ=255)
                                                                                                                                                                                                               aaTGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGcc                               >  1:150743/1‑251 (MQ=255)
                                                                                                                                                                                                               aaTGGGCGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATAA‑AGATCCTGCAAAACGATCGGGACCCCGGAtactagcc                               >  2:648860/1‑244 (MQ=255)
                                                                                                                                                                                                                     cGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACt                         <  1:787439/251‑1 (MQ=255)
                                                                                                                                                                                                                     cGTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACt                         <  1:668609/251‑1 (MQ=255)
                                                                                                                                                                                                                      gTATCCACATACTGAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTg                        >  1:437214/1‑251 (MQ=255)
                                                                                                                                                                                                                               tACTGAATTCCTG‑GGCTGTTAACTATTCCGGCCATCGGGCAAGACCGTTCTGACAAAGGGAAAGTGACGCGGTGGATGCAACTCGAACAAAAGTCGATAATGACTAGGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGGCTGGGGTTATATGCACGATCA‑AGATCCTGGAAAACGATCGAGACCGCGGATCATAGCCTAACCTGCGCAagag                <  2:175413/250‑1 (MQ=255)
                                                                                                                                                                                                                                aCTGAATTC‑TGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCgg                                        <  2:706581/225‑1 (MQ=255)
                                                                                                                                                                                                                                   gAATTC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGGAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtctt           <  1:272195/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTC‑TGTGGCTGATAACTCATCCTGCAATCGGGAAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACGCGAACAAAAGGCGATAATGACTAAGGCTGAAACATTCAGGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct         <  2:734447/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tC‑TGTGGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt       >  1:127808/1‑251 (MQ=255)
                                                                                                                                                                                                                                        c‑TGTGACTGATAGCTCATCCTGCAATCGTGCAAGACACTGCTGCCAAAGCGAAAGGAACACGGCGCACTCCACTCGGACAAAAGTCGATAATGAATCAGGATGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTG‑CATACCGAACCGCTGTCTGAGGTTCTATGCCCGATCCCACA‑CCTGCAAAACGATCGGGTCCCCGGCACATAGACTAACCTGCGCatgccatcttctttt      >  2:298953/1‑237 (MQ=255)
                                                                                                                                                                                                                                           gtgGCTGATAAATCATCCTGCAACCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCGGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGGCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTTCGCAAGATATCTTCTGt       <  2:866337/248‑1 (MQ=255)
                                                                                                                                                                                                                                            tgGCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GACCCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct   >  2:484701/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gCTGATAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc  >  1:397009/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
GGGGCCGCGACGTTGCTCGTCACTGCCGCTTGTGGCGTTTGCGTCACCGGTTTGGTGCCGACTTCAAAACGCAGCTGTGGGGCATCCGCTCCGCAGAAACTGGTTAGCAGTCCATTGATATTATTAAGGTACTTGTCCCGTACCCAATCGAGGACAAAACGGTTTGGCGCGTACAGGGCCAGCGTGTTATCGCTCAGTTCCGCCTGCAATGGGCGTATCCACATACTGAATTC‑TGTGGCTGGTAACTCATCCTGCAATCGGGCAAGACACTGCTGCCAAAGCGAAAGTGACACGGCGGACTCCACTCGAACAAAAGTCGATAATGACTAAGGCTGAAACATTCATGATTGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCAC‑GAACCGCTGTCTGCGGTTATATGCCCGATCA‑AGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTC  >  NC_000913/3883437‑3883923

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