Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 2,073,848 Δ2 bp intergenic (+106/+161) lysV → / ← xapR tRNA‑Lys/DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,073,8470G.99.0% 448.3 / ‑3.0 99intergenic (+105/+163)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (32/66);  total (32/67)
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
*AP0123062,073,8480T.100.0% 451.8 / NA 99intergenic (+106/+162)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (32/67);  total (32/67)

Only 100 of 105 total aligned reads displayed.

CTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGTGCTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTT  >  AP012306/2073698‑2073987
                                                                                                                                                     ||                                                                                                                                           
cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCg                                                                                                                                              >  2:3173399/1‑150 (MQ=32)
cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCg                                                                                                                                              >  1:3540561/1‑150 (MQ=32)
  tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcg                                                                                                                                            >  1:881022/1‑148 (MQ=11)
  tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACAGAGCCGCTCgcg                                                                                                                                            >  2:2722452/1‑148 (MQ=11)
                                                                                                    gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGa                                       >  2:3548956/1‑150 (MQ=0)
                                                                                                    gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGa                                       >  1:1946800/1‑150 (MQ=0)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      >  1:2972802/1‑150 (MQ=2)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      >  1:1839638/1‑150 (MQ=2)
                                                                                                     aaCCTGCAGCAGGGTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      >  2:36660/1‑150 (MQ=11)
                                                                                                      aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     >  2:2744994/1‑150 (MQ=255)
                                                                                                      aCCTGCAGCAGGTTGGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     <  1:1663821/150‑1 (MQ=11)
                                                                                                       ccTGCATCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    >  1:3527912/1‑150 (MQ=11)
                                                                                                       ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    >  1:3521952/1‑150 (MQ=255)
                                                                                                       ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    >  1:2814475/1‑150 (MQ=255)
                                                                                                       ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    <  1:1878104/150‑1 (MQ=255)
                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  2:1902947/150‑1 (MQ=255)
                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  2:3306943/150‑1 (MQ=255)
                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   >  1:1893852/1‑150 (MQ=255)
                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   <  1:1088226/150‑1 (MQ=255)
                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCCATTTGCGAAAtt                                   <  2:4031656/150‑1 (MQ=11)
                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  1:968127/150‑1 (MQ=255)
                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:703231/150‑1 (MQ=255)
                                                                                                           cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGcc                                >  2:2942310/1‑150 (MQ=255)
                                                                                                           cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGcc                                <  2:275519/150‑1 (MQ=255)
                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAGGCGAGTCATCCTGCACGACCCAACAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               >  2:107621/1‑150 (MQ=11)
                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               >  2:3655421/1‑150 (MQ=255)
                                                                                                             gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTg                              <  2:3670354/150‑1 (MQ=255)
                                                                                                             gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTg                              >  2:279188/1‑150 (MQ=255)
                                                                                                             gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTg                              <  1:645799/150‑1 (MQ=255)
                                                                                                             gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTg                              <  2:1191135/150‑1 (MQ=255)
                                                                                                              cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             <  2:2077315/150‑1 (MQ=255)
                                                                                                              cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             <  1:2172937/150‑1 (MQ=255)
                                                                                                              cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             <  2:2831389/150‑1 (MQ=255)
                                                                                                              cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa                             >  1:988302/1‑150 (MQ=255)
                                                                                                               agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  2:2141898/150‑1 (MQ=255)
                                                                                                                ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  1:1338384/150‑1 (MQ=255)
                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          <  1:3413015/150‑1 (MQ=255)
                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          <  1:1727089/150‑1 (MQ=255)
                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          >  1:500752/1‑150 (MQ=255)
                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          <  1:2145345/150‑1 (MQ=255)
                                                                                                                   ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  1:1839688/150‑1 (MQ=255)
                                                                                                                   ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  1:3741878/150‑1 (MQ=255)
                                                                                                                    tgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCt                       >  2:1892870/1‑150 (MQ=255)
                                                                                                                    tgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCt                       <  2:760766/150‑1 (MQ=255)
                                                                                                                    tgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCt                       <  2:462620/150‑1 (MQ=255)
                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCTACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  1:3824708/150‑1 (MQ=14)
                                                                                                                      agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     >  1:3621889/1‑150 (MQ=255)
                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  2:3603126/150‑1 (MQ=255)
                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  2:2334767/150‑1 (MQ=255)
                                                                                                                        ttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                   <  1:368054/150‑1 (MQ=255)
                                                                                                                        ttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                   <  2:3947137/150‑1 (MQ=255)
                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  1:140121/150‑1 (MQ=255)
                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  >  1:2008783/1‑150 (MQ=255)
                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 <  2:3267175/150‑1 (MQ=255)
                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 <  2:1762611/150‑1 (MQ=255)
                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 >  1:681081/1‑150 (MQ=255)
                                                                                                                           agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGAGGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGGGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                >  1:2612527/1‑150 (MQ=17)
                                                                                                                           agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                <  1:2647443/150‑1 (MQ=255)
                                                                                                                           agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                <  1:1308427/150‑1 (MQ=255)
                                                                                                                           agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                <  2:3207091/150‑1 (MQ=255)
                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              <  1:3690863/150‑1 (MQ=255)
                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              >  1:3527756/1‑150 (MQ=255)
                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              >  2:665254/1‑150 (MQ=255)
                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  2:3508429/150‑1 (MQ=255)
                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  1:598221/150‑1 (MQ=255)
                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  1:3611797/150‑1 (MQ=255)
                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  1:3900420/150‑1 (MQ=255)
                                                                                                                              gcAGCGAAACAACGGAGCCGCTGGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  1:3638022/150‑1 (MQ=17)
                                                                                                                               cAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            <  2:1391294/150‑1 (MQ=255)
                                                                                                                               cAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            >  1:3735200/1‑150 (MQ=255)
                                                                                                                                aGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGa           <  1:1239916/150‑1 (MQ=255)
                                                                                                                                aGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGa           <  1:1978284/150‑1 (MQ=255)
                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  2:1853385/150‑1 (MQ=255)
                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          >  1:362640/1‑150 (MQ=255)
                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          >  1:3713479/1‑150 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:3807241/150‑1 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        >  2:2076070/1‑150 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        >  2:3178871/1‑150 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:3574087/150‑1 (MQ=255)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  1:3418519/150‑1 (MQ=255)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  1:1275447/150‑1 (MQ=255)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  2:2197219/150‑1 (MQ=255)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  1:1013892/150‑1 (MQ=255)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  2:3451048/150‑1 (MQ=255)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       <  2:255740/150‑1 (MQ=255)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       >  1:804287/1‑150 (MQ=255)
                                                                                                                                     aacaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactact      <  2:2650057/150‑1 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  2:2645766/150‑1 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  1:2983998/150‑1 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  2:364098/150‑1 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  1:2552830/150‑1 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     >  1:1983007/1‑150 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     >  2:2966580/1‑150 (MQ=255)
                                                                                                                                       caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    >  2:115765/1‑150 (MQ=255)
                                                                                                                                        aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   >  1:1842315/1‑150 (MQ=255)
                                                                                                                                        aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  1:3772764/150‑1 (MQ=255)
                                                                                                                                        aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  1:1118670/150‑1 (MQ=255)
                                                                                                                                        aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  1:1037029/150‑1 (MQ=255)
                                                                                                                                         acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  1:2817519/150‑1 (MQ=255)
                                                                                                                                         acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  2:63576/150‑1 (MQ=255)
                                                                                                                                                     ||                                                                                                                                           
CTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGTGCTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTT  >  AP012306/2073698‑2073987

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