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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273792,9570CA 5.0% 1012.9 / 20.6 257V316V (GTC→GTAmodCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (110/134);  new base A (11/2);  total (121/136)
Fisher's exact test for biased strand distribution p-value = 7.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 269 total aligned reads displayed.

CTGCAAGCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGCTG  >  CP009273/792808‑793105
                                                                                                                                                     |                                                                                                                                                    
cTGCAAGCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTc                                                                                                                                                      >  1:2374838/1‑150 (MQ=255)
  gCAAGCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTAATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTcgg                                                                                                                                                    >  2:3278425/1‑150 (MQ=255)
     aGCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTcggcgg                                                                                                                                                 >  1:1005120/1‑150 (MQ=255)
      gCTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGt                                                                                                                                                >  1:45149/1‑150 (MQ=255)
       cTGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTa                                                                                                                                               <  1:3415847/150‑1 (MQ=255)
        tGCGCTACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTaa                                                                                                                                              >  1:3917218/1‑150 (MQ=255)
             tACGTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGACGTGGAACTGGAAGTCGGCGGTAAAACGc                                                                                                                                         >  1:2116865/1‑150 (MQ=255)
                gTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGACGTGGAACTGGAAGTCGGCGGTAAAACGCtgt                                                                                                                                      >  1:1630376/1‑150 (MQ=255)
                gTATCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCtgt                                                                                                                                      <  1:604371/150‑1 (MQ=255)
                   tCCGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACACCCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTggg                                                                                                                                   <  2:3986118/150‑1 (MQ=255)
                    ccGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGgc                                                                                                                                  <  2:432110/150‑1 (MQ=255)
                    ccGCATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGgc                                                                                                                                  >  1:3192333/1‑150 (MQ=255)
                       cATTCAGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGcgcg                                                                                                                               >  2:1618205/1‑150 (MQ=255)
                         ttCATGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTa                                                                                                                             >  1:2753806/1‑150 (MQ=255)
                            aGGCTTCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCa                                                                                                                          >  2:2592741/1‑150 (MQ=255)
                                ttCCGATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGAGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGccc                                                                                                                      >  2:3367153/1‑150 (MQ=255)
                                    gATGTTTCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTgg                                                                                                                  <  2:2327422/150‑1 (MQ=255)
                                         ttCTCTGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAg                                                                                                             <  1:2266020/150‑1 (MQ=255)
                                           ctctGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGTGCAAAATTTGTTAATTGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAggg                                                                                                           <  1:3000128/150‑1 (MQ=255)
                                            tctGGTTTTACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGa                                                                                                          <  2:683128/150‑1 (MQ=255)
                                                 ttttACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCAGTGGGCCAGGGATGAAc                                                                                                     <  1:3269778/150‑1 (MQ=255)
                                                  tttACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACt                                                                                                    <  2:3893891/150‑1 (MQ=255)
                                                   ttACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTg                                                                                                   <  1:4048828/150‑1 (MQ=255)
                                                    tACAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTgg                                                                                                  >  2:3017419/1‑150 (MQ=255)
                                                      cAACCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCCAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCg                                                                                                >  2:2356200/1‑150 (MQ=255)
                                                       aaCCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGa                                                                                               <  2:1743566/150‑1 (MQ=255)
                                                        aCCGCCGCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGAt                                                                                              <  2:3193352/150‑1 (MQ=255)
                                                         ccgccgCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATc                                                                                             >  1:352722/1‑150 (MQ=255)
                                                         ccgccgCAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATc                                                                                             >  1:760719/1‑150 (MQ=255)
                                                             cgcAGCAAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAAc                                                                                         <  2:661642/150‑1 (MQ=255)
                                                              gcagcaAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAAcc                                                                                        <  2:1962731/150‑1 (MQ=255)
                                                              gcagcaAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAAcc                                                                                        <  2:380019/150‑1 (MQ=255)
                                                               cagcaAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCt                                                                                       <  1:3047264/150‑1 (MQ=255)
                                                               cagcaAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCt                                                                                       >  1:1671236/1‑150 (MQ=255)
                                                                agcaAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGCCGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTg                                                                                      <  1:992167/150‑1 (MQ=255)
                                                                  caAACCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGc                                                                                    <  2:2314423/150‑1 (MQ=255)
                                                                   aaaCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATATTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGACGATCAAACCTGGcc                                                                                   >  1:1391376/1‑150 (MQ=255)
                                                                   aaaCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGcc                                                                                   >  2:2493214/1‑150 (MQ=255)
                                                                    aaCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCt                                                                                  >  1:2074938/1‑150 (MQ=255)
                                                                     aCCAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCtg                                                                                 <  1:2686799/150‑1 (MQ=255)
                                                                       cAGCATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCtgtg                                                                               <  1:1495895/150‑1 (MQ=255)
                                                                          cATTCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCt                                                                            >  2:3081250/1‑150 (MQ=255)
                                                                            ttCGTAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGt                                                                          <  2:1498776/150‑1 (MQ=255)
                                                                              cGTAACGTATTGCGGGCAAATGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAc                                                                        <  1:3570880/150‑1 (MQ=255)
                                                                                tAACGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAcgc                                                                      >  2:3704176/1‑150 (MQ=255)
                                                                                 aaCGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAcgcg                                                                     <  2:277117/150‑1 (MQ=255)
                                                                                 aaCGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTAcgcg                                                                     >  2:971734/1‑150 (MQ=255)
                                                                                   cGTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCa                                                                   <  2:125731/150‑1 (MQ=255)
                                                                                    gTATTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCaa                                                                  >  2:1875978/1‑150 (MQ=255)
                                                                                      aTTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAAt                                                                <  2:2935896/150‑1 (MQ=255)
                                                                                      aTTGCGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAAt                                                                <  2:392725/150‑1 (MQ=255)
                                                                                          cGGGCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTaaa                                                            >  2:587420/1‑150 (MQ=255)
                                                                                             gCAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgt                                                         <  1:2820552/150‑1 (MQ=255)
                                                                                              cAAAAGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtg                                                        >  1:1102105/1‑150 (MQ=255)
                                                                                                aaaGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtg                                                      >  2:3155961/1‑150 (MQ=255)
                                                                                                  aGTTGTTAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGACGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTc                                                    <  1:125335/150‑1 (MQ=255)
                                                                                                   gttgttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCg                                                   <  2:1377278/150‑1 (MQ=255)
                                                                                                    ttgttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCACATTAAAAGTGTGTCGa                                                  >  2:783024/1‑150 (MQ=255)
                                                                                                     tgttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGAt                                                 <  2:1883544/150‑1 (MQ=255)
                                                                                                      gttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATa                                                >  1:1578209/1‑150 (MQ=255)
                                                                                                      gttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATa                                                <  1:1651132/150‑1 (MQ=255)
                                                                                                       ttAATAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATaa                                               <  2:2978783/150‑1 (MQ=255)
                                                                                                         aaTAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAAcc                                             <  2:1189082/150‑1 (MQ=255)
                                                                                                         aaTAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAAcc                                             <  2:4117357/150‑1 (MQ=255)
                                                                                                         aaTAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAAcc                                             <  2:3802566/150‑1 (MQ=255)
                                                                                                          aTAGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCg                                            <  1:3554106/150‑1 (MQ=255)
                                                                                                            aGTTATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGcc                                          <  2:2836849/150‑1 (MQ=255)
                                                                                                               tATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGa                                       >  2:3115768/1‑150 (MQ=255)
                                                                                                               tATGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGa                                       >  1:1468130/1‑150 (MQ=255)
                                                                                                                 tGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGAtt                                     >  2:1321088/1‑150 (MQ=255)
                                                                                                                 tGACGACAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGAtt                                     <  1:2311919/150‑1 (MQ=255)
                                                                                                                  gacgacAACGGCCAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                    >  1:3874822/1‑150 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 >  1:3689567/1‑150 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 >  2:1830135/1‑150 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 >  2:630701/1‑150 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 <  2:2145242/150‑1 (MQ=255)
                                                                                                                     gacAACGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAAt                                 <  1:2836392/150‑1 (MQ=255)
                                                                                                                        aaCGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAg                              >  1:966378/1‑150 (MQ=255)
                                                                                                                        aaCGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAg                              <  1:2033633/150‑1 (MQ=255)
                                                                                                                        aaCGGCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAg                              >  1:2483139/1‑150 (MQ=255)
                                                                                                                           ggCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTg                           >  2:899005/1‑150 (MQ=255)
                                                                                                                            gCCAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTTATTTAATCAGGTgg                          <  1:17596/150‑1 (MQ=255)
                                                                                                                             ccAGGTGGAAGTGGAACTGGAAGTAGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGc                         >  1:1818046/1‑150 (MQ=255)
                                                                                                                              cAGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCt                        <  2:1311229/150‑1 (MQ=255)
                                                                                                                               aGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTCCGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                       >  2:547361/1‑150 (MQ=255)
                                                                                                                               aGGTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                       >  1:938479/1‑150 (MQ=255)
                                                                                                                                ggTGGAAGTGGAACTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtat                      >  2:3579399/1‑150 (MQ=255)
                                                                                                                                 gtggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                     <  1:31771/150‑1 (MQ=255)
                                                                                                                                  tggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATaa                    <  1:2246671/150‑1 (MQ=255)
                                                                                                                                  tggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATaa                    <  2:2274488/150‑1 (MQ=255)
                                                                                                                                   ggaagtggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATaaa                   >  1:3899340/1‑150 (MQ=255)
                                                                                                                                        tggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAAc              <  2:2390908/150‑1 (MQ=255)
                                                                                                                                         ggaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACt             >  1:1950234/1‑150 (MQ=255)
                                                                                                                                         ggaaCTGGAAGTCGGCGGTAAAACGCTGTGATCGCGAAACAGCACCTGGGCCAGTGATGAACTGGCAATAAAACCTGGCCTGTGGCTGTACGCGCAAATAAAAAGTGATTAGAAAACCGCCTGAATAAATCAGGAGTCTATAAAAGAACt             >  2:2063603/1‑150 (MQ=255)
                                                                                                                                          gaaCTGGAAGTCGGCGGTATAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCATACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGAGAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTg            <  2:755064/150‑1 (MQ=255)
                                                                                                                                          gaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTg            <  1:762571/150‑1 (MQ=255)
                                                                                                                                          gaaCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTg            <  1:1411829/150‑1 (MQ=255)
                                                                                                                                            aCTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTggg          >  1:360808/1‑150 (MQ=255)
                                                                                                                                             cTGGAAGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc         >  1:2804384/1‑150 (MQ=255)
                                                                                                                                                  aGTCGGCGGTAAAACGCTGTGGGCGCGTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc    <  1:3791795/150‑1 (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: