Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,142,138 A→G 6.8% S97S (AGT→AGC arcD‑I ← arginine/ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,142,1380AG 6.8% 293.0 / 9.1 88S97S (AGT→AGCarcD‑Iarginine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base A (45/37);  new base G (3/3);  total (48/40)
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.63e-01

ACAGGCAGATCAGGGCGAACAGCACCAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGCCCAGCCGATCAGTACGGCCCCAACCCCTGCGCTGGCGGCCATGTTTTGTGGCAGCGAGA  >  AE015451/1141909‑1142349
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                   
aCAGGCAGATCAGGGCGAACAGCACCAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATc                                                                                                                                                                                                               <  1:804124/236‑1 (MQ=255)
       gATCAGGGCGAACAGCACCAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAg                                                                                                                                                                                                         <  1:853273/235‑1 (MQ=255)
        aTCAGGGCGAACAGCACCAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAg                                                                                                                                                                                                         <  2:485885/234‑1 (MQ=255)
           aGGGCGAACAGCACCAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGcccc                                                                                                                                                                                                     >  1:385460/1‑235 (MQ=255)
                  agacagccAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAg                                                                                                                                                                                             <  1:630539/230‑1 (MQ=255)
                     gCACCAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAgg                                                                                                                                                                                            >  2:870827/1‑234 (MQ=255)
                                 ccacTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAg                                                                                                                                                                              >  1:825210/1‑236 (MQ=255)
                                   acTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGATGGGTGTGTTGCCTTCGCCGAACACCGGGAAGTAGAAGCCCAGGGTGCTGAACAGCAGCACGAAATAGCCGACGTTGCCCAGCCAGGCGCTGATCCAGTAGCCCCAGGCCGAGGAAAAGCCCATGTACTc                                                                                                                                                                            <  1:194951/236‑1 (MQ=14)
                                      ttGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGc                                                                                                                                                                                                        >  2:630536/1‑205 (MQ=255)
                                       gggCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCGTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAg                                                                                                                                                                                                         <  2:131617/202‑1 (MQ=255)
                                                      ccAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAagcca                                                                                                                                                                  >  2:709112/1‑227 (MQ=255)
                                                           ttGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGc                                                                                                                                                                          >  2:573830/1‑214 (MQ=255)
                                                             gatgaAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCAc                                                                                                                                                                                                                    >  1:393409/1‑170 (MQ=255)
                                                             gatgaAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGgcatacgcat                                                                                                                                                  >  1:558347/1‑236 (MQ=255)
                                                             gatgaAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGgcatacgc                                                                                                                                                    <  2:558344/234‑1 (MQ=255)
                                                               tgaAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAg                                                                                                                                                                                                         <  2:738755/179‑1 (MQ=255)
                                                                      ggcTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAg                                                                                                                                                                                             >  1:845934/1‑184 (MQ=255)
                                                                           ccTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTGCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAg                                                                                                                                                                                             <  2:594851/179‑1 (MQ=255)
                                                                            cTTGAGGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTg                                                                                                                                                            <  2:268327/211‑1 (MQ=255)
                                                                              tGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCTAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATAGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCAc                                                                                                                                                                                                                    >  2:69459/1‑153 (MQ=255)
                                                                                 taccGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCACATAGTCGCCGAACCAGCCCTGgcataagca                                                                                                                                                   >  1:717083/3‑214 (MQ=255)
                                                                                 tGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATc                                                                                                                                                                                                               >  2:387133/1‑155 (MQ=255)
                                                                                  accGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTa                                                                                                                                                                                                          >  1:481544/2‑159 (MQ=255)
                                                                                  accGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTa                                                                                                                                                                                                          <  2:481541/158‑1 (MQ=255)
                                                                                       gcAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGAGGTTGCGCAGCCAGGCACTGATACAGTAGc                                                                                                                                                                                                        >  2:500204/1‑156 (MQ=255)
                                                                                         aGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGccgc                                                                                                                                                                                                     >  1:899558/1‑155 (MQ=255)
                                                                                         aGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGt                                                                                                                                >  1:422543/1‑226 (MQ=255)
                                                                                         aGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGGAGcc                                                                                                                                                                                                       >  2:747071/1‑155 (MQ=255)
                                                                                              cAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTgg                                                                                                                 <  1:709115/236‑1 (MQ=255)
                                                                                                  aaGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAg                                                                                                             <  1:682598/236‑1 (MQ=255)
                                                                                                     ggCACCGCCCACAACAGGATCGACGCGCCAATGATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGgcgtaggcat                                                                                                                                                  <  2:251619/195‑1 (MQ=255)
                                                                                                      gCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCCTGGCATACGCATACACCCCCCCATCCAg                                                                                                                                  >  1:229061/1‑211 (MQ=255)
                                                                                                      gCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCGTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGa                                                                                                                                                                       <  2:18803/174‑1 (MQ=255)
                                                                                                      gCACCCCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAgg                                                                                                                                                                                            <  1:401069/153‑1 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTc                                                                                                                                                                            >  1:619836/1‑167 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTc                                                                                                                                                                            <  2:619833/167‑1 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTgg                                                                                                                                                           <  2:152044/184‑1 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTgg                                                                                                                                                           >  1:152045/1‑184 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTgg                                                                                                                 >  1:566193/1‑226 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGc                                                                                                                >  1:418584/1‑227 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTg                                                                                                       >  2:476246/1‑236 (MQ=255)
                                                                                                        aCCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTg                                                                                                       >  1:426536/1‑236 (MQ=255)
                                                                                                          cGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAg                                                                                                                                                                              >  1:415306/1‑163 (MQ=255)
                                                                                                             ccACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAac                                                                                                   <  1:389794/235‑1 (MQ=255)
                                                                                                               acaacaGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCa                                                                                                                                   >  2:697788/1‑201 (MQ=255)
                                                                                                               acaacaGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACg                                                                                                >  1:495473/1‑236 (MQ=255)
                                                                                                                caacaGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCATTGGCATACGCATACACCCCACCATCCATGTCgg                                                                                                                             >  1:738014/1‑206 (MQ=255)
                                                                                                                  acaGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAg                                                                                             <  1:387136/236‑1 (MQ=255)
                                                                                                                         cGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGc                                                                                                                                                                                      >  2:617717/1‑140 (MQ=255)
                                                                                                                            cgcgCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATc                                                                                                                                                                                                               <  2:638081/112‑1 (MQ=255)
                                                                                                                              cgcCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCAc                                                                                                                                                                                                                    >  1:104155/1‑105 (MQ=255)
                                                                                                                              cgcCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCAc                                                                                                                                                                                                                    <  2:104155/105‑1 (MQ=255)
                                                                                                                              cgcCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGc                                                                                                                >  1:304296/1‑205 (MQ=255)
                                                                                                                              cgcCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACAGGaa                                                                                              >  2:842894/1‑223 (MQ=255)
                                                                                                                              cgcCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGaa                                                                                              <  1:842897/223‑1 (MQ=255)
                                                                                                                              cgcCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGAGCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGt                                                                                                          <  1:234652/211‑1 (MQ=255)
                                                                                                                               gcCAATTATGGCGGCTGGGGTGTTGACTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGGCGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAg                                                                                  <  2:400949/234‑1 (MQ=255)
                                                                                                                                 cAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAac                                                                                                   >  2:513926/1‑215 (MQ=255)
                                                                                                                                 cAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAgg                                                                                       <  2:422540/227‑1 (MQ=255)
                                                                                                                                   aTTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACAcc                                                                                                                                             <  1:75397/171‑1 (MQ=255)
                                                                                                                                   aTTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACAcc                                                                                                                                             >  2:75397/1‑171 (MQ=255)
                                                                                                                                         gcggcTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGtt                                                                                                                   <  1:1328/191‑1 (MQ=255)
                                                                                                                                         gcggcTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTc                                                                                                         <  2:304293/201‑1 (MQ=255)
                                                                                                                                              tGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCTACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGa                                                                                                     >  1:473679/1‑200 (MQ=255)
                                                                                                                                                gggTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATc                                                               >  1:860010/1‑236 (MQ=255)
                                                                                                                                                   tgtTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGcc                                                            <  1:38672/236‑1 (MQ=255)
                                                                                                                                                   tgtTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGcc                                                            <  1:498616/236‑1 (MQ=255)
                                                                                                                                                         cTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGCCCAGcc                                                       <  1:697791/235‑1 (MQ=255)
                                                                                                                                                                gAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAg                                                                                             >  1:388817/1‑190 (MQ=255)
                                                                                                                                                                     ccGGGAAGTAGAAGCCCAGGGTGCTGAACAGCAGCACGAAATAGCCGACGTTGCCCAGCCAGGCGCTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAg                                                                                             >  2:275196/2‑185 (MQ=255)
                                                                                                                                                                     ccGGGAAGTAGAAGCCCAGGGTGCTGAACAGCAGCACGAAATAGCCGACGTTGCCCAGCCAGGCGCTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAg                                                                                             <  1:275199/184‑1 (MQ=255)
                                                                                                                                                                                  gCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGCCCAGCCGATCAGTACGGCCCCAACCCCTgcga                             >  2:645962/1‑235 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATc                                                                                                                                                                                                               <  2:860546/50‑1 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATc                                                                                                                                                                                                               >  1:860549/1‑50 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATGTAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAg                                                                                                                                  >  2:458424/1‑127 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCa                                                                                                                                   <  2:685031/126‑1 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCa                                                                                                                                   >  1:685034/1‑126 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACCGCGGTGATCCCCCAGc                                                        >  2:682192/1‑201 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAATAGCCGACGTTGCCCAGCCAGGCGCTGATCCAGTAGCCCCAGGCCGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAg                                                                                             <  1:167687/164‑1 (MQ=255)
                                                                                                                                                                                          tGCTGAACAGCAGCACGAAATAGCCGACGTTGCCCAGCCAGGCGCTGATCCAGTAGCCCCAGGCCGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAg                                                                                             >  2:167686/1‑164 (MQ=255)
                                                                                                                                                                                           gCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGcc                                                                                                               <  2:668990/145‑1 (MQ=255)
                                                                                                                                                                                            cTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGTCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTg                                                                       >  2:658339/1‑184 (MQ=255)
                                                                                                                                                                                                       caCGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGa                                                                 <  2:406331/179‑1 (MQ=255)
                                                                                                                                                                                                       caCGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGa                                                                 >  1:406334/1‑179 (MQ=255)
                                                                                                                                                                                                          gAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGCCCAGCCGATCAGTACGGCCCCAACCCCTgcgc                             <  1:41114/212‑1 (MQ=255)
                                                                                                                                                                                                              tAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGCCCAGCCGATCAGTACGGCCCCAACCCCTGCGCTGGCGGCCATGTTTTGTGGCAGCgaga  <  1:932075/235‑1 (MQ=255)
                                                                                                                                                                                                                      gTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGGCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCg                                                                     >  2:288383/1‑160 (MQ=255)
                                                                                                                                                                                                                                  ggCGCTGATCCAGTAGCCCCAGGCCGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCAGGCCGACTGCGGTGATCGCCCAGCCGATCAGTACggg                                           >  1:152538/1‑173 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                   
ACAGGCAGATCAGGGCGAACAGCACCAGCGGCACCACTTTGGCCACAGTGGTGACCAAGTTGATGAAGGCGGCTTCCTTGATGCCGCGCAGTACCAGGAAGTGCACCGCCCACAACAGGATCGACGCGCCAATTATGGCGGCTGGGGTGTTGCCTTCACCGAAGATCGGGAAGAAGTAGCCCAGGGTGCTGAACAGCAGCACGAAGTAGCCGACGTTGCCCAGCCAGGCACTGATCCAGTAGCCCCAGGCTGAGGAAAAGCCCATATAGTCGCCGAAGCCAGCCTTGGCATACGCATACACCCCACCATCCAGGTCGGGCTTGCGGTTGGCCAGGGTCTGGAACACGAAGGCCAGGGTCAGCATGCCGACTGCGGTGATCGCCCAGCCGATCAGTACGGCCCCAACCCCTGCGCTGGCGGCCATGTTTTGTGGCAGCGAGA  >  AE015451/1141909‑1142349

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