Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 602,835 G→A A293T (GCT→ACT)  pheP → phenylalanine:H(+) symporter PheP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913602,8350GA97.1% 196.9 / ‑1.9 70A293T (GCT→ACT) phePphenylalanine:H(+) symporter PheP
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (34/32);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-02

GATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGG  >  NC_000913/602591‑603081
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
gATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCTGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTg                                                                                                                                                                                                                                                  <  1:272277/251‑1 (MQ=255)
gATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCt                                                                                                                                                                                                                                                   <  1:404122/250‑1 (MQ=255)
     ttgcGCTGGCGGTAATTATGTTCCCCGTCGGCGGTCTGGAGCTGATGGGGAGTACTGCCGCTGAAGCGCGCGATCCGGAAAAAATCATTCCAAAAGCGGTAAATCAGGTGGGGTATCGCAGCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAACCTCGACAGCAACGTGGTAACTTCTGCGCt                                                                                                                                                                                                                                              <  2:647937/247‑1 (MQ=255)
      gTCGCGGGCGGTAATTATGTTCTCCCTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTg                                                                                                                                                                                                                                             <  2:503443/250‑1 (MQ=255)
          gtGGCGGTAATGATGTTCTCCTTCGGCGGGCTGGAGCTGATTGGGATTACTGCTGCTGAAGCGCGCGATCCGGAAAAAAGCATGCCAAAATCGGTAAAGCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGTTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCGTGACAGCAACGTGGTAACTTCTGCGCTGACCtt                                                                                                                                                                                                                                        <  2:579035/250‑1 (MQ=255)
            ggcgttatgtctGTTCTCCTTCGGGGGTCTGTAGCGTATTGTGATTGCTGTCGCTGAAGCGCGCGATTCGGAAAAAAGCATTCCAAAAGCGTTAAAGCAGGTCGTGTTTCGCATCCTTCTGTTTGACATCGGTTCCCTGGTGGTTGTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACATGAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTATCTTCTGCGTTGATCTTc                                                                                                                                                                                                                                       <  2:198860/239‑1 (MQ=255)
                   aTTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCt                                                                                                                                                                                                                               <  1:364589/251‑1 (MQ=255)
                     tATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTgg                                                                                                                                                                                                                             >  1:591718/1‑251 (MQ=255)
                      aTGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCGGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTgg                                                                                                                                                                                                                             <  2:406923/250‑1 (MQ=255)
                         ttCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGtcgc                                                                                                                                                                                                                         >  1:454526/1‑248 (MQ=255)
                             ccTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCg                                                                                                                                                                                                                     >  1:508177/1‑251 (MQ=255)
                               ttCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGATCACTGGGGGTTTTACTGGCGATCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGGGGAAACTTCTGCGCTGAAATTCGACATTCTGGTAACATAGCt                                                                                                                                                                                                                   >  2:37344/1‑251 (MQ=255)
                                                  ttGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGCCAGCACCGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCCTCGCTGTCAGTGtatacaagaggg                                                                                                                                                                                                >  2:173236/1‑243 (MQ=255)
                                                     ggATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGtt                                                                                                                                                                                             <  1:350309/251‑1 (MQ=255)
                                                     ggATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGtt                                                                                                                                                                                             >  1:396091/1‑251 (MQ=255)
                                                          aCTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTAct                                                                                                                                                                                         <  1:135640/250‑1 (MQ=255)
                                                           cgggcGCTGAAGCGCGCGAGCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATGCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTActct                                                                                                                                                                                       <  2:283328/247‑1 (MQ=255)
                                                             gcagCTGAAGCGCGAGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCGGCTGTTTTACATGGGTTCACTGGTGGGTTTACTGGCGCTCTAGCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGGCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTa                                                                                                                                                                                      <  2:222669/247‑1 (MQ=255)
                                                             gCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAAGCAGGTGGTGTATCGCATCCTGCTGGTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATATCAGCGGGGTTTActct                                                                                                                                                                                       <  2:396091/249‑1 (MQ=255)
                                                               cGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACa                                                                                                                                                                                   <  1:401336/251‑1 (MQ=255)
                                                                       cgcgcgATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATg                                                                                                                                                                           <  1:555455/251‑1 (MQ=255)
                                                                             aTCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGttt                                                                                                                                                                     >  1:536090/1‑251 (MQ=255)
                                                                                  gAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGATTACTCTAACAGCCGCAAGCTGTTTGGCCt                                                                                                                                                                >  2:205783/1‑251 (MQ=255)
                                                                                      aaaGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGGATGCTGTTTGGCCTTTCt                                                                                                                                                            >  2:343818/1‑251 (MQ=255)
                                                                                             tCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAgg                                                                                                                                                     >  1:666672/1‑251 (MQ=255)
                                                                                                       gTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCAATAATCTCGACAGCACCGTGGTACCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACCGCGGGGGTTACTCTAACAGCCGCATGCGGTTTGGCCTTTCTGTGCAGGGTAATgcgac                                                                                                                                           >  2:388643/1‑249 (MQ=255)
                                                                                                       gTAAATCAGGAGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGTCAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTACCAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGcc                                                                                                                                           >  2:68437/1‑251 (MQ=255)
                                                                                                               ggtggtGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGtttt                                                                                                                                   >  2:8526/1‑251 (MQ=255)
                                                                                                                     gtATCGCAGCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACt                                                                                                                              <  1:173236/250‑1 (MQ=255)
                                                                                                                          gCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGTTAATGCGCCGAAGTTTTTGACTCGCGTc                                                                                                                        >  1:305135/1‑251 (MQ=255)
                                                                                                                           cATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCTCGTTTGTGATGATTTTCCATACTCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCCTTCTGGCAACATCGCTGTCAGTGTATAACAGCGCGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGACGTTTTTGACTCCCGTcc                                                                                                                       >  1:77509/1‑250 (MQ=255)
                                                                                                                               ctgctgTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTTTGTGCAGGGTAATGCGCCGTAGGTTTTGACTCGCGTAAGCCg                                                                                                                   >  2:422868/1‑251 (MQ=255)
                                                                                                                                tgctgTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGt                                                                                                                  >  1:498614/1‑251 (MQ=255)
                                                                                                                                     ttttACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCgg                                                                                                             >  1:308071/1‑251 (MQ=255)
                                                                                                                                       ttACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATgcgc                                                                                                                                            <  1:503058/218‑1 (MQ=255)
                                                                                                                                       ttACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATgcgc                                                                                                                                            >  2:503058/1‑218 (MQ=255)
                                                                                                                                         acatcggtgctattgtggtTTTTCTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGGGAGGATTTTCCACAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCGTGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTGCTCTAACAGCCGCATGCGGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTcgttgtgt                                                                                                          <  2:657366/236‑7 (MQ=255)
                                                                                                                                              ggTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGAt                                                                                                    <  1:10146/251‑1 (MQ=255)
                                                                                                                                                      ggtggtTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGGGCCGATTAACTCGc                                                                                            >  1:312027/1‑251 (MQ=255)
                                                                                                                                                        tggtTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTg                                                                                          >  1:146151/1‑251 (MQ=255)
                                                                                                                                                         ggtTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTcg                                                                                                                <  1:629558/228‑1 (MQ=255)
                                                                                                                                                         ggtTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTcg                                                                                                                >  2:629558/1‑228 (MQ=255)
                                                                                                                                                             ttACTGGCGCTCTAGCCGTGGGTGGAAGTGAAATCCAACAGTATCCCGTTTGTGATGATTTTCCAAAATCTCGACAGCAACGTGGTAACTTCGGCGCTGAACGTCGTCAGTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCAGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGc                                                                                      <  1:301535/250‑1 (MQ=255)
                                                                                                                                                                  ggCGCTTTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTc                                                                                  <  1:446412/249‑1 (MQ=255)
                                                                                                                                                                        ctATCCGTGGGTGGAAGTGAAAGCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCTCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGc                                                                            <  1:4856/249‑1 (MQ=255)
                                                                                                                                                                                  gtggAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGTCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCt                                                                >  1:349470/1‑251 (MQ=255)
                                                                                                                                                                                      aaGTGAAAACCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGGCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCtggtg                                                            <  2:498614/251‑1 (MQ=255)
                                                                                                                                                                                               ccAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCt                                                                                                                                                                >  1:491812/1‑142 (MQ=255)
                                                                                                                                                                                               ccAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCt                                                                                                                                                                <  2:491812/142‑1 (MQ=255)
                                                                                                                                                                                               ccAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATc                                                   >  1:480781/1‑251 (MQ=255)
                                                                                                                                                                                                 aaCAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATAAATTCGCTGGAGGTGTTAATCa                                                  >  2:579653/1‑250 (MQ=255)
                                                                                                                                                                                                  aCAGGAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAAc                                                <  1:629594/251‑1 (MQ=255)
                                                                                                                                                                                                   cAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACt                                               <  1:8179/251‑1 (MQ=255)
                                                                                                                                                                                                       aGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATc                                            <  1:461405/250‑1 (MQ=255)
                                                                                                                                                                                                          ccGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATctgct                                        <  1:458104/251‑1 (MQ=255)
                                                                                                                                                                                                               tgtgATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCg                                    <  1:531508/250‑1 (MQ=255)
                                                                                                                                                                                                                gtgATGATTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGGTTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGcc                                  >  2:105022/1‑250 (MQ=255)
                                                                                                                                                                                                                     gcttttCCATCATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTCATCAACTATCTGCTGCCGCacaacg                             >  2:449196/3‑246 (MQ=255)
                                                                                                                                                                                                                      atTTTCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGATGCCGCAAAAAGc                            >  1:137894/1‑251 (MQ=255)
                                                                                                                                                                                                                        tttCCATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGGGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGt                          >  2:19526/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccATAATCTCCACAGCAACGTGGTAACTTCTGCGCTGAACTTCGCCATTCTGGTATCATCGCTGTCAGTGTATAATAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGGGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACGTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTg                       <  2:552742/251‑1 (MQ=255)
                                                                                                                                                                                                                            cATAATCTCGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTg                       >  1:552742/1‑250 (MQ=255)
                                                                                                                                                                                                                                 tctcGACAGCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGGGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGGGGTGTTAATCACCTATCTGCTGCAGCCAGTAGCGTTTGGGctgc                 >  2:9066/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ctcGACAGCAACGTGGTAACTTATGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTAATCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATAACTCCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTctgct                >  2:558444/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aCAGCAACGTGGTAACTGCTGCCCTGAACTTCTTCATTCTGGTAGCATCCCTGTCAGTGTATAACGGCTGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTGTGACTCGCGTCAGCCGTCGCGGTTTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATg            <  2:70610/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCAACGTTGTAACTTCTGCGCTTAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGTGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGgcg         <  1:356448/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCAACGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGAt                                                                         >  1:615819/1‑187 (MQ=255)
                                                                                                                                                                                                                                           aaCGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCt       <  1:8526/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aaCGTGGTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGCGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGAGTTTGGTCTGCTGATGgcgcg       >  2:143054/1‑250 (MQ=255)
                                                                                                                                                                                                                                                ggTAACTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGATTGGTCTGCTGAGGGCGCtggggg  >  2:25396/1‑248 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
GATTTTGTCGCTGGCGGTAATTATGTTCTCCTTCGGCGGTCTGGAGCTGATTGGGATTACTGCCGCTGAAGCGCGCGATCCGGAAAAAAGCATTCCAAAAGCGGTAAATCAGGTGGTGTATCGCATCCTGCTGTTTTACATCGGTTCACTGGTGGTTTTACTGGCGCTCTATCCGTGGGTGGAAGTGAAATCCAACAGTAGCCCGTTTGTGATGATTTTCCATAATCTCGACAGCAACGTGGTAGCTTCTGCGCTGAACTTCGTCATTCTGGTAGCATCGCTGTCAGTGTATAACAGCGGGGTTTACTCTAACAGCCGCATGCTGTTTGGCCTTTCTGTGCAGGGTAATGCGCCGAAGTTTTTGACTCGCGTCAGCCGTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGG  >  NC_000913/602591‑603081

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