Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,135,168 G→T 100% V416L (GTA→TTA)  pflD → putative glycine radical domain‑containing pyruvate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,135,1680GT98.0% 631.3 / 3.6 200V416L (GTA→TTA) pflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base T (88/108);  total (88/112)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 200 total aligned reads displayed.

CGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCGTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGTTTACTGGAGC  >  CP009273/4135023‑4135312
                                                                                                                                                 |                                                                                                                                                
cGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTa                                                                                                                                                >  1:698293/1‑148 (MQ=255)
 gCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAg                                                                                                                                               <  1:3707289/148‑1 (MQ=255)
     tGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTggg                                                                                                                                           >  1:3008708/1‑148 (MQ=255)
     tGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTggg                                                                                                                                           >  1:3008715/1‑148 (MQ=255)
     tGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTggg                                                                                                                                           <  1:1704785/148‑1 (MQ=255)
        tCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCGTAGTGGGCtg                                                                                                                                        <  2:3757487/148‑1 (MQ=255)
         cGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCtgt                                                                                                                                       >  1:2951210/1‑148 (MQ=255)
          gACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCtgtg                                                                                                                                      <  1:2044616/148‑1 (MQ=255)
          gACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCtgtg                                                                                                                                      <  1:2018178/148‑1 (MQ=255)
           acacGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTAGTCCTTAGTGGGCtgtgt                                                                                                                                     <  2:2930600/148‑1 (MQ=255)
           acacGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCCTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCtgtgt                                                                                                                                     >  2:1294106/1‑148 (MQ=255)
                  gTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGACTTa                                                                                                                              >  1:446546/1‑148 (MQ=255)
                  gTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTa                                                                                                                              >  2:1651383/1‑148 (MQ=255)
                   ttCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTAt                                                                                                                             <  2:340815/148‑1 (MQ=255)
                   ttCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTAt                                                                                                                             <  1:3139263/148‑1 (MQ=255)
                   ttCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTAt                                                                                                                             >  1:2897150/1‑148 (MQ=255)
                     ccTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCt                                                                                                                           <  2:2681020/148‑1 (MQ=255)
                        gatgaAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTAtt                                                                                                                        >  1:1624510/1‑148 (MQ=255)
                        gatgaAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTAtt                                                                                                                        >  1:1474836/1‑148 (MQ=255)
                          tgaAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTcc                                                                                                                      <  2:1251528/148‑1 (MQ=255)
                          tgaAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTcc                                                                                                                      >  1:957222/1‑148 (MQ=255)
                            aaaaCCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGATGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCg                                                                                                                    >  2:1628813/1‑148 (MQ=255)
                              aaCCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAAGTATCTATTCCCGGc                                                                                                                  >  2:844806/1‑148 (MQ=255)
                               aCCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCa                                                                                                                 >  2:1593050/1‑148 (MQ=255)
                                ccgccgAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAg                                                                                                                >  2:2042322/1‑148 (MQ=255)
                                 cgccgAAACCATTCGCTTCGGTACCGGTATTCCGCAAATATTTAACGATGAAGTGGTGGTGCCAGCGTTACTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGa                                                                                                               <  1:3658299/148‑1 (MQ=255)
                                      aaaCCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTa                                                                                                          <  1:3638397/148‑1 (MQ=255)
                                      aaaCCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTa                                                                                                          <  1:2020157/148‑1 (MQ=255)
                                       aaCCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTAc                                                                                                         <  1:1720289/148‑1 (MQ=255)
                                         ccATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACgg                                                                                                       >  2:1347820/1‑148 (MQ=255)
                                          cATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGc                                                                                                      <  1:2429047/148‑1 (MQ=255)
                                          cATTCGCTTCGATACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGc                                                                                                      <  1:2583643/148‑1 (MQ=255)
                                              cGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGc                                                                                                  <  2:1451615/148‑1 (MQ=255)
                                              cGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGc                                                                                                  <  2:1451247/148‑1 (MQ=255)
                                              cGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTGTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAAGTAGCTATTCCCGGCAGAACCTACGGCGTGc                                                                                                  <  2:3202563/148‑1 (MQ=255)
                                               gCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGAGGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCa                                                                                                 <  1:2271545/148‑1 (MQ=255)
                                                 ttCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATg                                                                                               <  2:3538597/148‑1 (MQ=255)
                                                  tCGGTACCGGTATTCCGCATAACTTTAACGATGAAGTTGTGGTGCCAGCGTTCCTCTACCGTGGCGTTTCGCTGGAAGAGGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTTGGGCTTGCATGa                                                                                              <  2:2852011/148‑1 (MQ=255)
                                                  tCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGa                                                                                              <  1:1518743/148‑1 (MQ=255)
                                                  tCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGa                                                                                              >  1:2476882/1‑148 (MQ=255)
                                                   cGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGAc                                                                                             <  2:2350305/148‑1 (MQ=255)
                                                    ggTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACa                                                                                            <  2:3595455/148‑1 (MQ=255)
                                                    ggTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACa                                                                                            <  2:956822/148‑1 (MQ=255)
                                                        ccGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATcgc                                                                                        >  2:3252081/1‑148 (MQ=255)
                                                         cGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATcgcg                                                                                       <  2:2957873/148‑1 (MQ=255)
                                                         cGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATcgcg                                                                                       <  2:2956942/148‑1 (MQ=255)
                                                         cGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATcgcg                                                                                       <  2:1849150/148‑1 (MQ=255)
                                                                ccGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTa                                                                                <  1:3369453/148‑1 (MQ=255)
                                                                ccGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTa                                                                                <  2:1302553/148‑1 (MQ=255)
                                                                   cAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATc                                                                             <  2:283219/148‑1 (MQ=255)
                                                                    aaaTCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATct                                                                            >  1:1548560/1‑148 (MQ=255)
                                                                     aaTCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctg                                                                           <  2:2214923/148‑1 (MQ=255)
                                                                     aaTCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctg                                                                           >  1:3404877/1‑148 (MQ=255)
                                                                     aaTCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctg                                                                           <  1:2138405/148‑1 (MQ=255)
                                                                     aaTCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctg                                                                           <  2:2764182/148‑1 (MQ=255)
                                                                       tCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctgct                                                                         >  1:2594033/1‑148 (MQ=255)
                                                                       tCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctgct                                                                         >  2:3658825/1‑148 (MQ=255)
                                                                        cTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGACTTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctgctg                                                                        >  1:1781454/1‑148 (MQ=255)
                                                                        cTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATctgctg                                                                        >  1:3675305/1‑148 (MQ=255)
                                                                         tttAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGa                                                                       >  1:792601/1‑148 (MQ=255)
                                                                          ttAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGaa                                                                      <  1:633887/148‑1 (MQ=255)
                                                                          ttAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGaa                                                                      >  1:1858121/1‑148 (MQ=255)
                                                                            aaCGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAg                                                                    <  1:2326135/148‑1 (MQ=255)
                                                                              cGATGAAGTGGTGGTGCCAGCGTTCCACAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTg                                                                  >  1:2222965/1‑148 (MQ=255)
                                                                                aTGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGAt                                                                >  1:3507133/1‑148 (MQ=255)
                                                                                aTGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGAt                                                                >  1:1823881/1‑148 (MQ=255)
                                                                                   aaGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGa                                                             >  1:1969249/1‑148 (MQ=255)
                                                                                    aGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGaa                                                            <  2:3727264/148‑1 (MQ=255)
                                                                                      tggtggtgCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAAt                                                          <  2:1239584/148‑1 (MQ=255)
                                                                                      tggtggtgCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAAt                                                          <  1:1593050/148‑1 (MQ=255)
                                                                                      tggtggtgCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAAt                                                          >  2:47158/1‑148 (MQ=255)
                                                                                          ggtgCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATctgc                                                      <  1:3747446/148‑1 (MQ=255)
                                                                                           gtgCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATctgcc                                                     <  1:652230/148‑1 (MQ=255)
                                                                                                    ttCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGaaa                                            >  2:1801197/1‑148 (MQ=255)
                                                                                                       cTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATg                                         >  2:3101185/1‑148 (MQ=255)
                                                                                                        tCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGa                                        <  2:3777007/148‑1 (MQ=255)
                                                                                                        tCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGa                                        >  1:890604/1‑148 (MQ=255)
                                                                                                        tCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGa                                        <  2:1515952/148‑1 (MQ=255)
                                                                                                         cAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGaa                                       <  2:525329/148‑1 (MQ=255)
                                                                                                         cAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGaa                                       >  1:1151759/1‑148 (MQ=255)
                                                                                                             cGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCa                                   <  1:412272/148‑1 (MQ=255)
                                                                                                             cGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCa                                   <  1:3315215/148‑1 (MQ=255)
                                                                                                             cGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCa                                   <  1:3315440/148‑1 (MQ=255)
                                                                                                               tGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAAt                                 <  1:1995704/148‑1 (MQ=255)
                                                                                                                   gTTTCGCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAAtgctg                             <  1:1294379/148‑1 (MQ=255)
                                                                                                                        gCTGGAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTg                        <  1:835935/148‑1 (MQ=255)
                                                                                                                         cTGGAAGAGGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATGATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGa                       <  2:1982140/148‑1 (MQ=255)
                                                                                                                           ggAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACt                     >  2:2573544/1‑148 (MQ=255)
                                                                                                                           ggAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACt                     >  2:748338/1‑148 (MQ=255)
                                                                                                                           ggAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACt                     >  2:1335372/1‑148 (MQ=255)
                                                                                                                           ggAAGATGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACt                     >  2:1552365/1‑148 (MQ=255)
                                                                                                                                aTGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGa                >  1:577578/1‑148 (MQ=255)
                                                                                                                                 tGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGaa               >  2:112491/1‑148 (MQ=255)
                                                                                                                                 tGCGCGCGACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGaa               <  2:577578/148‑1 (MQ=255)
                                                                                                                                    gcgcgACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGt            <  2:698293/148‑1 (MQ=255)
                                                                                                                                     cgcgACTATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGtt           <  2:633276/148‑1 (MQ=255)
                                                                                                                                           tATTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGTTTACTgg     >  2:3275578/1‑148 (MQ=255)
                                                                                                                                            aTTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGTTTACTGGa    >  2:3581428/1‑148 (MQ=255)
                                                                                                                                            aTTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGTTTACTGGa    >  2:3581391/1‑148 (MQ=255)
                                                                                                                                            aTTCCTTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGTTTACTGGa    >  1:2900972/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTGGAAGATGCGCGCGACTATTCCGTAGTGGGCTGTGTGGAATTATCTATTCCCGGCAGAACCTACGGCTTGCATGACATCGCGATGTTTAATCTGCTGAAAGTGATGGAAATCTGCCTGCATGAAAATGAAGGCAATGCTGCGCTGACTTATGAAGGTTTACTGGAGC  >  CP009273/4135023‑4135312

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: