Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,991,234 T→A 6.5% L231* (TTG→TAG)  yigB → 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,991,2340TA 6.5% 916.5 / 42.4 263L231* (TTG→TAG) yigB5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis
Reads supporting (aligned to +/- strand):  ref base T (106/139);  new base A (5/12);  total (111/152)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 264 total aligned reads displayed.

AGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAA  >  CP009273/3991087‑3991381
                                                                                                                                                   |                                                                                                                                                   
agaTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTg                                                                                                                                                    >  2:4091574/1‑149 (MQ=255)
 gaTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCCTCTGGAAATTTCCCGGTTgg                                                                                                                                                   >  2:2367102/1‑149 (MQ=255)
  aTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGc                                                                                                                                                  >  1:4465254/1‑149 (MQ=255)
  aTCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGc                                                                                                                                                  >  2:2094913/1‑149 (MQ=255)
    cTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAt                                                                                                                                                >  1:4387851/1‑149 (MQ=255)
       acttGTTGGGGACGATCTCACCTCTGACGTGGGTGGGGCAATTCGCAGCGGAATGCTGGCTTGTTGGATCTGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAtctc                                                                                                                                             <  1:3191854/146‑1 (MQ=255)
       aCATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAtctc                                                                                                                                             <  1:1182313/149‑1 (MQ=255)
       aCATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCAtctc                                                                                                                                             <  2:3138560/149‑1 (MQ=255)
        cATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATctct                                                                                                                                            <  1:3022533/149‑1 (MQ=255)
        cATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATctct                                                                                                                                            <  1:474835/149‑1 (MQ=255)
        cATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATctct                                                                                                                                            <  2:2179105/149‑1 (MQ=255)
         aTGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTg                                                                                                                                           >  2:909349/1‑149 (MQ=255)
         aTGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTg                                                                                                                                           >  1:1662381/1‑149 (MQ=255)
          tGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGa                                                                                                                                          <  2:2192091/149‑1 (MQ=255)
          tGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGa                                                                                                                                          >  1:1980828/1‑149 (MQ=255)
                ggACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGc                                                                                                                                    <  2:2765134/149‑1 (MQ=255)
                ggACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGc                                                                                                                                    <  2:1114913/149‑1 (MQ=255)
                 gACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCt                                                                                                                                   >  2:2018787/1‑149 (MQ=255)
                  aCGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTg                                                                                                                                  <  1:1421922/149‑1 (MQ=255)
                     aTCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGata                                                                                                                               <  2:4406364/149‑1 (MQ=255)
                     aTCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGata                                                                                                                               >  2:551513/1‑149 (MQ=255)
                     aTCTCACCACTGACGTGGGTGGGGCAAGTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGata                                                                                                                               <  2:2175556/149‑1 (MQ=255)
                       ctcACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGAtata                                                                                                                             >  1:754968/1‑149 (MQ=255)
                         caccacTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAAt                                                                                                                           <  1:504573/149‑1 (MQ=255)
                          accacTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATc                                                                                                                          <  2:2788128/149‑1 (MQ=255)
                          accacTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATc                                                                                                                          >  2:716363/1‑149 (MQ=255)
                           ccacTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCa                                                                                                                         >  2:1894653/1‑149 (MQ=255)
                                 ccGTGGGGGGGGCAAGGCGCAGCGGCATGCAGGCGTGTTGGATCAGACCGGAAACTGGCGATCTGATGCAAACCTGGGACAGCCGTTTACGGCCGCATCTGGAAATTTCCCGGTTGGCATCTCGGACCTCGCTGATATAATCAGCAAAt                                                                                                                   <  2:3962653/148‑1 (MQ=255)
                                 aCGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAAt                                                                                                                   <  1:1117767/149‑1 (MQ=255)
                                 aCGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAAt                                                                                                                   <  1:4249789/149‑1 (MQ=255)
                                  cGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATc                                                                                                                  <  1:3698734/149‑1 (MQ=255)
                                    tgggtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTg                                                                                                                <  2:3736165/149‑1 (MQ=255)
                                     gggtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATACTCAGCCAATCTGt                                                                                                               >  2:3997291/1‑149 (MQ=255)
                                       gtgggGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGtat                                                                                                             <  2:1499206/149‑1 (MQ=255)
                                          gggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGtatata                                                                                                          >  1:2027511/1‑149 (MQ=255)
                                          gggCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGtatata                                                                                                          >  2:3488437/1‑149 (MQ=255)
                                              aaTTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATAcc                                                                                                      <  2:4238571/149‑1 (MQ=255)
                                               aTTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGTAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATAccc                                                                                                     >  1:4014752/1‑149 (MQ=255)
                                                ttCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCa                                                                                                    <  1:121292/149‑1 (MQ=255)
                                                ttCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCa                                                                                                    <  1:3392828/149‑1 (MQ=255)
                                                ttCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCa                                                                                                    >  1:4667587/1‑149 (MQ=255)
                                                ttCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCa                                                                                                    >  1:2450589/1‑149 (MQ=255)
                                                   gCAGCTGAATGCAGGCTTGTTGGATCAGACCGGAATTTGGCGATCTGTTGCAATCCTGGGTCAGCCGTTTACTGCCGCTTCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCt                                                                                                 <  1:2672760/149‑1 (MQ=255)
                                                   gCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCt                                                                                                 <  1:3447160/149‑1 (MQ=255)
                                                     aGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCttt                                                                                               >  1:3986003/1‑149 (MQ=255)
                                                     aGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCttt                                                                                               <  2:1004450/149‑1 (MQ=255)
                                                               aGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAgg                                                                                     >  1:1748868/1‑149 (MQ=255)
                                                                 gCTTGTTGGATCAGACCGGAAACTGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGc                                                                                   <  1:459835/149‑1 (MQ=255)
                                                                   ttgttgGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGt                                                                                 >  2:4476080/1‑149 (MQ=255)
                                                                     gttgGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTg                                                                               <  1:1868330/149‑1 (MQ=255)
                                                                      ttgGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTgc                                                                              <  2:1317046/149‑1 (MQ=255)
                                                                      ttgGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTgc                                                                              >  2:3062545/1‑149 (MQ=255)
                                                                        ggATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTgcgc                                                                            >  1:3196592/1‑149 (MQ=255)
                                                                        ggATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTgcgc                                                                            >  1:4559548/1‑149 (MQ=255)
                                                                         gATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCt                                                                           <  2:3210636/149‑1 (MQ=255)
                                                                          aTCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCtt                                                                          >  2:3580985/1‑149 (MQ=255)
                                                                                     aaaTGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTCTATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAAc                                                               >  2:1687194/1‑149 (MQ=255)
                                                                                      aaTGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAAcc                                                              >  2:4538955/1‑149 (MQ=255)
                                                                                      aaTGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAAcc                                                              >  1:3726950/1‑149 (MQ=255)
                                                                                         ggCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGGTGGCATCTCTGACCTCGCTGATATACTCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAt                                                           <  1:4374554/149‑1 (MQ=255)
                                                                                           cGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAttt                                                         <  1:2166893/149‑1 (MQ=255)
                                                                                           cGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAttt                                                         <  1:1053180/149‑1 (MQ=255)
                                                                                            gATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTAtttt                                                        >  1:910228/1‑149 (MQ=255)
                                                                                                tGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTAcg                                                    >  2:1793609/1‑149 (MQ=255)
                                                                                                    gCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACgcggc                                                <  1:2989689/149‑1 (MQ=255)
                                                                                                    gCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACgcggc                                                <  1:2571546/149‑1 (MQ=255)
                                                                                                      aaaCCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGcggcgg                                              >  2:4305736/1‑149 (MQ=255)
                                                                                                      aaaCCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGcggcgg                                              <  1:1300894/149‑1 (MQ=255)
                                                                                                      aaaCCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGAGATAATCAGCAAATCGGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGcggcgg                                              <  2:2604649/149‑1 (MQ=255)
                                                                                                      aaaCCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGcggcgg                                              <  1:829192/149‑1 (MQ=255)
                                                                                                          cTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGcc                                          <  1:2259729/149‑1 (MQ=255)
                                                                                                          cTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGcc                                          >  2:4285024/1‑149 (MQ=255)
                                                                                                           tGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCa                                         >  1:354448/1‑149 (MQ=255)
                                                                                                            gggACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCaa                                        >  1:2053216/1‑149 (MQ=255)
                                                                                                                cAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGa                                    <  2:3927680/149‑1 (MQ=255)
                                                                                                                  gCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACg                                  >  1:286892/1‑149 (MQ=255)
                                                                                                                    cGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGtt                                >  1:1493560/1‑149 (MQ=255)
                                                                                                                      tttACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTc                              <  2:2873281/149‑1 (MQ=255)
                                                                                                                      tttACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTc                              >  2:2492661/1‑149 (MQ=255)
                                                                                                                       ttACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCt                             <  1:2759303/149‑1 (MQ=255)
                                                                                                                       ttACTGCCGCATCTGGAAATGTCCCGGGTGGCATCTCTGACCTCGCTGATATACTCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCt                             <  2:4306695/149‑1 (MQ=255)
                                                                                                                        tACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCtt                            <  2:4431458/149‑1 (MQ=255)
                                                                                                                         aCTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTa                           >  1:2931825/1‑149 (MQ=255)
                                                                                                                          cTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTAc                          <  1:2404247/149‑1 (MQ=255)
                                                                                                                          cTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTAc                          <  1:62799/149‑1 (MQ=255)
                                                                                                                            gCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCt                        >  2:3322733/1‑149 (MQ=255)
                                                                                                                              cGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGc                      <  2:3059368/149‑1 (MQ=255)
                                                                                                                                 ctCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCg                   >  1:1533382/2‑149 (MQ=255)
                                                                                                                                 aTCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCCACCTCTTTGTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCg                   >  2:3833178/1‑149 (MQ=255)
                                                                                                                                  tCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGa                  <  2:502751/149‑1 (MQ=255)
                                                                                                                                     ggAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAg               <  1:4343605/149‑1 (MQ=255)
                                                                                                                                        aaggtCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCt            <  1:2190641/145‑1 (MQ=255)
                                                                                                                                         aTTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCACCCTATTGTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCtt           >  2:3454155/1‑149 (MQ=255)
                                                                                                                                           ttCCCGGTGGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTaa         <  1:2605129/149‑1 (MQ=255)
                                                                                                                                             cccGGTAGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTGTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATg       <  2:4667587/149‑1 (MQ=255)
                                                                                                                                              ccGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGa      >  2:41993/1‑149 (MQ=255)
                                                                                                                                              ccGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTACTTACCTGCTCGACAGCCTTAATGa      >  2:2038067/1‑149 (MQ=255)
                                                                                                                                               cGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGAc     >  1:4425321/1‑149 (MQ=255)
                                                                                                                                               cGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGAc     <  1:3254225/149‑1 (MQ=255)
                                                                                                                                                  ttGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACaaa  <  1:3782774/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGATCTTACATGTTGGGGACGATCTCACCACTGACGTGGGTGGGGCAATTCGCAGCGGAATGCAGGCTTGTTGGATCAGACCGGAAAATGGCGATCTGATGCAAACCTGGGACAGCCGTTTACTGCCGCATCTGGAAATTTCCCGGTTGGCATCTCTGACCTCGCTGATATAATCAGCAAATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAA  >  CP009273/3991087‑3991381

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