Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 788,902 G→T 94.5% 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,9020GT94.5% 793.7 / 35.5 254R115S (CGT→AGT) modFfused molybdate transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base G (8/6);  new base T (109/131);  total (117/137)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 256 total aligned reads displayed.

CGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGT  >  CP009273/788755‑789043
                                                                                                                                                   |                                                                                                                                             
cGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAActcg                                                                                                                                             <  2:230301/150‑1 (MQ=255)
    cGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGc                                                                                                                                         >  1:3509792/1‑150 (MQ=255)
    cGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGc                                                                                                                                         <  2:1455999/150‑1 (MQ=255)
      gCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAGAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCAt                                                                                                                                       <  1:2484842/150‑1 (MQ=255)
       cTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATc                                                                                                                                      >  2:4326291/1‑150 (MQ=255)
         catcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCt                                                                                                                                    >  1:1273431/1‑150 (MQ=255)
            caAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTa                                                                                                                                 <  1:2064713/150‑1 (MQ=255)
               gAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACtt                                                                                                                              <  1:586281/150‑1 (MQ=255)
               gAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACtt                                                                                                                              <  2:469708/150‑1 (MQ=255)
                    aacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATcc                                                                                                                         <  2:3975637/150‑1 (MQ=255)
                     acaaGTCAGGCTCCGAAATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCt                                                                                                                        <  1:434740/150‑1 (MQ=255)
                      caaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTg                                                                                                                       <  2:2215676/150‑1 (MQ=255)
                       aaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGa                                                                                                                      >  1:3971427/1‑150 (MQ=255)
                        aGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGaa                                                                                                                     >  1:4072644/1‑150 (MQ=255)
                         gTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAAt                                                                                                                    <  2:2010913/150‑1 (MQ=255)
                            aGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAt                                                                                                                 <  2:3465429/150‑1 (MQ=255)
                            aGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAt                                                                                                                 >  1:2790238/1‑150 (MQ=255)
                              gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                               <  1:4609171/150‑1 (MQ=255)
                              gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                               <  2:3630170/150‑1 (MQ=255)
                               cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtctc                                                                                                              <  1:3712205/150‑1 (MQ=255)
                                      aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                       <  1:3067314/150‑1 (MQ=255)
                                      aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                       >  1:209108/1‑150 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                    >  2:2235164/1‑150 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                    >  2:1866528/1‑150 (MQ=255)
                                            gcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGt                                                                                                 >  2:790338/1‑150 (MQ=255)
                                             ccTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTc                                                                                                >  2:3125026/1‑150 (MQ=255)
                                                gACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCgg                                                                                             >  1:2627020/1‑150 (MQ=255)
                                                 acacAGCCGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgt                                                                                            <  1:1563029/150‑1 (MQ=255)
                                                  cacaGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtg                                                                                           >  2:1257225/1‑150 (MQ=255)
                                                  cacaGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtg                                                                                           <  1:84460/150‑1 (MQ=255)
                                                  cacaGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtg                                                                                           <  1:683565/150‑1 (MQ=255)
                                                     agcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCa                                                                                        >  2:3495571/1‑150 (MQ=255)
                                                     agcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCa                                                                                        <  1:2828953/150‑1 (MQ=255)
                                                       cagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATc                                                                                      >  2:4104342/1‑150 (MQ=255)
                                                        agGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCt                                                                                     >  2:4437738/1‑150 (MQ=255)
                                                         gggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                    >  2:1118211/1‑150 (MQ=255)
                                                              tttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCa                                                                               >  1:1420392/1‑150 (MQ=255)
                                                                     tctcGCCAGTGGAAAGGTATTTACAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGGGTCATCTTCGCCAGGGCCGa                                                                        <  2:598455/150‑1 (MQ=255)
                                                                     tctcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGa                                                                        <  1:3488199/150‑1 (MQ=255)
                                                                       tcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGc                                                                      <  1:2920208/150‑1 (MQ=255)
                                                                       tcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGc                                                                      >  2:2176902/1‑150 (MQ=255)
                                                                        cGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCa                                                                     <  2:2895249/150‑1 (MQ=255)
                                                                        cGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCa                                                                     >  2:1765706/1‑150 (MQ=255)
                                                                         gCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCat                                                                    <  2:2472756/150‑1 (MQ=255)
                                                                         gCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCat                                                                    <  1:3076866/150‑1 (MQ=255)
                                                                           cAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat                                                                  >  1:1716696/1‑150 (MQ=255)
                                                                            aGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATc                                                                 <  2:2720999/150‑1 (MQ=255)
                                                                               ggAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGGGATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgt                                                              >  1:1356746/1‑150 (MQ=255)
                                                                                gAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtg                                                             <  2:2832190/150‑1 (MQ=255)
                                                                                 aaaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtgt                                                            <  1:3376398/150‑1 (MQ=255)
                                                                                  aaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGtt                                                           >  1:1873856/1‑150 (MQ=255)
                                                                                    ggTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTAt                                                         <  2:1211125/150‑1 (MQ=255)
                                                                                      tATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTc                                                       <  1:2122773/150‑1 (MQ=255)
                                                                                       aTTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTcc                                                      >  1:3204689/1‑150 (MQ=255)
                                                                                        tttAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCg                                                     <  1:2620318/150‑1 (MQ=255)
                                                                                        tttAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCg                                                     <  2:1055272/150‑1 (MQ=255)
                                                                                        tttAAAGCGTCGGTCGAGGAGGGCGGAAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCg                                                     >  2:1751575/1‑150 (MQ=255)
                                                                                          tAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCt                                                   <  2:3285232/150‑1 (MQ=255)
                                                                                          tAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCt                                                   <  2:3515428/150‑1 (MQ=255)
                                                                                          tAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCt                                                   <  2:2008040/150‑1 (MQ=255)
                                                                                              gCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCAACCTTTACTGCAACCTGAATGATCTAAGCCGTAGTGCGTCCGTTGTCATCTTCGCCAGGGCCGAGCATATCGTTGTTATTCCGCTGCCa                                               >  1:4034652/1‑150 (MQ=255)
                                                                                               cGTCGGTCGAGGGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAt                                              >  1:1062797/1‑150 (MQ=255)
                                                                                                gtcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAtt                                             <  1:3572086/150‑1 (MQ=255)
                                                                                                gtcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAtt                                             >  2:4498630/1‑150 (MQ=255)
                                                                                                 tcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtc                                            <  2:784684/150‑1 (MQ=255)
                                                                                                 tcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtc                                            >  1:4292710/1‑150 (MQ=255)
                                                                                                  cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcg                                           <  1:738681/150‑1 (MQ=255)
                                                                                                    gtcgAGGAGGGCGGTTATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCAATtcgtc                                         <  2:4148069/150‑1 (MQ=255)
                                                                                                      cgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGc                                       >  2:913228/1‑150 (MQ=255)
                                                                                                      cgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGc                                       >  1:881859/1‑150 (MQ=255)
                                                                                                       gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                      >  2:3838430/1‑150 (MQ=255)
                                                                                                       gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                      <  1:3599994/150‑1 (MQ=255)
                                                                                                       gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                      <  1:3310014/150‑1 (MQ=255)
                                                                                                        aggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTg                                     <  2:4340396/150‑1 (MQ=255)
                                                                                                         ggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGa                                    <  1:3168072/150‑1 (MQ=255)
                                                                                                         ggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGa                                    >  2:1524978/1‑150 (MQ=255)
                                                                                                                   aaTACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTg                          >  1:4587633/1‑150 (MQ=255)
                                                                                                                     tACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCCGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCa                        <  1:2416698/150‑1 (MQ=255)
                                                                                                                         gAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTg                    <  1:3447657/150‑1 (MQ=255)
                                                                                                                          aaCTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGc                   >  2:2330721/1‑150 (MQ=255)
                                                                                                                           aCTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCt                  <  2:4335446/150‑1 (MQ=255)
                                                                                                                           aCTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCt                  <  2:1836201/150‑1 (MQ=255)
                                                                                                                            ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTc                 <  1:2353416/150‑1 (MQ=255)
                                                                                                                               ctgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGaa              <  2:4353908/150‑1 (MQ=255)
                                                                                                                                tgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAg             <  2:682647/150‑1 (MQ=255)
                                                                                                                                tgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAg             >  1:380330/1‑150 (MQ=255)
                                                                                                                                tgcGCCAGTTGCATGCAACGCGGTGCATCCTTTACCTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAg             <  1:886474/150‑1 (MQ=255)
                                                                                                                                 gcgcCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAgg            >  2:3065832/1‑150 (MQ=255)
                                                                                                                                  cgcCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGCGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGga           <  2:788401/150‑1 (MQ=255)
                                                                                                                                     cAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagag        <  2:1273431/150‑1 (MQ=255)
                                                                                                                                      aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagaga       >  1:1065598/1‑150 (MQ=255)
                                                                                                                                      aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagaga       <  2:1397125/150‑1 (MQ=255)
                                                                                                                                      aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagaga       >  1:4264741/1‑150 (MQ=255)
                                                                                                                                       gTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc      <  1:3235666/150‑1 (MQ=255)
                                                                                                                                       gTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc      <  1:2678104/150‑1 (MQ=255)
                                                                                                                                        ttGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACg     >  1:4529468/1‑150 (MQ=255)
                                                                                                                                        ttGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACg     >  2:766844/1‑150 (MQ=255)
                                                                                                                                         tgcatgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGa    >  1:907900/1‑150 (MQ=255)
                                                                                                                                           catgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGt  >  2:2544377/1‑150 (MQ=255)
                                                                                                                                           catgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGt  >  1:3417283/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGT  >  CP009273/788755‑789043

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: