Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,006,954 T→A 13.3% L181Q (CTG→CAG)  metE → 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,006,9540TA13.3% 1105.1 / 140.4 354L181Q (CTG→CAG) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (169/138);  new base A (21/26);  total (190/164)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 360 total aligned reads displayed.

TAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACG  >  CP009273/4006805‑4007100
                                                                                                                                                     |                                                                                                                                                  
 aaaGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTg                                                                                                                                                   <  1:2844977/150‑1 (MQ=255)
 aaaGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTg                                                                                                                                                   >  1:2469329/1‑150 (MQ=255)
   aGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCagag                                                                                                                                                 >  1:694487/1‑150 (MQ=255)
       cAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCctg                                                                                                                                             <  1:2207325/150‑1 (MQ=255)
           aGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGa                                                                                                                                         <  1:4387843/150‑1 (MQ=255)
                  aCTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACAt                                                                                                                                  <  1:9195/150‑1 (MQ=255)
                     gACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCt                                                                                                                               <  1:3790467/150‑1 (MQ=255)
                      aCCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTg                                                                                                                              >  2:2691876/1‑150 (MQ=255)
                        cTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGcc                                                                                                                            <  2:4196678/150‑1 (MQ=255)
                        cTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCTGcc                                                                                                                            >  2:2502995/1‑150 (MQ=255)
                         tGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCg                                                                                                                           >  2:3340151/1‑150 (MQ=255)
                          ggACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCgg                                                                                                                          <  1:1004951/150‑1 (MQ=255)
                           gACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGt                                                                                                                         <  1:1509434/150‑1 (MQ=255)
                            aCGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGtt                                                                                                                        >  2:2456289/1‑150 (MQ=255)
                            aCGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGtt                                                                                                                        >  1:625955/1‑150 (MQ=255)
                            aCGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGtt                                                                                                                        <  2:1093205/150‑1 (MQ=255)
                             cGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGttt                                                                                                                       >  1:2697068/1‑150 (MQ=255)
                               cAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTAt                                                                                                                     <  2:1354424/150‑1 (MQ=255)
                               cAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCTGCCGGTTTAt                                                                                                                     <  1:569376/150‑1 (MQ=255)
                                 gctgctGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCa                                                                                                                   >  1:712809/1‑150 (MQ=255)
                                  ctgctgGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAg                                                                                                                  <  1:1706649/150‑1 (MQ=255)
                                      tgGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAg                                                                                                              >  1:4640632/1‑150 (MQ=255)
                                           gAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTg                                                                                                         <  1:3568641/150‑1 (MQ=255)
                                             aGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGTTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGc                                                                                                       >  1:3472062/1‑150 (MQ=255)
                                                      ggcgctggcgctGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGc                                                                                              >  2:141123/1‑150 (MQ=255)
                                                         gctggcgctggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGTGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGaa                                                                                           <  2:4032325/150‑1 (MQ=255)
                                                           tggcgctggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAAc                                                                                         >  1:1878780/1‑150 (MQ=255)
                                                           tggcgctggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAAc                                                                                         <  2:1013998/150‑1 (MQ=255)
                                                            ggcgctggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAAcg                                                                                        <  2:1497318/150‑1 (MQ=255)
                                                              cgctggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAAcgcg                                                                                      >  1:3239399/1‑150 (MQ=255)
                                                              cgctggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAAcgcg                                                                                      <  2:2271616/150‑1 (MQ=255)
                                                                ctggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGc                                                                                    >  2:2888227/1‑150 (MQ=255)
                                                                 tggGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCa                                                                                   <  2:817051/150‑1 (MQ=255)
                                                                  gggCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCAt                                                                                  <  2:2342186/150‑1 (MQ=255)
                                                                   ggCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATc                                                                                 >  2:2585849/1‑150 (MQ=255)
                                                                      cacaAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAg                                                                              <  1:3008829/150‑1 (MQ=255)
                                                                        caAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTg                                                                            >  1:954258/1‑150 (MQ=255)
                                                                        caAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTg                                                                            >  2:4830808/1‑150 (MQ=255)
                                                                        caAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTg                                                                            <  1:4550159/150‑1 (MQ=255)
                                                                          aGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGt                                                                         >  2:3198721/1‑150 (MQ=255)
                                                                             tGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTAc                                                                       >  2:3898904/1‑150 (MQ=255)
                                                                                aaCCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGa                                                                    >  2:4596359/1‑150 (MQ=255)
                                                                                aaCCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGa                                                                    <  2:4543660/150‑1 (MQ=255)
                                                                                 aCCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGAt                                                                   <  1:405985/150‑1 (MQ=255)
                                                                                  ccTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGAtt                                                                  <  1:2632183/150‑1 (MQ=255)
                                                                                   cTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATtg                                                                 <  2:3948367/150‑1 (MQ=255)
                                                                                   cTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATtg                                                                 <  2:584259/150‑1 (MQ=255)
                                                                                        ctgctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGaa                                                            >  1:470119/1‑150 (MQ=255)
                                                                                        ctgctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGaa                                                            <  1:3156695/150‑1 (MQ=255)
                                                                                        ctgctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGaa                                                            <  1:498278/150‑1 (MQ=255)
                                                                                         tgctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAAc                                                           >  1:1086444/1‑150 (MQ=255)
                                                                                         tgctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAAc                                                           >  2:2702567/1‑150 (MQ=255)
                                                                                            tgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCcg                                                        <  2:4691213/150‑1 (MQ=255)
                                                                                            tgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCcg                                                        >  2:445083/1‑150 (MQ=255)
                                                                                            tgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCcg                                                        >  1:1310402/1‑150 (MQ=255)
                                                                                            tgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCcg                                                        <  2:3672863/150‑1 (MQ=255)
                                                                                             ggggCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCcgc                                                       >  2:2415790/1‑150 (MQ=255)
                                                                                                  cGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTgg                                                  >  1:4857598/1‑150 (MQ=255)
                                                                                                   ggTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGt                                                 >  1:2556618/1‑150 (MQ=255)
                                                                                                    gTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTa                                                <  2:1927457/150‑1 (MQ=255)
                                                                                                     ttACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTAc                                               >  2:301713/1‑150 (MQ=255)
                                                                                                       aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTg                                             >  2:2054636/1‑150 (MQ=255)
                                                                                                       aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTg                                             >  1:3084063/1‑150 (MQ=255)
                                                                                                        ccTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTgg                                            <  1:2934871/150‑1 (MQ=255)
                                                                                                        ccTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTgg                                            >  2:370720/1‑150 (MQ=255)
                                                                                                         cTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGa                                           >  2:787365/1‑150 (MQ=255)
                                                                                                          tggctgtggctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGaa                                          >  1:233616/1‑150 (MQ=255)
                                                                                                            gctgtggctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACt                                        <  2:1280413/150‑1 (MQ=255)
                                                                                                            gctgtggctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACt                                        >  2:1237821/1‑150 (MQ=255)
                                                                                                             ctgtggctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTa                                       <  2:655441/150‑1 (MQ=255)
                                                                                                              tgtggctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTAc                                      >  2:2658563/1‑150 (MQ=255)
                                                                                                              tgtggctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTAc                                      >  2:1862531/1‑150 (MQ=255)
                                                                                                                tggctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACca                                    <  1:1767373/150‑1 (MQ=255)
                                                                                                                  gctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACcaca                                  <  2:2929818/150‑1 (MQ=255)
                                                                                                                  gctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACcaca                                  >  2:4376288/1‑150 (MQ=255)
                                                                                                                  gctgGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACcaca                                  <  2:3085277/150‑1 (MQ=255)
                                                                                                                        gAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTg                            >  2:4140484/1‑150 (MQ=255)
                                                                                                                        gAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACACGCGCTGGTACTGGAACTACCACAGGCGTg                            >  1:521175/1‑150 (MQ=255)
                                                                                                                         aaaGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTgg                           >  2:773586/1‑150 (MQ=255)
                                                                                                                          aaGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGc                          >  1:3170310/1‑150 (MQ=255)
                                                                                                                          aaGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGc                          >  1:3710401/1‑150 (MQ=255)
                                                                                                                          aaGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGc                          >  1:738998/1‑150 (MQ=255)
                                                                                                                           aGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCt                         <  2:2741069/150‑1 (MQ=255)
                                                                                                                           aGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCt                         >  2:1336876/1‑150 (MQ=255)
                                                                                                                            gTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTg                        <  1:1242299/150‑1 (MQ=255)
                                                                                                                             tGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTgg                       >  2:1688953/1‑150 (MQ=255)
                                                                                                                              gAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGa                      <  2:502157/150‑1 (MQ=255)
                                                                                                                              gAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGa                      <  2:3807614/150‑1 (MQ=255)
                                                                                                                              gAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGa                      >  2:4261505/1‑150 (MQ=255)
                                                                                                                              gAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGa                      >  2:264289/1‑150 (MQ=255)
                                                                                                                                aaGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACg                    >  1:3369049/1‑150 (MQ=255)
                                                                                                                                aaGGTGAACAATTTGACCGCCAGAGCCTGCTGAACGACATTCGGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACg                    <  2:1763540/150‑1 (MQ=255)
                                                                                                                                    tGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATa                <  2:354637/150‑1 (MQ=255)
                                                                                                                                    tGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATa                >  2:2172037/1‑150 (MQ=255)
                                                                                                                                         aaTTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAAc           <  2:3599330/150‑1 (MQ=255)
                                                                                                                                          aTTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAAcc          >  1:2028347/1‑150 (MQ=255)
                                                                                                                                                ccGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTTCTGGCAGATCTGGCGAAACGCGGCATCGAGTGTGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCTTGGCTGGACGCATACAAACCCGCTTa    <  1:4184285/150‑1 (MQ=255)
                                                                                                                                                ccGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTa    <  2:2137551/150‑1 (MQ=255)
                                                                                                                                                 cGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTac   >  2:698168/1‑150 (MQ=255)
                                                                                                                                                 cGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTac   <  2:1634448/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGACGAGGCGCTGGCGCTGGGCCACAAGGTGAAACCTGTGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAATTTGACCGCCTGAGCCTGCTGAACGACATTCTGCCGGTTTATCAGCAAGTGCTGGCAGAACTGGCGAAACGCGGCATCGAGTGGGTACAGATTGATGAACCCGCGCTGGTACTGGAACTACCACAGGCGTGGCTGGACGCATACAAACCCGCTTACG  >  CP009273/4006805‑4007100

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: