Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,993 A→C 6.8% I125L (ATT→CTT)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,9930AC 6.8% 179.7 / 3.6 58I125L (ATT→CTT) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base C (2/2);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

CACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGCCAATGAGTTGTTAAC  >  NC_000913/1147700‑1148285
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                    
cACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                                                                                                                                                                                  <  1:597770/298‑1 (MQ=255)
 aCTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGGTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAg                                                                                                                                                                                                                                                                                                 <  1:368488/298‑1 (MQ=255)
      aTTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCg                                                                                                                                                                                                                                                                                           <  1:298367/299‑1 (MQ=255)
         cTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                                                                                                                                                                      <  2:345423/285‑1 (MQ=255)
          tAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGgcgc                                                                                                                                                                                                                                                                                       >  1:5598/1‑299 (MQ=255)
          tAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCAGCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                                                                                                                                                                      <  2:106778/284‑1 (MQ=255)
               cGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAg                                                                                                                                                                                                                                                                                                 <  2:372742/284‑1 (MQ=255)
                 cAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAc                                                                                                                                                                                                                                                                                >  1:415557/1‑299 (MQ=255)
                    ctcactcTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCg                                                                                                                                                                                                                                                                                           <  2:501163/285‑1 (MQ=255)
                       actcTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTAAAGCCCCTGGAGGGGATTGAGCGTCCGgcg                                                                                                                                                                                                                                                                                        >  2:79968/1‑285 (MQ=255)
                         tcttctTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGgcgc                                                                                                                                                                                                                                                                                       >  2:236990/1‑284 (MQ=255)
                            tcttTTAGTCGGCAATTCCGGCGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACcac                                                                                                                                                                                                                                                                     <  1:201365/299‑1 (MQ=255)
                                       gCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTTGGAGTTCACTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAg                                                                                                                                                                                                                                                          >  1:268206/1‑299 (MQ=255)
                                          aTTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAACGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGc                                                                                                                                                                                                                                                       >  1:365981/1‑299 (MQ=255)
                                                            tcttccgatcTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                                                                                                                                                                                      <  1:226121/289‑1 (MQ=255)
                                                                        aaGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCCGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGACGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACAGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAACGGGCAAAACGGTGGTCCtt                                                                                                                                                                                                                                        >  2:198895/1‑284 (MQ=255)
                                                                                     acaacGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGATGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACCCAAGCCCCTGGAGGGGATTGAGCCTACGGAGCTGGTGACGGTATTACCACATCAGCAAAAGAGAAAAACGGGGGTCCTTGACTTATGGGCCa                                                                                                                                                                                                                           >  2:85331/1‑284 (MQ=255)
                                                                                     acaacGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGAc                                                                                                                                                                                                             <  1:530855/298‑1 (MQ=255)
                                                                                     acaacGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGAc                                                                                                                                                                                                             <  1:521212/298‑1 (MQ=255)
                                                                                                 tCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGGCGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTg                                                                                                                                                                                                 <  1:304379/298‑1 (MQ=255)
                                                                                                     cAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGc                                                                                                                                                                                                          <  2:304336/285‑1 (MQ=255)
                                                                                                                 gcgcAGTCAGTTATCGCCAGTGATGCCCGGCCTGCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                                                                                                                                                                               <  2:119609/284‑1 (MQ=255)
                                                                                                                   gcAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                                                                                                                                                                              >  1:116273/1‑299 (MQ=255)
                                                                                                                     agtcagtTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATggg                                                                                                                                                                             <  1:216675/298‑1 (MQ=255)
                                                                                                                                           ccGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAggtgg                                                                                                                                                      >  1:404805/1‑299 (MQ=255)
                                                                                                                                                      aaGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTccc                                                                                                                                           >  1:350001/1‑299 (MQ=255)
                                                                                                                                                      aaGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTccc                                                                                                                                           >  1:467400/1‑299 (MQ=255)
                                                                                                                                                        gCTAGCCGCGCCAGTCGTGGGAGGTCACTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGTTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCGTGACTTATGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAtgt                                                                                                                                                        <  2:86300/284‑1 (MQ=255)
                                                                                                                                                               gcgcCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCACCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTAGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGATCAGTTCTGGCTGCAGAGGTGGTGGaa                                                                                                                                                 >  2:275433/1‑284 (MQ=255)
                                                                                                                                                                     gTCGTGGGAGTTCAATGCGCGTGGCCCGGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCa                                                                                                                                          <  2:415665/285‑1 (MQ=255)
                                                                                                                                                                            gAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCa                                                                                                                     >  1:271424/1‑299 (MQ=255)
                                                                                                                                                                             aGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTcgc                                                                                                                                  >  2:315424/1‑285 (MQ=255)
                                                                                                                                                                             aGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCaa                                                                                                                    >  1:554959/1‑299 (MQ=255)
                                                                                                                                                                                       gcgTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGTCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTga                                                                                                           <  1:97606/298‑1 (MQ=255)
                                                                                                                                                                                         gTGGCCCTGGGGCGGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCt                                                                                                                       <  2:408048/284‑1 (MQ=255)
                                                                                                                                                                                            ctcttccgctcttggAAAGAAGTGCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAat                                                                                                                   <  2:62635/271‑1 (MQ=255)
                                                                                                                                                                                              ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGAGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGaagaagaag                                                                                                   >  1:134767/1‑299 (MQ=255)
                                                                                                                                                                                               cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGGCAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGt                                                                                                  >  1:362886/1‑299 (MQ=255)
                                                                                                                                                                                                      tGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGCAGTGGAAATTCCCAATACTCGCGTGGCGTTGCTCAATATTGGTgaa                                                                                                          >  2:478676/1‑284 (MQ=255)
                                                                                                                                                                                                             agaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACa                                                                                    <  1:187337/299‑1 (MQ=255)
                                                                                                                                                                                                                agGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGGGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTa                                                                                 <  1:312275/299‑1 (MQ=255)
                                                                                                                                                                                                                     gAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACACTGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAACTTCCCAATCCTCGCGTGGCGTTGCCCACTATTGGTGAAGAAGAAGTAAAGGGtc                                                                                          >  2:3706/1‑285 (MQ=255)
                                                                                                                                                                                                                              ccTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCa                                                                   >  1:85220/1‑299 (MQ=255)
                                                                                                                                                                                                                                  tgtCAGTGCCGGTAATACCGGGGCGCTGATTGGGCGGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCg                                                                             <  2:354262/285‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTg                                                             <  1:268890/298‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTa                                                         >  1:325559/1‑299 (MQ=255)
                                                                                                                                                                                                                                                 tACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATcc                                                <  1:91647/299‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATc                                                 >  2:433924/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                         aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTa                          >  1:269815/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                           aaTTATTACTCAAGCACCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTAGTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGCGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAAtt                                    <  2:236365/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ttaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAg                  >  1:312390/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                tgaCTCAAGCCCCTGGAGGGGATGGCGCGTCCGGCGCTGGTGACGGTATTACCACATCATCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTatat                                <  2:286210/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 cgcTCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGcc                >  1:393371/3‑299 (MQ=255)
                                                                                                                                                                                                                                                                                  gcTCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTAAGGGCTCAGTTCTGGCGGAAGAGGTGGTGGAGATTCCCAATCCTCGCGTGGCGTGGCTAAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACCGTCTTCGTGATGCCTAAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCgg                             >  2:379320/2‑285 (MQ=255)
                                                                                                                                                                                                                                                                                    tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAACCAATCCCTTCTATCAATTATCTCGGCt                           >  2:210284/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                      aaGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGCCAATGa           >  1:217101/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                        gacccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGACAATGAgt         >  1:31234/3‑299 (MQ=255)
                                                                                                                                                                                                                                                                                         ccccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCtt                      >  2:470303/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                              ggAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAg                  >  2:435662/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                                                aGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTTGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGTGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGCCAATGAGTTGTTAAc  <  1:386703/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                    
CACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGCCAATGAGTTGTTAAC  >  NC_000913/1147700‑1148285

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