Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 788,902 G→T 92.9% R115S (CGT→AGT)  modF ← fused molybdate transporter subunits of ABC superfamily: ATP‑binding components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273788,9020GT92.9% 736.2 / 40.6 226R115S (CGT→AGT) modFfused molybdate transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base G (6/10);  new base T (112/98);  total (118/108)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 226 total aligned reads displayed.

GAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGT  >  CP009273/788756‑789048
                                                                                                                                                  |                                                                                                                                                  
 aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                >  2:1326371/1‑150 (MQ=255)
 aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                <  2:2597142/150‑1 (MQ=255)
        catcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCt                                                                                                                                         <  1:2338018/150‑1 (MQ=255)
         atcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCtt                                                                                                                                        <  2:3992518/150‑1 (MQ=255)
         atcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCtt                                                                                                                                        <  2:3552731/150‑1 (MQ=255)
          tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCttt                                                                                                                                       <  2:1414772/150‑1 (MQ=255)
          tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCttt                                                                                                                                       >  2:4568556/1‑150 (MQ=255)
            aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTAc                                                                                                                                     >  1:3271610/1‑150 (MQ=255)
             aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACt                                                                                                                                    >  2:2545405/1‑150 (MQ=255)
             aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTATCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGGTGCATGCCACTCGGTGCATCCTTTCCt                                                                                                                                    >  2:3702383/1‑150 (MQ=255)
              gAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACtt                                                                                                                                   <  2:357222/150‑1 (MQ=255)
                 tCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCAt                                                                                                                                <  1:1113465/150‑1 (MQ=255)
                 tCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCAt                                                                                                                                <  1:3322825/150‑1 (MQ=255)
                   aacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATcc                                                                                                                              >  1:4926445/1‑150 (MQ=255)
                      aaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGa                                                                                                                           <  2:114460/150‑1 (MQ=255)
                            ggCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATACTGAATGAtc                                                                                                                     >  2:4875203/1‑150 (MQ=255)
                             gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                                    <  2:1544598/150‑1 (MQ=255)
                             gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGGATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                                    >  1:2245530/1‑150 (MQ=255)
                                ccGACATCAGCGCCTGACACAGCAGGGTTTGTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAg                                                                                                                 <  2:5110820/150‑1 (MQ=255)
                                 cGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGc                                                                                                                >  1:1822976/1‑150 (MQ=255)
                                  gACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGcc                                                                                                               <  1:923004/150‑1 (MQ=255)
                                   aCATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCg                                                                                                              >  1:3915572/1‑150 (MQ=255)
                                    cATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGt                                                                                                             <  1:1535833/150‑1 (MQ=255)
                                     aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            >  2:3535599/1‑150 (MQ=255)
                                     aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            >  1:2449109/1‑150 (MQ=255)
                                     aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            >  2:4851192/1‑150 (MQ=255)
                                     aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            <  1:1478492/150‑1 (MQ=255)
                                     aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            <  1:2845967/150‑1 (MQ=255)
                                     aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            <  1:2121314/150‑1 (MQ=255)
                                      tCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGACAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                                                                                           >  2:660385/1‑150 (MQ=255)
                                       cAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGt                                                                                                          <  1:553126/150‑1 (MQ=255)
                                        aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                         >  1:775046/1‑150 (MQ=255)
                                        aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                         >  1:1346193/1‑150 (MQ=255)
                                          cgcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCg                                                                                                       >  1:5062439/1‑150 (MQ=255)
                                              tGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCg                                                                                                   <  2:4532339/150‑1 (MQ=255)
                                                      cagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATc                                                                                           >  2:1526073/1‑150 (MQ=255)
                                                       agGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCt                                                                                          >  1:3486810/1‑150 (MQ=255)
                                                        gggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                         <  2:4634398/150‑1 (MQ=255)
                                                            ttttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGcc                                                                                     <  2:296662/150‑1 (MQ=255)
                                                            ttttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGcc                                                                                     <  1:1876174/150‑1 (MQ=255)
                                                             tttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCa                                                                                    <  2:3804237/150‑1 (MQ=255)
                                                                      tcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGc                                                                           >  2:4010070/1‑150 (MQ=255)
                                                                        gCCAGTGGAAAGGTATTTACAGCGTCGGTCGAGGCGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCat                                                                         <  2:3881338/150‑1 (MQ=255)
                                                                         ccAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCata                                                                        >  1:4455055/1‑150 (MQ=255)
                                                                         ccAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCata                                                                        >  1:3669652/1‑150 (MQ=255)
                                                                          cAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat                                                                       <  2:1408697/150‑1 (MQ=255)
                                                                          cAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat                                                                       >  1:1934320/1‑150 (MQ=255)
                                                                          cAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat                                                                       <  1:3070419/150‑1 (MQ=255)
                                                                            gTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCg                                                                     <  1:2180716/150‑1 (MQ=255)
                                                                                aaaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtgt                                                                 >  2:2794059/1‑150 (MQ=255)
                                                                                  aGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTa                                                               >  1:1689521/1‑150 (MQ=255)
                                                                                      aTTTAAAGCGTCGGTCGAGGAGGGCGGTAATGCCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCAGATCGGTGTTCTTcc                                                           >  2:3304979/1‑150 (MQ=255)
                                                                                      aTTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTcc                                                           >  2:2321494/1‑150 (MQ=255)
                                                                                      aTTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTcc                                                           >  2:4926344/1‑150 (MQ=255)
                                                                                          aaaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTg                                                       >  2:1352489/1‑150 (MQ=255)
                                                                                          aaaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTg                                                       >  2:4868076/1‑150 (MQ=255)
                                                                                           aaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGc                                                      >  2:703603/1‑150 (MQ=255)
                                                                                              cGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAt                                                   >  2:5191619/1‑150 (MQ=255)
                                                                                               gtcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAtt                                                  <  1:2310030/150‑1 (MQ=255)
                                                                                                tcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtc                                                 <  1:4845857/150‑1 (MQ=255)
                                                                                                 cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcg                                                <  2:1191646/150‑1 (MQ=255)
                                                                                                 cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcg                                                <  1:2470956/150‑1 (MQ=255)
                                                                                                  ggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                               >  2:4336718/1‑150 (MQ=255)
                                                                                                  ggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                               <  1:2677017/150‑1 (MQ=255)
                                                                                                  ggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                               <  2:1934320/150‑1 (MQ=255)
                                                                                                   gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtc                                              <  2:3197541/150‑1 (MQ=255)
                                                                                                   gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtc                                              <  1:2519223/150‑1 (MQ=255)
                                                                                                     cgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGc                                            <  2:2069429/150‑1 (MQ=255)
                                                                                                     cgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGc                                            >  2:2541353/1‑150 (MQ=255)
                                                                                                     cgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGc                                            >  2:1299121/1‑150 (MQ=255)
                                                                                                      gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                           >  1:1334656/1‑150 (MQ=255)
                                                                                                      gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                           >  1:3364873/1‑150 (MQ=255)
                                                                                                      gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                           >  1:3771565/1‑150 (MQ=255)
                                                                                                      gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                           >  1:4418937/1‑150 (MQ=255)
                                                                                                       aggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTg                                          >  1:3196415/1‑150 (MQ=255)
                                                                                                            ggCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAg                                     <  1:1899790/150‑1 (MQ=255)
                                                                                                              cGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGtt                                   >  2:5124307/1‑150 (MQ=255)
                                                                                                               ggTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGttt                                  >  2:51448/1‑150 (MQ=255)
                                                                                                                   aTACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGc                              <  1:2420890/150‑1 (MQ=255)
                                                                                                                   aTACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGc                              >  2:1760953/1‑150 (MQ=255)
                                                                                                                     aCCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCaa                            >  1:832433/1‑150 (MQ=255)
                                                                                                                      ccGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAAt                           >  1:1011058/1‑150 (MQ=255)
                                                                                                                           ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTc                      >  1:116220/1‑150 (MQ=255)
                                                                                                                             gctgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGa                    <  2:2141669/150‑1 (MQ=255)
                                                                                                                               tgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAg                  <  1:488577/150‑1 (MQ=255)
                                                                                                                               tgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAg                  <  1:4471588/150‑1 (MQ=255)
                                                                                                                                gcgcCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAgg                 >  2:1270270/1‑150 (MQ=255)
                                                                                                                                  gcCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgag               <  2:1815570/150‑1 (MQ=255)
                                                                                                                                  gcCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgag               <  2:1815547/150‑1 (MQ=255)
                                                                                                                                  gcCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgag               <  1:790063/150‑1 (MQ=255)
                                                                                                                                      gTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc           >  1:1986075/1‑150 (MQ=255)
                                                                                                                                      gTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc           <  1:4336718/150‑1 (MQ=255)
                                                                                                                                        tgcatgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGa         <  1:4668413/150‑1 (MQ=255)
                                                                                                                                        tgcatgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGa         <  1:5011796/150‑1 (MQ=255)
                                                                                                                                          catgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGt       <  1:4382289/150‑1 (MQ=255)
                                                                                                                                            tgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGa     >  2:808093/1‑150 (MQ=255)
                                                                                                                                            tgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGa     >  2:1091348/1‑150 (MQ=255)
                                                                                                                                               aaCTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATCTCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAtgt  >  1:2564030/1‑150 (MQ=255)
                                                                                                                                               aaCTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAtgt  <  2:225832/150‑1 (MQ=255)
                                                                                                                                               aaCTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAtgt  >  1:4039270/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGT  >  CP009273/788756‑789048

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: