Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,987 G→C 6.9% E123Q (GAG→CAG)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,9870GC 6.9% 197.5 / 6.5 58E123Q (GAG→CAG) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base C (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGCCAATGAGT  >  NC_000913/1147700‑1148278
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                   
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)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                   
CACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGAAATTCCCAATCCTCGCGTGGCGTTGCTCAATATTGGTGAAGAAGAAGTAAAGGGTCTCGACAGTATTCGGGATGCCTCAGCGGTGCTTAAAACAATCCCTTCTATCAATTATATCGGCTATCTTGAAGCCAATGAGT  >  NC_000913/1147700‑1148278

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: ‑