Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,628,141 A→G 10.8% E276G (GAG→GGG)  nadR → trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,1410AG10.8% 210.2 / 3.3 79E276G (GAG→GGG) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (20/48);  new base G (7/2);  total (29/50)
Fisher's exact test for biased strand distribution p-value = 7.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGG  >  NC_000913/4627905‑4628374
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           
tgATCCGAAACTTTCCTTTTTGAGTAACAGCGGTGCGCAGATCCTCGAAAACCCTTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATgaga                                                                                                                                                                                                                                           <  2:111263/238‑1 (MQ=255)
 gATCCGGAACGTTCCTTTAGGAGTATGATCGGTTCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATgaga                                                                                                                                                                                                                                           <  1:118389/238‑1 (MQ=255)
   tCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATc                                                                                                                                                                                                                                         >  1:860534/1‑238 (MQ=255)
   tCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACACCCAGTGCGGGGGAAAATGTCCGCGATTATGTCTTTTCACACCCCGGCGGGGATgagaaag                                                                                                                                                                                                                                        >  2:517722/1‑236 (MQ=255)
     cGAAACGTTCCTTTATGGGTATACGCGGTGCGAAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                                                                                                                                                                                                      <  2:480893/239‑1 (MQ=255)
     cGAAACGTGCTTTTTTGTGTAGCAGCGGTTCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAAGATATTCCTACGGAATTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                                                                                                                                                                                                      <  2:802188/239‑1 (MQ=255)
     cGAAACGTCCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                                                                                                                                                                                                      <  1:532841/239‑1 (MQ=255)
        aaCGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTg                                                                                                                                                                                                                                   <  1:370276/239‑1 (MQ=255)
         aCGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGc                                                                                                                                                                                                                                  <  1:203461/239‑1 (MQ=255)
             aCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGGGTTGGGATCGCAtgtcag                                                                                                                                                                                                                                >  1:724089/1‑232 (MQ=255)
                 ttATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCt                                                                                                                                                                                                                          <  1:888792/239‑1 (MQ=255)
                    tggttATCAGCGGGGCGCAGCTCCGCGAAACCCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGAACCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGAc                                                                                                                                                                                                                       <  2:556168/235‑1 (MQ=255)
                    tggatATCAGCGGTCCGCAGGTCCGGGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGAc                                                                                                                                                                                                                       <  1:338260/235‑1 (MQ=255)
                    tGAGGATCAGCGGGTCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATGCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGAc                                                                                                                                                                                                                       <  1:394947/239‑1 (MQ=255)
                         aTCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCCGGATTATGTCTTTTCACACCTCGGCGGGGATGGAATGGCATTGCAGTATTTTgctttagt                                                                                                                                                                                                                  >  1:874020/1‑232 (MQ=255)
                               ggTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATaaaa                                                                                                                                                                                                             <  1:717794/238‑1 (MQ=255)
                                   cgcAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGCATTCTGACTCCGTAaaaaaccg                                                                                                                                                                                                         >  1:354693/1‑234 (MQ=255)
                                    gcAGATCCGCGAAAACCCTGTCCGCTACTGGTAATATGTTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATTTCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATcgcg                                                                                                                                                                                                        <  1:57507/238‑1 (MQ=255)
                                          ccGCGAAACCCCGTTCCCGTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGcc                                                                                                                                                                                                <  1:779680/239‑1 (MQ=255)
                                              caaaaCCCGTTCCGTTTGGGGGAATACATTCCTACGGAAGGGAAGCCGTGTTTTGGGCTTACCGTGGGGATCCTGGGCGGCGTGTCGGGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACAACGATAAAAGCGCGCTGGGCCACGc                                                                                                                                                                                            <  2:150226/238‑1 (MQ=255)
                                                aatcccGTTCCGCTATTGGGAATATTTTCCTCCCGAAGGGAAGCCGTTTTTTGTGCGTACCGTTGCGGTTCTTGGCGGCGAGTCGAGCGGTAAATCCCCCCTGGTAAACAAACTTGCCAATATCTTCAACACCCCCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTc                                                                                                                                                                                          <  2:717319/236‑1 (MQ=255)
                                                aaccccTTTCCGCTACTGGGAATTTATTCCTACCGAAGTGAAGGCGTTTTGTGGGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTc                                                                                                                                                                                          <  2:167029/236‑1 (MQ=255)
                                                aaaCCCGTTCCGCTACTGGGAATCTTTTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTc                                                                                                                                                                                          <  1:879231/239‑1 (MQ=255)
                                                 agcccTTTCCGTTTGTGGGAATATTTGCCCACCGAGGGGAAGCCGTTTTTTGTGCGCACCGTGGCGATCCTTGGCGGGGGGTCGAGCGTTAAACCCCCCCTGGTAAACAAACCTGCCAATATCCTCAACCCCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCa                                                                                                                                                                                         <  2:596481/237‑1 (MQ=255)
                                                        tCCGCTACTGGGAATATTTTCCTACCGAAGTGAAGCCGGTTTTTGTGCGTACCGTGGCGATTCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACAt                                                                                                                                                                                   <  1:688599/238‑1 (MQ=255)
                                                             tACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTCCCGGGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGAtt                                                                                                                                                                              <  1:237621/238‑1 (MQ=255)
                                                               cTGGGAAATTTTTCCCTCCGAGGTGAAGCCGTTTTTTTGGCTTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCCCCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTg                                                                                                                                                                           <  2:860534/239‑1 (MQ=255)
                                                                 gggAATTTTTTCCTACCGAAGTGAAGCCGTTTTTTGGTCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATATCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCg                                                                                                                                                                         <  2:206509/239‑1 (MQ=255)
                                                                   gAATATATTCCTGACGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGACCCTTGGCGGCGAGTAGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACTCCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCgg                                                                                                                                                                        <  1:800381/238‑1 (MQ=255)
                                                                         atTCCTACCGAAGTGAAGCCGTTTTTTGGGCGTACCGTGGGGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAata                                                                                                                                                                 <  1:803504/239‑1 (MQ=255)
                                                                          ttCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAatat                                                                                                                                                                <  1:146669/239‑1 (MQ=255)
                                                                              tACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATAGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCCGGCTTCTGACTACGATAAAATCGCGCTGGGCCCCGCCCAATACATTGATTTTTCGGGGAAAATTGccc                                                                                                                                                            >  1:325307/1‑238 (MQ=255)
                                                                              tACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGCAAATCCACCCTGGTAACCAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGGGGATGGGTCGCATTTCCGGATTTTGACCTCCGTAAAAAAGCGCTGGGGCCACGCCCATTACTTTCTTTTTGGGGTaaaaatggca                                                                                                                                                             >  2:632000/1‑233 (MQ=255)
                                                                                  gAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCCGTATTCTGACTACGCTAAAATCCCGCTGGGCCACGCCCCATACATTGATTTTTCGGTGAAATATGCCAATaa                                                                                                                                                        >  1:7759/1‑239 (MQ=255)
                                                                                  gAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGGTAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTTCGGTGAAATATGCCTATaa                                                                                                                                                        >  1:744662/1‑239 (MQ=255)
                                                                                   aaGTGAAGCCGGTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATaaa                                                                                                                                                       <  1:604063/239‑1 (MQ=255)
                                                                                    aGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCCCCCTGGTAAACAAACTTGCCAATATCTTCAACCCCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATaaa                                                                                                                                                       <  1:142808/238‑1 (MQ=255)
                                                                                     gTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGGTGAGATCGCCTTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCCATaaata                                                                                                                                                     >  1:711277/1‑237 (MQ=255)
                                                                                        aaGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGCTTATGTCTTTTCACACCTCGGCGGGGATGAGATCGCATTGCCGTATTTTGACTACGCTAAAATCGCGCTGGGACACGCGCAAAACATTTTTTTTTCGGGGAAAAATGCCAAAAAAGTGGCa                                                                                                                                                 >  2:593816/1‑239 (MQ=255)
                                                                                              ttttttGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCAGTATTCTGAc                                                                                                                                                                                                                       >  1:788267/1‑165 (MQ=255)
                                                                                              ttttttGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGAc                                                                                                                                                                                                                       <  2:788267/165‑1 (MQ=255)
                                                                                                   tttgCGTCCCGTGGGTATCCTTGGGGGCGATTCGGTCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCGTCAACCCCCCCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATa                                                                                                                                       <  2:354693/237‑1 (MQ=255)
                                                                                                        gttccGTGGCGATCCTTGGGGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGAt                                                                                                                                  <  1:517722/236‑1 (MQ=255)
                                                                                                                                aGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGGTTACGATTTTTTTACCACTCAGGCGTTTTGc                                                                                                          >  1:710801/1‑239 (MQ=255)
                                                                                                                                   cGAGCGGTAAATCCCCCCTGGGAAACAAACTTGGCAATATCTTCAACTCCACCAGTGCGTGGGAATATGGCCGGGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCaa                                                                                                        <  1:293059/238‑1 (MQ=255)
                                                                                                                                    gggCGGTAATTCCCCCCTGGTAAACAAACTTGCCAATATCTTCAACCCCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCaaaa                                                                                                      <  2:267120/237‑1 (MQ=255)
                                                                                                                                           aaaTCCACCCTGGTAAACAAACTTGCCAATTTCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACg                                                                                                <  1:622323/238‑1 (MQ=255)
                                                                                                                                           aaaTCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTAAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACg                                                                                                <  1:804351/238‑1 (MQ=255)
                                                                                                                                            aaTCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCAGGATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTTTGGTGGAATATGCCAATAAAGGGGCATTTTTCGAAACCGCTTTTTTCCCCCACCAGGCGCTTTTCaaaaaaaaagaa                                                                                              >  2:497752/1‑232 (MQ=255)
                                                                                                                                            aaTCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGa                                                                                               <  1:694111/238‑1 (MQ=255)
                                                                                                                                                   ccTGGTAAACAAACTTGCCAATATCTTCAACGCCACCAGTGCGGTGGAATAGGGCCGCGATTATGTCGTTTCACAACTCGGCGGTGATGAGCTCGCATTGCTGTATTATGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCggg                                                                                       <  1:401619/239‑1 (MQ=255)
                                                                                                                                                      ggTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCCCAATACATTGATTTTACGGGGAAATATTCCAATAAAATGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTTTGCAAAAAGTACGAAGGGCGGGAAc                                                                                    >  1:693889/1‑239 (MQ=255)
                                                                                                                                                                 ttGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATACCATTGCAGGATTTTGACTACGATAAAATCGGGCTGGGCCACGCCTCATACATTGATTTTTCGGTGGAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTTGCACACCTCATGCGTTCTGCAAAAAGAACAAAGAGCGGGAACacaccttcggg                                                                         >  2:612078/1‑229 (MQ=255)
                                                                                                                                                                    ccAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGCTTTTGTCCCCACCCAGGCGTTTTGCCAAAAGAACGAAGGGCGGGAACAACCGTTCGTGCAg                                                                     >  1:825410/1‑239 (MQ=255)
                                                                                                                                                                       attttttCCACCCCCCCAGTGGGGTGGAATTTGGCCGCGATTAGGTCTTTTCACACCTGGGCGGGGATGAGATCGCATTGCATTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGgcg                                                                  <  2:711277/234‑1 (MQ=255)
                                                                                                                                                                       atatCTTCAACCCCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGgcg                                                                  <  1:556272/239‑1 (MQ=255)
                                                                                                                                                                       atatCTTCAACACCACCAGGGCGTGGGAATATGGCCGTGATTATGGCTGTTCACACCTCGGCGGGAATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGgcg                                                                  <  1:801709/239‑1 (MQ=255)
                                                                                                                                                                         atCTTCAACCCCCCCAGTGCGGGGGAATTTGGCCAGGATTTTGTCTTTTCACCCCGCGGCGGTGATGAGATCGCAATGCAGTATTCTGACTACGATAAAATCGCGCTGGGACACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCt                                                                <  2:450417/239‑1 (MQ=255)
                                                                                                                                                                           cTTCAACCCCACCAGGGCAGGGGAATATGGCCGGGATTATGGCTTTACACAGCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtga                                                              <  1:4852/239‑1 (MQ=255)
                                                                                                                                                                                acccccccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCCCCCCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattgat                                                         <  1:618210/233‑1 (MQ=255)
                                                                                                                                                                                    caccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAACACGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCGATAAAGGGGCATTTATCGTAACCGATTTTGTCCCAACTCAGGCGTTTAGCAAAAAGAACGAAGGGCGGGGAAACCCGTTCTAGGAGGGGCTGATTGATGaaa                                                     >  2:29139/1‑238 (MQ=255)
                                                                                                                                                                                      ccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGAATACGATAAAATCGCGCTGGGCCACGCTCAATAAATTGATTTTGCGGTGAAATATGCCAATAAAGATGCACTTTTCGATACCGATTTTGTTCCCACCCAGGCGTTCTTCAAAAAGTACGACGGGCGGGCAAATCCGTTCCTTCAGGCGCCGATTGATGaaaaa                                                   >  2:35983/1‑236 (MQ=255)
                                                                                                                                                                                         gTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGGGAAATATGCCAATAAAATGGCATTTATCGATACCGCTTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCGGATTGTTGAATACCGt                                                >  1:516376/1‑239 (MQ=255)
                                                                                                                                                                                           gCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGtt                                               <  1:678931/238‑1 (MQ=255)
                                                                                                                                                                                            cGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCCCTGGGCCACGCCCAATACATTGATTTTTCGGGGAAATAAGCCCATAAAGCGGCCTTTTTCGATACCCGTTTTTTCACCACCCAGGCGTTCTTCAAAAAGAAAGAAGGGGGGGAACCTACGTTCTTGGCAGGGCTGGTTGGAGAAAACCGcttg                                             >  2:164005/1‑238 (MQ=255)
                                                                                                                                                                                                 gAATTTGGCTGCGATTATGTCTTTTCACACCTAGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATc                                         <  1:468874/238‑1 (MQ=255)
                                                                                                                                                                                                    tatGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCCCGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTTTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGGACATCCGTTCGTGCCGGGCCTGGTTGATGAAAACCGTTTTGAACTggg                                     >  1:356716/1‑238 (MQ=255)
                                                                                                                                                                                                       ggCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGCGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCCCAATACTTTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGTTTTTTTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGGAAATCCGTTCGTGCCGGCGCTGCTTAATGAATACCTGTTCGGTCTGGTgt                                   >  1:780691/1‑237 (MQ=255)
                                                                                                                                                                                                        gCCGCGTTTTTTTCTTTTCCCACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATc                                 <  2:7759/239‑1 (MQ=255)
                                                                                                                                                                                                        gCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTCGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGGGCAGGCGCTGATTGATGAAAACCGGTTCGATCTGGGGAt                                  >  1:361521/1‑238 (MQ=255)
                                                                                                                                                                                                          cgggATTATGTCTTTTACAGCCTCGGCGGTGATGGGATCGTATTGAAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTAAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATcc                                <  1:272250/235‑1 (MQ=255)
                                                                                                                                                                                                           gcgAATATGTCTTTTCCCACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGAGAAAATCGCGCTGGGCCGCGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCctg                              <  1:632000/239‑1 (MQ=255)
                                                                                                                                                                                                           gagATTATGGCTTTCCACACCTCGGCGGTGATGAGGTCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCctg                              <  1:180458/237‑1 (MQ=255)
                                                                                                                                                                                                            cgATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGCATGGCGGGAACATCCGTTCGGGCAGGCGCTGATTGATGAAAACCGTTTCGAttctgaatcct                                >  1:185192/1‑226 (MQ=255)
                                                                                                                                                                                                               tttttttttttCACCCCGCGGCGGTTGTGAGATCCCATTGCAATATTCTGACGTAGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCctgctg                           <  2:459171/231‑1 (MQ=255)
                                                                                                                                                                                                                  tGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATAACCGTTTTGTCACCACTCAGGCGTTCTGCAAAAAGTAAGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATAACGTTTCGATCTGGTGAACCCGCTGgaaa                       >  2:202808/1‑237 (MQ=255)
                                                                                                                                                                                                                     cTTTTCCCACCTCGGCGGTGATGGGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGaaca                    <  1:457071/239‑1 (MQ=255)
                                                                                                                                                                                                                       tttCACACCTCGGCGGTGATGGGAGCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATAAATTGATTTTGCGGTGAAATATGCCAATAAAGGAGCATTTATCGATGACGATTTTGTCACCACTCAGGCGCTTTGCAAAAAGTACAATGGGCGGGAAAATCCGTTCCGGCAGCTGCGGGTTGTTGAAAACCGCTTGGACCTGGGGTGCCTGCGGGGGaaaaga                  >  2:354118/1‑235 (MQ=255)
                                                                                                                                                                                                                       tttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGGGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGGTCTGGGGATCCTTCTGGGGaacaaa                  >  1:329418/1‑238 (MQ=255)
                                                                                                                                                                                                                       ttgcaccCCTCGGCGGTGCGGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTAAATACATTGATTTTGCGGTTAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGaacaac                  <  1:317100/236‑1 (MQ=255)
                                                                                                                                                                                                                         tCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGGGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGGGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGGGATCCTTCTGGGGAacaaca                 >  1:204781/1‑238 (MQ=255)
                                                                                                                                                                                                                               ccTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTAAGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGAGGTGAAATATGCCAATAAAGTGGCATTTATCGGTACCGATTTTGTTACCACTCAGGCGATCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCATGCGCTGATTGATGCATACCATTTCGATCTGGTGATCCTGCTTGAGAGAAACACGcagc           >  1:571368/1‑236 (MQ=255)
                                                                                                                                                                                                                                         gtatggGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGTGGGTCCACGGTCAATACATGGATTTGGCGGTGAAATATGCCAATAAAGTGGAATTTATCGATACCGATTTGGTCCCCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCgg  <  1:148816/237‑1 (MQ=255)
                                                                                                                                                                                                                                          tatgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCgg  <  1:354118/237‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           
CGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGG  >  NC_000913/4627905‑4628374

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