Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 792,957 C→A 11.6% V316V (GTC→GTA modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273792,9570CA11.6% 1003.4 / 68.1 280V316V (GTC→GTAmodCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (113/134);  new base A (17/16);  total (130/150)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

CTGCAAGCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTG  >  CP009273/792808‑793105
                                                                                                                                                     |                                                                                                                                                    
   cAAGCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTaggc                                                                                                                                                   <  1:1613183/150‑1 (MQ=255)
        tGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTaa                                                                                                                                              <  2:2635421/150‑1 (MQ=255)
          cgcTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTaaaa                                                                                                                                            >  2:2166173/1‑150 (MQ=255)
           gcTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAAc                                                                                                                                           >  1:3080518/1‑150 (MQ=255)
           gcTACGTATCCGCAGTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAAc                                                                                                                                           <  2:3164080/150‑1 (MQ=255)
                     cGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGgcg                                                                                                                                 <  2:3555409/150‑1 (MQ=255)
                      gCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGgcgc                                                                                                                                >  1:2487239/1‑150 (MQ=255)
                        aTTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTTTTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGTAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGt                                                                                                                              <  1:1749406/150‑1 (MQ=255)
                        aTTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGt                                                                                                                              <  2:761969/150‑1 (MQ=255)
                        aTTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGt                                                                                                                              >  2:3758900/1‑150 (MQ=255)
                        aTTCAGGCTTCCGATGTTTCTCAGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGt                                                                                                                              >  2:878857/1‑150 (MQ=255)
                         ttCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTa                                                                                                                             >  2:3307259/1‑150 (MQ=255)
                         ttCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTa                                                                                                                             >  2:1865821/1‑150 (MQ=255)
                          tCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTAt                                                                                                                            <  2:3956918/150‑1 (MQ=255)
                           cAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATACGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGTAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATc                                                                                                                           >  2:1618187/1‑150 (MQ=255)
                            aGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCa                                                                                                                          <  2:3710811/150‑1 (MQ=255)
                                ttCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGccc                                                                                                                      <  2:2523955/150‑1 (MQ=255)
                                 tCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCg                                                                                                                     >  1:2087061/1‑150 (MQ=255)
                                 tCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCg                                                                                                                     >  2:1897702/1‑150 (MQ=255)
                                 tCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGACAGTTGGAAGTGGAACTGGAAGTCGGAGGTAAAACGATGTGGGCGCGTATCAGCCCg                                                                                                                     >  2:638871/1‑150 (MQ=255)
                                   cGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTg                                                                                                                   <  1:2343728/150‑1 (MQ=255)
                                   cGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTg                                                                                                                   <  2:3345662/150‑1 (MQ=255)
                                    gATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTgg                                                                                                                  >  1:1567144/1‑150 (MQ=255)
                                       gTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGcc                                                                                                               <  2:1980438/150‑1 (MQ=255)
                                          tctcTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAgg                                                                                                            <  1:2281430/150‑1 (MQ=255)
                                            tctGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGa                                                                                                          <  2:1572578/150‑1 (MQ=255)
                                              tGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATg                                                                                                        >  2:3429752/1‑150 (MQ=255)
                                              tGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATg                                                                                                        >  2:2754496/1‑150 (MQ=255)
                                               ggTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGa                                                                                                       >  1:1366782/1‑150 (MQ=255)
                                               ggTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGa                                                                                                       <  2:1137926/150‑1 (MQ=255)
                                                 ttttACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAAc                                                                                                     >  1:3698710/1‑150 (MQ=255)
                                                 ttttACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAAc                                                                                                     <  2:2354423/150‑1 (MQ=255)
                                                     acacCCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGc                                                                                                 <  1:4587618/150‑1 (MQ=255)
                                                      cAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCg                                                                                                <  2:4466019/150‑1 (MQ=255)
                                                        aCCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGAt                                                                                              >  2:1117368/1‑150 (MQ=255)
                                                         ccgccgCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATc                                                                                             <  2:2442961/150‑1 (MQ=255)
                                                         ccgccgCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATc                                                                                             <  1:1556939/150‑1 (MQ=255)
                                                           gccgcAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCaa                                                                                           <  1:2678266/150‑1 (MQ=255)
                                                                agcaAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTg                                                                                      <  1:2939624/150‑1 (MQ=255)
                                                                   aaaCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGcc                                                                                   <  2:2667849/150‑1 (MQ=255)
                                                                   aaaCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGcc                                                                                   <  1:285405/150‑1 (MQ=255)
                                                                    aaCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCt                                                                                  <  1:1510085/150‑1 (MQ=255)
                                                                     aCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCtg                                                                                 <  1:4614407/150‑1 (MQ=255)
                                                                     aCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCtg                                                                                 <  1:1958630/150‑1 (MQ=255)
                                                                      ccAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCtgt                                                                                <  1:653869/150‑1 (MQ=255)
                                                                       cAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCtgtg                                                                               <  2:575490/150‑1 (MQ=255)
                                                                       cAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCtgtg                                                                               <  1:49517/150‑1 (MQ=255)
                                                                        aGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTgg                                                                              <  2:1048432/150‑1 (MQ=255)
                                                                          cATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCt                                                                            <  1:1985749/150‑1 (MQ=255)
                                                                          cATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCt                                                                            <  2:1357534/150‑1 (MQ=255)
                                                                           aTTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTg                                                                           >  1:2876043/1‑150 (MQ=255)
                                                                            ttCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGt                                                                          >  2:1771109/1‑150 (MQ=255)
                                                                             tCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTa                                                                         <  1:3268486/150‑1 (MQ=255)
                                                                               gTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAcg                                                                       >  1:101870/1‑150 (MQ=255)
                                                                               gTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAcg                                                                       <  1:3922408/150‑1 (MQ=255)
                                                                               gTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAcg                                                                       <  1:4045557/150‑1 (MQ=255)
                                                                               gTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAcg                                                                       <  1:3771996/150‑1 (MQ=255)
                                                                                     tATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCaca                                                                 >  1:3481738/1‑148 (MQ=255)
                                                                                      aTTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAAt                                                                >  1:835245/1‑150 (MQ=255)
                                                                                        tGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTa                                                              <  2:1499424/150‑1 (MQ=255)
                                                                                           gggCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTaaaa                                                           >  2:259797/1‑150 (MQ=255)
                                                                                             gCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgt                                                         >  2:3296137/1‑150 (MQ=255)
                                                                                                 aaGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtgt                                                     >  1:2350896/1‑150 (MQ=255)
                                                                                                     tgttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGAt                                                 <  2:4227184/150‑1 (MQ=255)
                                                                                                     tgttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGAt                                                 <  2:3819942/150‑1 (MQ=255)
                                                                                                        tAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAAc                                              >  1:1294541/1‑150 (MQ=255)
                                                                                                             gTTATGACGACAACGGCCAGGTGGAAGTGGNACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCt                                         >  1:4284124/1‑150 (MQ=255)
                                                                                                              ttATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTg                                        >  2:671774/1‑150 (MQ=255)
                                                                                                              ttATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTg                                        >  1:3468192/1‑150 (MQ=255)
                                                                                                                aTGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGAt                                      <  1:992298/150‑1 (MQ=255)
                                                                                                                  gacgacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                    <  1:2261364/150‑1 (MQ=255)
                                                                                                                   acgacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                   >  1:2765309/1‑150 (MQ=255)
                                                                                                                   acgacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                   <  2:1157072/150‑1 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 >  1:3939807/1‑150 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 >  1:4141082/1‑150 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 <  2:3365859/150‑1 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 >  1:1160345/1‑150 (MQ=255)
                                                                                                                      acAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                                <  2:802953/150‑1 (MQ=255)
                                                                                                                       cAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCa                               <  2:1208716/150‑1 (MQ=255)
                                                                                                                        aaCGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAg                              >  1:3334799/1‑150 (MQ=255)
                                                                                                                        aaCGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAg                              >  1:3549940/1‑150 (MQ=255)
                                                                                                                         aCGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAgg                             <  1:3758900/150‑1 (MQ=255)
                                                                                                                           ggCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTg                           >  1:2783791/1‑150 (MQ=255)
                                                                                                                           ggCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTg                           >  1:337630/1‑150 (MQ=255)
                                                                                                                             ccAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGc                         >  1:2222081/1‑150 (MQ=255)
                                                                                                                               aGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                       >  2:2305962/1‑150 (MQ=255)
                                                                                                                                  tggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATaa                    >  1:3641913/1‑150 (MQ=255)
                                                                                                                                  tggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATaa                    >  1:1095919/1‑150 (MQ=255)
                                                                                                                                   ggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATaaa                   >  1:4504180/1‑150 (MQ=255)
                                                                                                                                     aagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATg                 <  2:3548166/150‑1 (MQ=255)
                                                                                                                                       gtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGaa               <  1:1202423/150‑1 (MQ=255)
                                                                                                                                        tggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAAc              >  1:2531611/1‑150 (MQ=255)
                                                                                                                                         ggaaCTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACt             >  2:1571485/1‑150 (MQ=255)
                                                                                                                                            aCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTggg          <  1:1865821/150‑1 (MQ=255)
                                                                                                                                              tGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCa        <  1:3208118/150‑1 (MQ=255)
                                                                                                                                              tGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCa        <  1:885340/150‑1 (MQ=255)
                                                                                                                                               ggAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCaa       <  1:1992544/150‑1 (MQ=255)
                                                                                                                                               ggAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCaa       <  1:3721480/150‑1 (MQ=255)
                                                                                                                                               ggAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCaa       >  2:3231270/1‑150 (MQ=255)
                                                                                                                                                    tCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATACATGAACTGGGCAATGCTg  >  2:4327387/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCAAGCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTG  >  CP009273/792808‑793105

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: