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.100.0% 521.9 / NA 118intergenic (+105/+163)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (35/83);  total (35/83)
*AP0123062,073,8480T.99.2% 517.3 / ‑2.2 118intergenic (+106/+162)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (35/82);  minor base A (0/1);  total (35/83)
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 129 total aligned reads displayed.

CTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGTGCTGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTT  >  AP012306/2073698‑2073987
                                                                                                                                                     ||                                                                                                                                           
cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCg                                                                                                                                              <  2:3435180/150‑1 (MQ=32)
cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCg                                                                                                                                              >  2:1253144/1‑150 (MQ=32)
 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             >  2:1906739/1‑149 (MQ=14)
 ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgc                                                                                                                                             >  1:2442338/1‑149 (MQ=14)
   ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                                                           >  2:960242/1‑147 (MQ=1)
   ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                                                           >  2:4388719/1‑147 (MQ=1)
                                                                                                    gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACACCGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGa                                       <  1:1002218/150‑1 (MQ=1)
                                                                                                    gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGa                                       >  2:3741514/1‑150 (MQ=0)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      >  1:1166677/1‑150 (MQ=2)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  2:4582982/150‑1 (MQ=2)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  2:42930/150‑1 (MQ=2)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  2:4275011/150‑1 (MQ=2)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      >  2:1865114/1‑150 (MQ=2)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      >  1:4283769/1‑150 (MQ=2)
                                                                                                     aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa                                      <  1:2714714/150‑1 (MQ=2)
                                                                                                      aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa                                     <  2:4576822/150‑1 (MQ=255)
                                                                                                       ccTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt                                    >  2:4265421/1‑150 (MQ=255)
                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   >  2:2240104/1‑150 (MQ=255)
                                                                                                        cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt                                   >  2:4537262/1‑150 (MQ=255)
                                                                                                         tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTg                                  >  1:3631117/1‑150 (MQ=255)
                                                                                                         tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTg                                  <  1:3304006/150‑1 (MQ=255)
                                                                                                         tGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTg                                  >  1:3070479/1‑150 (MQ=255)
                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:4143260/150‑1 (MQ=255)
                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:2073069/150‑1 (MQ=255)
                                                                                                          gcagcaGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                 <  2:3074861/150‑1 (MQ=255)
                                                                                                           cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGcc                                >  2:4182700/1‑150 (MQ=255)
                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               <  1:1562798/150‑1 (MQ=255)
                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               >  1:2991602/1‑150 (MQ=255)
                                                                                                            agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt                               <  1:1759752/150‑1 (MQ=255)
                                                                                                             gcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTg                              >  2:1301735/1‑150 (MQ=255)
                                                                                                               agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  1:21258/150‑1 (MQ=255)
                                                                                                               agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            <  2:53189/150‑1 (MQ=255)
                                                                                                               agGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGAt                            >  1:2525117/1‑150 (MQ=255)
                                                                                                                ggTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATg                           <  1:3357239/150‑1 (MQ=255)
                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          <  1:2824326/150‑1 (MQ=255)
                                                                                                                 gTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgc                          >  1:3646630/1‑150 (MQ=255)
                                                                                                                  tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         >  2:2135017/1‑150 (MQ=255)
                                                                                                                  tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         <  1:3423567/150‑1 (MQ=255)
                                                                                                                  tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         <  1:3150809/150‑1 (MQ=255)
                                                                                                                  tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         <  1:200904/150‑1 (MQ=255)
                                                                                                                  tttGAGTTGAGCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACCGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcg                         >  1:3798871/1‑150 (MQ=14)
                                                                                                                   ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  1:3849874/150‑1 (MQ=255)
                                                                                                                   ttgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATgcgc                        <  2:830529/150‑1 (MQ=255)
                                                                                                                    tgagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCt                       <  1:1204205/150‑1 (MQ=255)
                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGTCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  2:2049217/150‑1 (MQ=14)
                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      <  1:4251015/150‑1 (MQ=255)
                                                                                                                     gagttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                      >  1:1600642/1‑150 (MQ=255)
                                                                                                                      agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                     >  1:1301224/1‑150 (MQ=255)
                                                                                                                      agttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCCTc                     <  2:1585687/150‑1 (MQ=14)
                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    >  2:2622124/1‑150 (MQ=255)
                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  2:4346761/150‑1 (MQ=255)
                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  2:514296/150‑1 (MQ=255)
                                                                                                                       gttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCa                    <  1:2943119/150‑1 (MQ=255)
                                                                                                                        ttgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAc                   <  1:2465004/150‑1 (MQ=255)
                                                                                                                         tgagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACt                  <  1:1169642/150‑1 (MQ=255)
                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 <  2:526790/150‑1 (MQ=255)
                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 <  1:1366158/150‑1 (MQ=255)
                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 <  1:4486835/150‑1 (MQ=255)
                                                                                                                          gagCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACtt                 >  2:288806/1‑150 (MQ=255)
                                                                                                                           agCGCAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTa                <  2:2693146/150‑1 (MQ=255)
                                                                                                                            gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  1:2960393/150‑1 (MQ=255)
                                                                                                                            gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  2:1306491/150‑1 (MQ=255)
                                                                                                                            gcgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAg               <  2:4238145/150‑1 (MQ=255)
                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              <  2:3859442/150‑1 (MQ=255)
                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              <  1:2279374/150‑1 (MQ=255)
                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              <  1:3966988/150‑1 (MQ=255)
                                                                                                                             cgcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGc              <  1:967009/150‑1 (MQ=255)
                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  2:841796/150‑1 (MQ=255)
                                                                                                                              gcAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCa             <  2:3195612/150‑1 (MQ=255)
                                                                                                                               cAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            <  1:1432135/150‑1 (MQ=255)
                                                                                                                               cAGCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAg            >  2:4201311/1‑150 (MQ=255)
                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  1:487021/150‑1 (MQ=255)
                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  2:3222337/150‑1 (MQ=255)
                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  2:3547626/150‑1 (MQ=255)
                                                                                                                                 gCGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGac          <  1:3314430/150‑1 (MQ=255)
                                                                                                                                  cGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGact         <  2:2567342/150‑1 (MQ=255)
                                                                                                                                  cGAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGact         <  2:480165/150‑1 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  2:1201221/150‑1 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:2603291/150‑1 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  2:1657157/150‑1 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:932620/150‑1 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  2:3851934/150‑1 (MQ=255)
                                                                                                                                   gAAACAACGGAGCCGCT‑‑AGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGacta        <  1:3087617/150‑1 (MQ=17)
                                                                                                                                    aaaCAACGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactac       >  1:1860825/1‑150 (MQ=255)
                                                                                                                                     aacaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactact      <  2:1431467/150‑1 (MQ=255)
                                                                                                                                     aacaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactact      >  2:1413550/1‑150 (MQ=255)
                                                                                                                                     aacaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactact      <  1:3156272/150‑1 (MQ=255)
                                                                                                                                     aacaacGGAGACGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGactact      >  1:2164420/1‑150 (MQ=17)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  2:4404554/150‑1 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  1:3558486/150‑1 (MQ=255)
                                                                                                                                      acaacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGActacta     <  2:4005137/150‑1 (MQ=255)
                                                                                                                                       caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCCCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  1:4341376/150‑1 (MQ=18)
                                                                                                                                       caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  1:681643/150‑1 (MQ=255)
                                                                                                                                       caacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAt    <  2:2749409/150‑1 (MQ=255)
                                                                                                                                        aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  1:3529696/150‑1 (MQ=255)
                                                                                                                                        aacGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAtt   <  2:1658064/150‑1 (MQ=255)
                                                                                                                                         acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  >  1:3472823/1‑150 (MQ=255)
                                                                                                                                         acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  2:1005416/150‑1 (MQ=255)
                                                                                                                                         acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  2:3943570/150‑1 (MQ=255)
                                                                                                                                         acGGAGCCGCTC‑‑GCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTAttt  <  1:3641634/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 < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: